OOFEM  2.4
OOFEM.org - Object Oriented Finite Element Solver
oofem::NLTransientTransportProblem Class Reference

This class represents nonlinear transient transport problem. More...

#include <nltransienttransportproblem.h>

+ Inheritance diagram for oofem::NLTransientTransportProblem:
+ Collaboration diagram for oofem::NLTransientTransportProblem:

Public Member Functions

 NLTransientTransportProblem (int i, EngngModel *_master)
 Constructor. More...
 
virtual ~NLTransientTransportProblem ()
 Destructor. More...
 
virtual TimeStepgiveNextStep ()
 Returns next time step (next to current step) of receiver. More...
 
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 IRResultType initializeFrom (InputRecord *ir)
 Initializes receiver according to object description in input reader. More...
 
virtual const char * giveInputRecordName () const
 
virtual const char * giveClassName () const
 Returns class name of the receiver. More...
 
virtual fMode giveFormulation ()
 Indicates type of non linear computation (total or updated formulation). 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 void updateDofUnknownsDictionary (DofManager *dman, TimeStep *tStep)
 Updates necessary values in Dofs unknown dictionaries. More...
 
virtual int giveCurrentNumberOfIterations ()
 Returns number of iterations that was required to reach equilibrium - used for adaptive step length in staggered problem. More...
 
- Public Member Functions inherited from oofem::NonStationaryTransportProblem
 NonStationaryTransportProblem (int i, EngngModel *_master)
 Constructor. More...
 
virtual ~NonStationaryTransportProblem ()
 Destructor. More...
 
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 TimeStepgiveSolutionStepWhenIcApply (bool force=false)
 Returns the solution step when Initial Conditions (IC) apply. More...
 
virtual NumericalMethodgiveNumericalMethod (MetaStep *mStep)
 Returns reference to receiver's numerical method. More...
 
virtual int checkConsistency ()
 Allows programmer to test some receiver's internal data, before computation begins. More...
 
virtual int requiresUnknownsDictionaryUpdate ()
 Allows to change number of equations during solution. More...
 
virtual bool requiresEquationRenumbering (TimeStep *)
 Returns true if equation renumbering is required for given solution step. More...
 
FunctiongiveDtFunction ()
 Returns time function for time step increment. More...
 
double giveDeltaT (int n)
 Returns the time step length for given step number n, initial step is number 0. More...
 
double giveDiscreteTime (int n)
 Returns time for time step number n (array discreteTimes must be specified) More...
 
- Public Member Functions inherited from oofem::StationaryTransportProblem
 StationaryTransportProblem (int i, EngngModel *_master)
 Constructor. More...
 
virtual ~StationaryTransportProblem ()
 Destructor. More...
 
virtual void updateComponent (TimeStep *tStep, NumericalCmpn cmpn, Domain *d)
 Updates components mapped to numerical method if necessary during solution process. More...
 
virtual FieldPtr giveField (FieldType key, TimeStep *)
 Returns the smart pointer to requested field, Null otherwise. More...
 
- Public Member Functions inherited from oofem::EngngModel
 EngngModel (int i, EngngModel *_master=NULL)
 Constructor. More...
 
virtual ~EngngModel ()
 Destructor. More...
 
 EngngModel (const EngngModel &)=delete
 
EngngModeloperator= (const EngngModel &)=delete
 
DomaingiveDomain (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 ErrorEstimatorgiveDomainErrorEstimator (int n)
 Service for accessing ErrorEstimator corresponding to particular domain. More...
 
virtual MaterialInterfacegiveMaterialInterface (int n)
 Returns material interface representation for given 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...
 
EngngModelgiveMasterEngngModel ()
 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...
 
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...
 
ProblemCommunicatorgiveProblemCommunicator (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...
 
MetaStepgiveCurrentMetaStep ()
 Returns current meta step. More...
 
virtual TimeStepgiveCurrentStep (bool force=false)
 Returns current time step. More...
 
virtual TimeStepgivePreviousStep (bool force=false)
 Returns previous time step. More...
 
TimeStepgenerateNextStep ()
 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...
 
MetaStepgiveMetaStep (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...
 
ExportModuleManagergiveExportModuleManager ()
 Returns receiver's export module manager. More...
 
EngngModelTimergiveTimer ()
 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 int forceEquationNumbering (int i)
 Forces equation renumbering on given domain. More...
 
virtual int forceEquationNumbering ()
 Forces equation renumbering on all domains associated to engng model. More...
 
virtual ParallelContextgiveParallelContext (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 void printDofOutputAt (FILE *stream, Dof *iDof, TimeStep *tStep)
 DOF printing routine. 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...
 
EngngModelContextgiveContext ()
 Context requesting service. More...
 
virtual int giveNumberOfSlaveProblems ()
 Returns number of slave problems. More...
 
virtual EngngModelgiveSlaveProblem (int i)
 Returns i-th slave problem. 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 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 LoadBalancergiveLoadBalancer ()
 Returns reference to receiver's load balancer. More...
 
virtual LoadBalancerMonitorgiveLoadBalancerMonitor ()
 Returns reference to receiver's load balancer monitor. More...
 
void initParallel ()
 Request domain rank and problem size. More...
 
EngngModelgiveEngngModel ()
 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 Types

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

Protected Member Functions

virtual void updateInternalState (TimeStep *tStep)
 Updates IP values on elements. More...
 
virtual void applyIC (TimeStep *tStep)
 This function is normally called at the first time to project initial conditions to previous (0^th) solution vector. More...
 
void createPreviousSolutionInDofUnknownsDictionary (TimeStep *tStep)
 
void assembleAlgorithmicPartOfRhs (FloatArray &rhs, const UnknownNumberingScheme &s, TimeStep *tStep)
 
- Protected Member Functions inherited from oofem::NonStationaryTransportProblem
virtual void assembleDirichletBcRhsVector (FloatArray &answer, TimeStep *tStep, ValueModeType mode, const UnknownNumberingScheme &s, Domain *d)
 Assembles part of RHS due to Dirichlet boundary conditions. More...
 
virtual void copyUnknownsInDictionary (ValueModeType mode, TimeStep *fromTime, TimeStep *toTime)
 Copy unknowns in DOF's from previous to current position. More...
 
- 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

double rtol
 
int nsmax
 
nlttp_ModeType NR_Mode
 
int MANRMSteps
 
int currentIterations
 
- Protected Attributes inherited from oofem::NonStationaryTransportProblem
StateCounterType internalVarUpdateStamp
 Contains last time stamp of internal variable update. More...
 
LinSystSolverType solverType
 
std::unique_ptr< SparseLinearSystemNMlinSolver
 
FloatArray bcRhs
 Right hand side vector from boundary conditions. More...
 
double initT
 Initial time from which the computation runs. Default is zero. More...
 
double deltaT
 Length of time step. More...
 
double alpha
 
int lumpedCapacityStab
 If set then stabilization using lumped capacity will be used. More...
 
int dtFunction
 Associated time function for time step increment. More...
 
FloatArray discreteTimes
 Specified times where the problem is solved. More...
 
bool changingProblemSize
 Determines if there are change in the problem size (no application/removal of Dirichlet boundary conditions). More...
 
- Protected Attributes inherited from oofem::StationaryTransportProblem
SparseMtrxType sparseMtrxType
 
std::unique_ptr< PrimaryFieldUnknownsField
 This field stores solution vector. For fixed size of problem, the PrimaryField is used, for growing/decreasing size, DofDistributedPrimaryField applies. More...
 
std::unique_ptr< SparseMtrxconductivityMatrix
 
FloatArray internalForces
 
FloatArray eNorm
 
SparseNonLinearSystemNMnMethod
 Numerical method used to solve the problem. More...
 
bool keepTangent
 
- 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< MetaStepmetaStepList
 List of problem metasteps. More...
 
std::unique_ptr< TimeStepstepWhenIcApply
 Solution step when IC (initial conditions) apply. More...
 
std::unique_ptr< TimeStepcurrentStep
 Current time step. More...
 
std::unique_ptr< TimeSteppreviousStep
 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
 
ExportModuleManagerexportModuleManager
 Export module manager. More...
 
InitModuleManagerinitModuleManager
 Initialization module manager. More...
 
problemMode pMode
 Domain mode. More...
 
problemScale pScale
 Multiscale mode. More...
 
time_t startTime
 Solution start time. More...
 
EngngModelmaster
 Master e-model; if defined receiver is in maintained (slave) mode. More...
 
EngngModelContextcontext
 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...
 
ErrorEstimatordefaultErrEstimator
 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...
 
CommunicatorBuffcommBuff
 Common Communicator buffer. More...
 
ProblemCommunicatorcommunicator
 Communicator. More...
 
ProblemCommunicatornonlocCommunicator
 NonLocal Communicator. Necessary when nonlocal constitutive models are used. More...
 
std::vector< ParallelContextparallelContextList
 List where parallel contexts are stored. More...
 
bool suppressOutput
 Flag for suppressing output to file. More...
 
std::string simulationDescription
 
LoadBalancerlb
 Load Balancer. More...
 
LoadBalancerMonitorlbm
 
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...
 

Detailed Description

This class represents nonlinear transient transport problem.

The problem can be growing/decreasing, signalized by flag "changingproblemsize" in the problem description. The solution is stored in UnknownsField, which can obtain/ project solution from/to DOFs (nodes). If the problem keeps the same equation numbers, solution is taken from UnknownsField without any projection, which is more efficient. See the matlibmanual for solution strategy of balance equations and the solution algorithm.

Todo:
Documentation errors (there is no "UnknownsField" used here).

Definition at line 59 of file nltransienttransportproblem.h.

Member Enumeration Documentation

Enumerator
nrsolverModifiedNRM 
nrsolverFullNRM 
nrsolverAccelNRM 

Definition at line 62 of file nltransienttransportproblem.h.

Constructor & Destructor Documentation

oofem::NLTransientTransportProblem::NLTransientTransportProblem ( int  i,
EngngModel _master = NULL 
)

Constructor.

Definition at line 50 of file nltransienttransportproblem.C.

oofem::NLTransientTransportProblem::~NLTransientTransportProblem ( )
virtual

Destructor.

Definition at line 55 of file nltransienttransportproblem.C.

Member Function Documentation

void oofem::NLTransientTransportProblem::applyIC ( TimeStep tStep)
protectedvirtual

This function is normally called at the first time to project initial conditions to previous (0^th) solution vector.

Parameters
tStepPrevious solution step.

Reimplemented from oofem::NonStationaryTransportProblem.

Definition at line 282 of file nltransienttransportproblem.C.

References oofem::NonStationaryTransportProblem::applyIC(), oofem::EngngModel::giveDomain(), oofem::Domain::giveElements(), oofem::TransportElement::updateInternalState(), and oofem::Element::updateYourself().

Referenced by giveCurrentNumberOfIterations(), and solveYourselfAt().

void oofem::NLTransientTransportProblem::createPreviousSolutionInDofUnknownsDictionary ( TimeStep tStep)
protected
virtual const char* oofem::NLTransientTransportProblem::giveClassName ( ) const
inlinevirtual

Returns class name of the receiver.

Reimplemented from oofem::NonStationaryTransportProblem.

Definition at line 85 of file nltransienttransportproblem.h.

virtual int oofem::NLTransientTransportProblem::giveCurrentNumberOfIterations ( )
inlinevirtual

Returns number of iterations that was required to reach equilibrium - used for adaptive step length in staggered problem.

Reimplemented from oofem::EngngModel.

Definition at line 90 of file nltransienttransportproblem.h.

References applyIC(), assembleAlgorithmicPartOfRhs(), createPreviousSolutionInDofUnknownsDictionary(), currentIterations, and updateInternalState().

virtual fMode oofem::NLTransientTransportProblem::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::NonStationaryTransportProblem.

Definition at line 86 of file nltransienttransportproblem.h.

References giveUnknownDictHashIndx(), oofem::EngngModel::nonLinFormulation, and updateDofUnknownsDictionary().

virtual const char* oofem::NLTransientTransportProblem::giveInputRecordName ( ) const
inlinevirtual
TimeStep * oofem::NLTransientTransportProblem::giveNextStep ( )
virtual
double oofem::NLTransientTransportProblem::giveUnknownComponent ( ValueModeType  ,
TimeStep ,
Domain ,
Dof  
)
virtual
int oofem::NLTransientTransportProblem::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::NonStationaryTransportProblem.

Definition at line 327 of file nltransienttransportproblem.C.

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

Referenced by giveFormulation(), and giveUnknownComponent().

IRResultType oofem::NLTransientTransportProblem::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::NonStationaryTransportProblem.

Definition at line 61 of file nltransienttransportproblem.C.

References _IFT_NLTransientTransportProblem_manrmsteps, _IFT_NLTransientTransportProblem_nsmax, _IFT_NLTransientTransportProblem_rtol, oofem::NonStationaryTransportProblem::initializeFrom(), IR_GIVE_FIELD, IR_GIVE_OPTIONAL_FIELD, oofem::IRRT_OK, MANRMSteps, NR_Mode, nrsolverAccelNRM, nrsolverModifiedNRM, nsmax, and rtol.

void oofem::NLTransientTransportProblem::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::NonStationaryTransportProblem.

Definition at line 97 of file nltransienttransportproblem.C.

References oofem::NonStationaryTransportProblem::alpha, applyIC(), oofem::EngngModel::assemble(), assembleAlgorithmicPartOfRhs(), oofem::EngngModel::assembleVectorFromDofManagers(), oofem::EngngModel::assembleVectorFromElements(), oofem::NonStationaryTransportProblem::changingProblemSize, oofem::classFactory, oofem::FloatArray::computeNorm(), oofem::StationaryTransportProblem::conductivityMatrix, oofem::NonStationaryTransportProblem::copyUnknownsInDictionary(), oofem::ClassFactory::createSparseMtrx(), currentIterations, oofem::EngngModel::giveCurrentMetaStep(), oofem::EngngModel::giveDomain(), oofem::TimeStep::giveMetaStepNumber(), oofem::TimeStep::giveNumber(), oofem::EngngModel::giveNumberOfDomainEquations(), oofem::NonStationaryTransportProblem::giveNumericalMethod(), oofem::TimeStep::givePreviousStep(), oofem::TimeStep::giveSolutionStateCounter(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::incrementStateCounter(), oofem::TimeStep::isTheFirstStep(), oofem::NonStationaryTransportProblem::linSolver, oofem::NonStationaryTransportProblem::lumpedCapacityStab, MANRMSteps, NR_Mode, nrsolverAccelNRM, nrsolverFullNRM, nsmax, OOFEM_ERROR, OOFEM_LOG_INFO, OOFEM_LOG_RELEVANT, oofem::FloatArray::resize(), rtol, oofem::StationaryTransportProblem::sparseMtrxType, oofem::EngngModel::stepWhenIcApply, oofem::StationaryTransportProblem::UnknownsField, updateInternalState(), and oofem::FloatArray::zero().

void oofem::NLTransientTransportProblem::updateDofUnknownsDictionary ( DofManager ,
TimeStep  
)
virtual

Updates necessary values in Dofs unknown dictionaries.

See also
EngngModel::requiresUnknownsDictionaryUpdate
Dof::updateUnknownsDictionary

Reimplemented from oofem::EngngModel.

Definition at line 345 of file nltransienttransportproblem.C.

References oofem::FloatArray::at(), and oofem::StationaryTransportProblem::UnknownsField.

Referenced by giveFormulation(), and updateInternalState().

void oofem::NLTransientTransportProblem::updateInternalState ( TimeStep tStep)
protectedvirtual
void oofem::NLTransientTransportProblem::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::NonStationaryTransportProblem.

Definition at line 314 of file nltransienttransportproblem.C.

References oofem::NonStationaryTransportProblem::alpha, oofem::TimeStep::giveIntrinsicTime(), oofem::TimeStep::giveTargetTime(), oofem::TimeStep::giveTimeIncrement(), oofem::TimeStep::setIntrinsicTime(), and oofem::NonStationaryTransportProblem::updateYourself().

Member Data Documentation

int oofem::NLTransientTransportProblem::currentIterations
protected

Definition at line 68 of file nltransienttransportproblem.h.

Referenced by giveCurrentNumberOfIterations(), and solveYourselfAt().

int oofem::NLTransientTransportProblem::MANRMSteps
protected

Definition at line 67 of file nltransienttransportproblem.h.

Referenced by initializeFrom(), and solveYourselfAt().

nlttp_ModeType oofem::NLTransientTransportProblem::NR_Mode
protected

Definition at line 66 of file nltransienttransportproblem.h.

Referenced by initializeFrom(), and solveYourselfAt().

int oofem::NLTransientTransportProblem::nsmax
protected

Definition at line 65 of file nltransienttransportproblem.h.

Referenced by initializeFrom(), and solveYourselfAt().

double oofem::NLTransientTransportProblem::rtol
protected

Definition at line 64 of file nltransienttransportproblem.h.

Referenced by initializeFrom(), and solveYourselfAt().


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:38 for OOFEM by doxygen 1.8.11 written by Dimitri van Heesch, © 1997-2011