geos::FlowSolverBase Class Reference

#include <FlowSolverBase.hpp>

Inheritance diagram for geos::FlowSolverBase:

Classes
class	BCMessage
	Class used for displaying boundary warning message. More...
struct	viewKeyStruct

Public Member Functions
	FlowSolverBase (const string &name, Group *const parent)
	main constructor for Group Objects More...
	FlowSolverBase ()=delete
	deleted default constructor
	FlowSolverBase (FlowSolverBase const &)=delete
	deleted copy constructor
	FlowSolverBase (FlowSolverBase &&)=default
	default move constructor
FlowSolverBase &	operator= (FlowSolverBase const &)=delete
	deleted assignment operator
FlowSolverBase &	operator= (FlowSolverBase &&)=delete
	deleted move operator
virtual void	registerDataOnMesh (Group &MeshBodies) override
	Register wrappers that contain data on the mesh objects. More...
void	prepareStencilWeights (DomainPartition &domain) const
	Prepare the stencil weights by removing the contribution of the hydraulic aperture before the aperture is updated. More...
void	updateStencilWeights (DomainPartition &domain) const
	Update the stencil weights by adding the contribution of the hydraulic aperture after the aperture is updated. More...
void	enableFixedStressPoromechanicsUpdate ()
void	enableJumpStabilization ()
void	updatePorosityAndPermeability (CellElementSubRegion &subRegion) const
virtual void	updatePorosityAndPermeability (SurfaceElementSubRegion &subRegion) const
virtual void	saveSequentialIterationState (DomainPartition &domain) override
	Utility function to save the iteration state (useful for sequential simulations) More...
void	findMinMaxElevationInEquilibriumTarget (DomainPartition &domain, std::map< string, localIndex > const &equilNameToEquilId, arrayView1d< real64 > const &maxElevation, arrayView1d< real64 > const &minElevation) const
	For each equilibrium initial condition, loop over all the target cells and compute the min/max elevation. More...
void	computeSourceFluxSizeScalingFactor (real64 const &time, real64 const &dt, DomainPartition &domain, std::map< string, localIndex > const &bcNameToBcId, arrayView1d< globalIndex > const &bcAllSetsSize) const
	For each source flux boundary condition, loop over all the target cells and sum the owned cells. More...
integer &	isThermal ()
virtual units::Unit	getMassUnit () const
void	allowNegativePressure ()
	Function to activate the flag allowing negative pressure.
void	setKeepVariablesConstantDuringInitStep (bool const keepVariablesConstantDuringInitStep)
	Utility function to keep the flow variables during a time step (used in poromechanics simulations) More...
virtual bool	checkSequentialSolutionIncrements (DomainPartition &domain) const override
	Check if the solution increments are ok to use. More...
void	enableLaggingFractureStencilWeightsUpdate ()
real64	sumAquiferFluxes (BoundaryStencil const &stencil, AquiferBoundaryCondition::KernelWrapper const &aquiferBCWrapper, ElementViewConst< arrayView1d< real64 const > > const &pres, ElementViewConst< arrayView1d< real64 const > > const &presOld, ElementViewConst< arrayView1d< real64 const > > const &gravCoef, real64 const &timeAtBeginningOfStep, real64 const &dt)
Public Member Functions inherited from geos::PhysicsSolverBase
	PhysicsSolverBase (string const &name, Group *const parent)
	Constructor for PhysicsSolverBase. More...
	PhysicsSolverBase (PhysicsSolverBase &&)=default
	Move constructor for PhysicsSolverBase.
virtual	~PhysicsSolverBase () override
	Destructor for PhysicsSolverBase.
	PhysicsSolverBase ()=delete
	Deleted constructor.
	PhysicsSolverBase (PhysicsSolverBase const &)=delete
	Deleted copy constructor.
PhysicsSolverBase &	operator= (PhysicsSolverBase const &)=delete
	Deleted copy assignment operator.
PhysicsSolverBase &	operator= (PhysicsSolverBase &&)=delete
	Deleted move assignment operator.
virtual string	getCatalogName () const =0
virtual void	initialize_postMeshGeneration () override
	Initialization tasks after mesh generation is completed.
void	generateMeshTargetsFromTargetRegions (Group const &meshBodies)
	Generate mesh targets from target regions. More...
virtual void	cleanup (real64 const time_n, integer const cycleNumber, integer const eventCounter, real64 const eventProgress, DomainPartition &domain) override
	Called as the code exits the main run loop. More...
virtual bool	execute (real64 const time_n, real64 const dt, integer const cycleNumber, integer const eventCounter, real64 const eventProgress, DomainPartition &domain) override
	Main extension point of executable targets. More...
ParallelMatrix &	getSystemMatrix ()
	Getter for system matrix. More...
ParallelMatrix const &	getSystemMatrix () const
	Getter for system rhs vector. More...
ParallelVector &	getSystemRhs ()
	Getter for system rhs vector. More...
ParallelVector const &	getSystemRhs () const
	Getter for system rhs vector. More...
ParallelVector &	getSystemSolution ()
	Getter for system solution vector. More...
ParallelVector const &	getSystemSolution () const
	Getter for system solution vector. More...
DofManager &	getDofManager ()
	Getter for degree-of-freedom manager. More...
DofManager const &	getDofManager () const
	Getter for degree-of-freedom manager. More...
CRSMatrix< real64, globalIndex > &	getLocalMatrix ()
	Getter for local matrix. More...
CRSMatrixView< real64 const, globalIndex const >	getLocalMatrix () const
	Getter for local matrix. More...
virtual real64	solverStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
	entry function to perform a solver step More...
virtual real64	setNextDt (real64 const &currentDt, DomainPartition &domain)
	function to set the next time step size More...
virtual real64	setNextDtBasedOnNewtonIter (real64 const &currentDt)
	function to set the next time step size based on Newton convergence More...
virtual real64	setNextDtBasedOnStateChange (real64 const &currentDt, DomainPartition &domain)
	function to set the next dt based on state change More...
virtual real64	setNextDtBasedOnCFL (real64 const &currentDt, DomainPartition &domain)
	function to set the next dt based on state change More...
virtual real64	explicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
	Entry function for an explicit time integration step. More...
virtual real64	nonlinearImplicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
	Function for a nonlinear implicit integration step. More...
virtual bool	lineSearch (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, ParallelVector &rhs, ParallelVector &solution, real64 const scaleFactor, real64 &lastResidual)
	Function to perform line search. More...
virtual bool	lineSearchWithParabolicInterpolation (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, ParallelVector &rhs, ParallelVector &solution, real64 const scaleFactor, real64 &lastResidual, real64 &residualNormT)
	Function to perform line search using a parabolic interpolation to find the scaling factor. More...
virtual real64	linearImplicitStep (real64 const &time_n, real64 const &dt, integer const cycleNumber, DomainPartition &domain)
	Function for a linear implicit integration step. More...
virtual void	implicitStepSetup (real64 const &time_n, real64 const &dt, DomainPartition &domain)
	function to perform setup for implicit timestep More...
virtual void	setupDofs (DomainPartition const &domain, DofManager &dofManager) const
	Populate degree-of-freedom manager with fields relevant to this solver. More...
virtual void	setupSystem (DomainPartition &domain, DofManager &dofManager, CRSMatrix< real64, globalIndex > &localMatrix, ParallelVector &rhs, ParallelVector &solution, bool const setSparsity=true)
	Set up the linear system (DOF indices and sparsity patterns) More...
virtual void	assembleSystem (real64 const time, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs)
	function to assemble the linear system matrix and rhs More...
virtual void	applyBoundaryConditions (real64 const time, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs)
	apply boundary condition to system More...
void	debugOutputSystem (real64 const &time, integer const cycleNumber, integer const nonlinearIteration, ParallelMatrix const &matrix, ParallelVector const &rhs) const
	Output the assembled linear system for debug purposes. More...
void	debugOutputSolution (real64 const &time, integer const cycleNumber, integer const nonlinearIteration, ParallelVector const &solution) const
	Output the linear system solution for debug purposes. More...
virtual real64	calculateResidualNorm (real64 const &time, real64 const &dt, DomainPartition const &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localRhs)
	calculate the norm of the global system residual More...
virtual void	solveLinearSystem (DofManager const &dofManager, ParallelMatrix &matrix, ParallelVector &rhs, ParallelVector &solution)
	function to apply a linear system solver to the assembled system. More...
virtual bool	checkSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor)
	Function to check system solution for physical consistency and constraint violation. More...
virtual real64	scalingForSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution)
	Function to determine if the solution vector should be scaled back in order to maintain a known constraint. More...
virtual void	applySystemSolution (DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor, real64 const dt, DomainPartition &domain)
	Function to apply the solution vector to the state. More...
virtual bool	updateConfiguration (DomainPartition &domain)
	updates the configuration (if needed) based on the state after a converged Newton loop. More...
virtual void	outputConfigurationStatistics (DomainPartition const &domain) const
virtual void	resetConfigurationToBeginningOfStep (DomainPartition &domain)
	resets the configuration to the beginning of the time-step. More...
virtual bool	resetConfigurationToDefault (DomainPartition &domain) const
	resets the configuration to the default value. More...
virtual void	updateState (DomainPartition &domain)
	Recompute all dependent quantities from primary variables (including constitutive models) More...
virtual void	resetStateToBeginningOfStep (DomainPartition &domain)
	reset state of physics back to the beginning of the step. More...
virtual void	implicitStepComplete (real64 const &time, real64 const &dt, DomainPartition &domain)
	perform cleanup for implicit timestep More...
virtual real64	getTimestepRequest (real64 const) override
	getter for the next timestep size More...
real64	getTimestepRequest ()
	getter for the next timestep size More...
virtual Group *	createChild (string const &childKey, string const &childName) override
	creates a child group of of this PhysicsSolverBase instantiation More...
Timestamp	getSystemSetupTimestamp () const
	getter for the timestamp of the system setup More...
Timestamp	getMeshModificationTimestamp (DomainPartition &domain) const
	getter for the timestamp of the mesh modification on the mesh levels More...
void	setSystemSetupTimestamp (Timestamp timestamp)
	set the timestamp of the system setup More...
R1Tensor const	gravityVector () const
	return the value of the gravity vector specified in PhysicsSolverManager More...
LinearSolverParameters &	getLinearSolverParameters ()
	accessor for the linear solver parameters. More...
LinearSolverParameters const &	getLinearSolverParameters () const
	const accessor for the linear solver parameters. More...
NonlinearSolverParameters &	getNonlinearSolverParameters ()
	accessor for the nonlinear solver parameters. More...
NonlinearSolverParameters const &	getNonlinearSolverParameters () const
	const accessor for the nonlinear solver parameters. More...
virtual void	synchronizeNonlinearSolverParameters ()
	syncronize the nonlinear solver parameters.
localIndex	targetRegionIndex (string const &regionName) const
	Get position of a given region within solver's target region list. More...
template<typename LAMBDA >
void	forDiscretizationOnMeshTargets (Group const &meshBodies, LAMBDA &&lambda) const
	Loop over the target discretization on all mesh targets and apply callback. More...
template<typename LAMBDA >
void	forDiscretizationOnMeshTargets (Group &meshBodies, LAMBDA &&lambda) const
	Loop over the target discretization on all mesh targets and apply callback. More...
string	getDiscretizationName () const
	return the name of the discretization object More...
virtual bool	registerCallback (void *func, const std::type_info &funcType) final override
	function to set the value of m_assemblyCallback More...
SolverStatistics &	getSolverStatistics ()
	accessor for the solver statistics. More...
SolverStatistics const &	getSolverStatistics () const
	const accessor for the solver statistics. More...
map< std::pair< string, string >, array1d< string > > const &	getMeshTargets () const
	accessor for m_meshTargets More...
Public Member Functions inherited from geos::ExecutableGroup
virtual void	signalToPrepareForExecution (real64 const time_n, real64 const dt, integer const cycle, DomainPartition &domain)
	Inform the object that it expects to execute during the next timestep. More...
void	setTimesteppingBehavior (TimesteppingBehavior const timesteppingBehavior)
	Set the timestep behavior for a target. More...
TimesteppingBehavior	getTimesteppingBehavior () const
	Get the target's time step behavior. More...
Public Member Functions inherited from geos::dataRepository::Group
template<typename LOG_LEVEL_INFO >
std::enable_if_t< geos::is_log_level_info< LOG_LEVEL_INFO >, void >	addLogLevel ()
	Append a levelCondition and a log description to the description of the wrapped object given a log info struct. Must be called in constructor. More...
virtual void	reinit ()
	Performs re-initialization of certain variable depending on the solver being used.
	Group (string const &name, Group *const parent)
	Constructor. More...
	Group (string const &name, conduit::Node &rootNode)
	Constructor. More...
	Group (Group &&source)=default
	Move constructor. More...
virtual	~Group ()
	Destructor, deletes all Groups and Wrappers owned by this Group.
	Group ()=delete
	Deleted default constructor.
	Group (Group const &)=delete
	Deleted copy constructor.
Group &	operator= (Group const &)=delete
	Deleted copy assignment operator. More...
Group &	operator= (Group &&)=delete
	Deleted move assignment operator. More...
void	printDataHierarchy (integer indent=0) const
	Prints the data hierarchy recursively. More...
string	dumpInputOptions () const
string	dumpSubGroupsNames () const
string	dumpWrappersNames () const
template<typename T = Group>
T &	registerGroup (string const &name, std::unique_ptr< T > newObject)
	Register a new Group as a sub-group of current Group. More...
template<typename T = Group>
T &	registerGroup (string const &name, T *newObject)
	Register a new Group as a sub-group of current Group. More...
template<typename T = Group>
T &	registerGroup (string const &name)
	Register a new Group as a sub-group of current Group. More...
template<typename T = Group>
T &	registerGroup (subGroupMap::KeyIndex const &keyIndex)
	Register a new Group as a sub-group of current Group. More...
template<typename T = Group, typename TBASE = Group>
T &	registerGroup (string const &name, string const &catalogName)
	Register a new Group as a sub-group of current Group. More...
void	deregisterGroup (string const &name)
	Removes a child group from this group. More...
template<typename T = Group, typename KEY = void>
T *	getGroupPointer (KEY const &key)
	Return a pointer to a sub-group of the current Group. More...
template<typename T = Group, typename KEY = void>
T const *	getGroupPointer (KEY const &key) const
	Return a pointer to a sub-group of the current Group. More...
template<typename T = Group, typename KEY = void>
T &	getGroup (KEY const &key)
	Return a reference to a sub-group of the current Group. More...
template<typename T = Group, typename KEY = void>
T const &	getGroup (KEY const &key) const
	Return a reference to a sub-group of the current Group. More...
template<typename T = Group>
T &	getGroupByPath (string const &path)
	Retrieve a group from the hierarchy using a path. More...
template<typename T = Group>
T const &	getGroupByPath (string const &path) const
	Retrieve a group from the hierarchy using a path. More...
subGroupMap &	getSubGroups ()
	Get the subgroups object. More...
subGroupMap const &	getSubGroups () const
	Get the subgroups object. More...
localIndex	numSubGroups () const
	return the number of sub groups in this Group More...
std::vector< string >	getSubGroupsNames () const
template<typename T = Group>
bool	hasGroup (string const &name) const
	Check whether a sub-group exists. More...
template<typename T >
bool	hasSubGroupOfType () const
	Check whether a sub-group exists by type. More...
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void	forSubGroups (LAMBDA &&lambda)
	Apply the given functor to subgroups that can be casted to one of specified types. More...
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void	forSubGroups (LAMBDA &&lambda) const
	Apply the given functor to subgroups that can be casted to one of specified types. More...
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void	forSubGroupsIndex (LAMBDA &&lambda)
	Apply the given functor to subgroups that can be casted to one of specified types. More...
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LAMBDA >
void	forSubGroupsIndex (LAMBDA &&lambda) const
	Apply the given functor to subgroups that can be casted to one of specified types. More...
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LOOKUP_CONTAINER , typename LAMBDA >
void	forSubGroups (LOOKUP_CONTAINER const &subGroupKeys, LAMBDA &&lambda)
	Apply the given functor to subgroups that can be casted to one of specified types. More...
template<typename GROUPTYPE = Group, typename ... GROUPTYPES, typename LOOKUP_CONTAINER , typename LAMBDA >
void	forSubGroups (LOOKUP_CONTAINER const &subGroupKeys, LAMBDA &&lambda) const
	Apply the given functor to subgroups that can be casted to one of specified types. More...
template<typename LAMBDA >
void	forWrappers (LAMBDA &&lambda)
	Apply the given functor to wrappers. More...
template<typename LAMBDA >
void	forWrappers (LAMBDA &&lambda) const
	Apply the given functor to wrappers. More...
template<typename TYPE , typename ... TYPES, typename LAMBDA >
void	forWrappers (LAMBDA &&lambda)
	Apply the given functor to wrappers that can be cast to one of specified types. More...
template<typename TYPE , typename ... TYPES, typename LAMBDA >
void	forWrappers (LAMBDA &&lambda) const
	Apply the given functor to wrappers that can be cast to one of specified types. More...
void	initialize ()
	Run initialization functions on this and all subgroups. More...
virtual void	initializationOrder (string_array &order)
	Sets the initialization order for sub-Groups. More...
void	initializePostInitialConditions ()
	Initialization routine to be called after calling ApplyInitialConditions(). More...
void	postRestartInitializationRecursive ()
	Initialization routine to be called after calling reading a restart file. More...
void	processInputFileRecursive (xmlWrapper::xmlDocument &xmlDocument, xmlWrapper::xmlNode &targetNode)
	Recursively read values using ProcessInputFile() from the input file and put them into the wrapped values for this group. Also add the includes content to the xmlDocument when Include nodes are encountered. More...
void	processInputFileRecursive (xmlWrapper::xmlDocument &xmlDocument, xmlWrapper::xmlNode &targetNode, xmlWrapper::xmlNodePos const &nodePos)
	Same as processInputFileRecursive(xmlWrapper::xmlDocument &, xmlWrapper::xmlNode &) but allow to reuse an existing xmlNodePos. More...
void	postInputInitializationRecursive ()
	Recursively call postInputInitialization() to apply post processing after reading input values.
template<typename T , typename TBASE = T>
Wrapper< TBASE > &	registerWrapper (string const &name, wrapperMap::KeyIndex::index_type *const rkey=nullptr)
	Create and register a Wrapper around a new object. More...
template<typename T , typename TBASE = T>
Wrapper< TBASE > &	registerWrapper (Group::wrapperMap::KeyIndex const &viewKey)
	Create and register a Wrapper around a new object. More...
template<typename T >
Wrapper< T > &	registerWrapper (string const &name, std::unique_ptr< T > newObject)
	Register a Wrapper around a given object and take ownership. More...
template<typename T >
Wrapper< T > &	registerWrapper (string const &name, T *newObject)
	Register a Wrapper around an existing object, does not take ownership of the object. More...
WrapperBase &	registerWrapper (std::unique_ptr< WrapperBase > wrapper)
	Register and take ownership of an existing Wrapper. More...
void	deregisterWrapper (string const &name)
	Removes a Wrapper from this group. More...
void	generateDataStructureSkeleton (integer const level)
	Build a complete datastructure for schema generation. More...
virtual void	expandObjectCatalogs ()
	Expand any catalogs in the data structure.
virtual void	setSchemaDeviations (xmlWrapper::xmlNode schemaRoot, xmlWrapper::xmlNode schemaParent, integer documentationType)
	Inform the schema generator of any deviations between the xml and GEOS data structures. More...
virtual void	registerDataOnMeshRecursive (Group &meshBodies)
	Calls RegisterDataOnMesh() recursively. More...
virtual localIndex	packSize (string_array const &wrapperNames, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
	Get the size required to pack a list of wrappers. More...
virtual localIndex	packSize (string_array const &wrapperNames, arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
	Get the size required to pack a list of indices within a list of wrappers. More...
localIndex	packSize (arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
	Get the size required to pack a list of indices for all registered wrappers. More...
virtual localIndex	pack (buffer_unit_type *&buffer, string_array const &wrapperNames, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
	Pack a list of wrappers to a buffer. More...
virtual localIndex	pack (buffer_unit_type *&buffer, string_array const &wrapperNames, arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
	Pack a list of indices within a list of wrappers. More...
localIndex	pack (buffer_unit_type *&buffer, arrayView1d< localIndex const > const &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events) const
	Pack a list of indices for all registered wrappers. More...
virtual localIndex	unpack (buffer_unit_type const *&buffer, arrayView1d< localIndex > &packList, integer const recursive, bool onDevice, parallelDeviceEvents &events, MPI_Op op=MPI_REPLACE)
	Unpack a buffer. More...
template<typename KEY >
WrapperBase const &	getWrapperBase (KEY const &key) const
	Return a reference to a WrapperBase stored in this group. More...
template<typename KEY >
WrapperBase &	getWrapperBase (KEY const &key)
	Return a reference to a WrapperBase stored in this group. More...
indexType	getWrapperIndex (string const &name) const
wrapperMap const &	wrappers () const
	Get access to the internal wrapper storage. More...
wrapperMap &	wrappers ()
	Get access to the internal wrapper storage. More...
indexType	numWrappers () const
	Return the number of wrappers. More...
std::vector< string >	getWrappersNames () const
template<typename LOOKUP_TYPE >
bool	hasWrapper (LOOKUP_TYPE const &lookup) const
	Check if a wrapper exists. More...
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > const &	getWrapper (LOOKUP_TYPE const &index) const
	Retrieve a Wrapper stored in this group. More...
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > &	getWrapper (LOOKUP_TYPE const &index)
	Retrieve a Wrapper stored in this group. More...
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > const *	getWrapperPointer (LOOKUP_TYPE const &index) const
	Retrieve a Wrapper stored in this group. More...
template<typename T , typename LOOKUP_TYPE >
Wrapper< T > *	getWrapperPointer (LOOKUP_TYPE const &index)
	Retrieve a Wrapper stored in this group. More...
template<typename T , typename LOOKUP_TYPE >
GEOS_DECLTYPE_AUTO_RETURN	getReference (LOOKUP_TYPE const &lookup) const
	Look up a wrapper and get reference to wrapped object. More...
template<typename T , typename LOOKUP_TYPE >
T &	getReference (LOOKUP_TYPE const &lookup)
	Look up a wrapper and get reference to wrapped object. More...
virtual void	resize (localIndex const newSize)
	Resize the group and all contained wrappers that resize with parent. More...
virtual void	reserve (indexType const newsize)
	Set the new capacity and reserve it in all wrappers that resize with parent. More...
localIndex	capacity () const
	Get the "capacity" of the group, which determines the capacity of resizable wrappers. More...
localIndex	size () const
	Get the "size" of the group, which determines the number of elements in resizable wrappers. More...
string const &	getName () const
	Get group name. More...
string	getPath () const
	Return the path of this Group in the data repository. Starts with '/' followed by the hierarchy of the children of the "Problem" in which the Group is. More...
DataContext const &	getDataContext () const
template<typename KEY >
DataContext const &	getWrapperDataContext (KEY key) const
Group &	getParent ()
	Access the group's parent. More...
Group const &	getParent () const
	Access the group's parent. More...
bool	hasParent () const
localIndex	getIndexInParent () const
	Get the group's index within its parent group. More...
localIndex	getSubGroupIndex (keyType const &key) const
	Get the index of a sub-Group within this group. More...
int	sizedFromParent () const
	Check whether this Group is resized when its parent is resized. More...
Group &	setSizedFromParent (int val)
	Set whether this wrapper is resized when its parent is resized. More...
RestartFlags	getRestartFlags () const
	Get flags that control restart output of this group. More...
void	setRestartFlags (RestartFlags flags)
	Set flags that control restart output of this group. More...
InputFlags	getInputFlags () const
	Get input flags for schema generation. More...
void	setInputFlags (InputFlags flags)
	Set input flags for schema generation. More...
conduit::Node &	getConduitNode ()
	Return the Conduit node object associated with this group. More...
conduit::Node const &	getConduitNode () const
	Return the Conduit node object associated with this group. More...
void	prepareToWrite ()
	Register the group and its wrappers with Conduit.
void	finishWriting ()
	Write the group and its wrappers into Conduit.
void	loadFromConduit ()
	Read the group and its wrappers from Conduit.
void	enableLogLevelInput ()
void	setLogLevel (integer const logLevel)
	Set verbosity level. More...
integer	getLogLevel () const

Static Public Member Functions
static string	coupledSolverAttributePrefix ()
	String used to form the solverName used to register single-physics solvers in CoupledSolver.
Static Public Member Functions inherited from geos::PhysicsSolverBase
static CatalogInterface::CatalogType &	getCatalog ()
	Get the singleton catalog for PhysicsSolverBase. More...
Static Public Member Functions inherited from geos::dataRepository::Group
template<typename T0 , typename T1 , typename ... CASTTYPES, typename CONTAINERTYPE , typename LAMBDA >
static bool	applyLambdaToContainer (CONTAINERTYPE &container, LAMBDA &&lambda)
	Apply a given functor to a container if the container can be cast to one of the specified types. More...
static CatalogInterface::CatalogType &	getCatalog ()
	Get the singleton catalog for this Group. More...

Protected Member Functions
virtual void	saveAquiferConvergedState (real64 const &time, real64 const &dt, DomainPartition &domain)
	Increment the cumulative flux from each aquifer. More...
virtual void	saveConvergedState (ElementSubRegionBase &subRegion) const
	Utility function to save the converged state. More...
virtual void	validatePoreVolumes (DomainPartition const &domain) const
	Helper function to compute/report the elements with small pore volumes. More...
virtual void	precomputeData (MeshLevel &mesh, arrayView1d< string const > const &regionNames)
virtual void	initializePreSubGroups () override
	Called by Initialize() prior to initializing sub-Groups.
virtual void	initializePostInitialConditionsPreSubGroups () override
	Called by InitializePostInitialConditions() prior to initializing sub-Groups.
virtual void	setConstitutiveNamesCallSuper (ElementSubRegionBase &subRegion) const override
	This function sets constitutive name fields on an ElementSubRegionBase, and calls the base function it overrides. More...
Protected Member Functions inherited from geos::PhysicsSolverBase
real64	eisenstatWalker (real64 const newNewtonNorm, real64 const oldNewtonNorm, LinearSolverParameters::Krylov const &krylovParams)
	Eisenstat-Walker adaptive tolerance. More...
Protected Member Functions inherited from geos::dataRepository::Group
virtual void	postInputInitialization ()
virtual void	initializePostSubGroups ()
	Called by Initialize() after to initializing sub-Groups.
virtual void	initializePostInitialConditionsPostSubGroups ()
	Called by InitializePostInitialConditions() after to initializing sub-Groups.
virtual void	postRestartInitialization ()
	Performs initialization required after reading from a restart file.

Protected Attributes
integer	m_numDofPerCell
	the number of Degrees of Freedom per cell
integer	m_isThermal
	flag to determine whether or not this is a thermal simulation
real64	m_inputTemperature
	the input temperature
integer	m_keepVariablesConstantDuringInitStep
	flag to freeze the initial state during initialization in coupled problems
bool	m_isFixedStressPoromechanicsUpdate
	enable the fixed stress poromechanics update of porosity
bool	m_isJumpStabilized
	enable pressure jump stabilzation for fixed-stress poromechanics
integer	m_allowNegativePressure
	flag if negative pressure is allowed
real64	m_maxAbsolutePresChange
	maximum (absolute) pressure change in a Newton iteration
real64	m_sequentialPresChange
	maximum (absolute) pressure change in a sequential iteration
real64	m_maxSequentialPresChange
real64	m_sequentialTempChange
	maximum (absolute) temperature change in a sequential iteration
real64	m_maxSequentialTempChange
Protected Attributes inherited from geos::PhysicsSolverBase
real64	m_cflFactor
	Courant–Friedrichs–Lewy factor for the timestep.
real64	m_maxStableDt
	maximum stable time step
real64	m_nextDt
	timestep of the next cycle
integer	m_numTimestepsSinceLastDtCut
	Number of cycles since last timestep cut.
string	m_discretizationName
	name of the FV discretization object in the data repository
DofManager	m_dofManager
	Data structure to handle degrees of freedom.
ParallelMatrix	m_matrix
	System matrix.
ParallelVector	m_rhs
	System right-hand side vector.
ParallelVector	m_solution
	System solution vector.
CRSMatrix< real64, globalIndex >	m_localMatrix
	Local system matrix and rhs.
std::unique_ptr< PreconditionerBase< LAInterface > >	m_precond
	Custom preconditioner for the "native" iterative solver.
integer	m_writeLinearSystem
	flag for debug output of matrix, rhs, and solution
LinearSolverParametersInput	m_linearSolverParameters
	Linear solver parameters.
LinearSolverResult	m_linearSolverResult
	Result of the last linear solve.
NonlinearSolverParameters	m_nonlinearSolverParameters
	Nonlinear solver parameters.
SolverStatistics	m_solverStatistics
	Solver statistics.
Timestamp	m_systemSetupTimestamp
	Timestamp of the last call to setup system.
std::function< void(CRSMatrix< real64, globalIndex >, array1d< real64 >) >	m_assemblyCallback
	Callback function for assembly step.
std::map< std::string, std::chrono::system_clock::duration >	m_timers
	Timers for the aggregate profiling of the solver.

Additional Inherited Members
Public Types inherited from geos::PhysicsSolverBase
using	CatalogInterface = dataRepository::CatalogInterface< PhysicsSolverBase, string const &, Group *const >
	Type alias for catalog interface used by this class. See CatalogInterface.
Public Types inherited from geos::ExecutableGroup
enum class	TimesteppingBehavior : integer { DeterminesTimeStepSize , DoesNotDetermineTimeStepSize }
	Timestepping type. More...
Public Types inherited from geos::dataRepository::Group
using	subGroupMap = MappedVector< Group, Group *, keyType, indexType >
	The template specialization of MappedVector to use for the collection of sub-Group objects.
using	wrapperMap = MappedVector< WrapperBase, WrapperBase *, keyType, indexType >
	The template specialization of MappedVector to use for the collection wrappers objects.
using	CatalogInterface = dataRepository::CatalogInterface< Group, string const &, Group *const >
	Type alias for catalog interface used by this class. See CatalogInterface.
Static Protected Member Functions inherited from geos::PhysicsSolverBase
template<typename CONSTITUTIVE_BASE_TYPE >
static string	getConstitutiveName (ElementSubRegionBase const &subRegion)
	Get the Constitutive Name object. More...
template<typename CONSTITUTIVE_BASE_TYPE >
static string	getConstitutiveName (ParticleSubRegionBase const &subRegion)
	Get the Constitutive Name object. More...
template<typename BASETYPE = constitutive::ConstitutiveBase, typename LOOKUP_TYPE >
static BASETYPE const &	getConstitutiveModel (dataRepository::Group const &dataGroup, LOOKUP_TYPE const &key)
	Get the Constitutive Model object. More...
template<typename BASETYPE = constitutive::ConstitutiveBase, typename LOOKUP_TYPE >
static BASETYPE &	getConstitutiveModel (dataRepository::Group &dataGroup, LOOKUP_TYPE const &key)
	Get the Constitutive Model object. More...

Detailed Description

Base class for finite volume fluid flow solvers. Provides some common features

Definition at line 37 of file FlowSolverBase.hpp.

Constructor & Destructor Documentation

FlowSolverBase()

geos::FlowSolverBase::FlowSolverBase	(	const string &	name,
		Group *const	parent
	)

main constructor for Group Objects

Parameters

name	the name of this instantiation of Group in the repository
parent	the parent group of this instantiation of Group

Member Function Documentation

checkSequentialSolutionIncrements()

virtual bool geos::FlowSolverBase::checkSequentialSolutionIncrements ( DomainPartition &  domain ) const

overridevirtual

Check if the solution increments are ok to use.

Parameters

domain

the domain partition

Returns

true if the solution increments are ok to use, false otherwise

Reimplemented from geos::PhysicsSolverBase.

Reimplemented in geos::CompositionalMultiphaseBase.

computeSourceFluxSizeScalingFactor()

void geos::FlowSolverBase::computeSourceFluxSizeScalingFactor	(	real64 const &	time,
		real64 const &	dt,
		DomainPartition &	domain,
		std::map< string, localIndex > const &	bcNameToBcId,
		arrayView1d< globalIndex > const &	bcAllSetsSize
	)		const

For each source flux boundary condition, loop over all the target cells and sum the owned cells.

Parameters

[in]	time	the time at the beginning of the time step
[in]	dt	the time step size
[in]	domain	the domain partition
[in]	bcNameToBcId	the map from the name of the boundary condition to the boundary condition index
[out]	bcAllSetsSize	the total number of owned cells for each source flux boundary condition

findMinMaxElevationInEquilibriumTarget()

void geos::FlowSolverBase::findMinMaxElevationInEquilibriumTarget	(	DomainPartition &	domain,
		std::map< string, localIndex > const &	equilNameToEquilId,
		arrayView1d< real64 > const &	maxElevation,
		arrayView1d< real64 > const &	minElevation
	)		const

For each equilibrium initial condition, loop over all the target cells and compute the min/max elevation.

Parameters

[in]	domain	the domain partition
[in]	equilNameToEquilId	the map from the name of the initial condition to the initial condition index (used in min/maxElevation)
[out]	maxElevation	the max elevation for each initial condition
[out]	minElevation	the min elevation for each initial condition

getMassUnit()

virtual units::Unit geos::FlowSolverBase::getMassUnit ( ) const

inlinevirtual

Returns

The unit in which we evaluate the amount of fluid per element (Mass or Mole).

Reimplemented in geos::CompositionalMultiphaseBase.

Definition at line 148 of file FlowSolverBase.hpp.

prepareStencilWeights()

void geos::FlowSolverBase::prepareStencilWeights

(

DomainPartition &

domain

)

const

Prepare the stencil weights by removing the contribution of the hydraulic aperture before the aperture is updated.

Parameters

[in]

domain

the domain partition

registerDataOnMesh()

virtual void geos::FlowSolverBase::registerDataOnMesh ( Group &  MeshBodies )

overridevirtual

Register wrappers that contain data on the mesh objects.

Parameters

MeshBodies

the group of mesh bodies

Reimplemented from geos::PhysicsSolverBase.

Reimplemented in geos::SinglePhaseHybridFVM, geos::SinglePhaseBase, geos::ReactiveCompositionalMultiphaseOBL, geos::ProppantTransport, geos::CompositionalMultiphaseHybridFVM, and geos::CompositionalMultiphaseBase.

saveAquiferConvergedState()

virtual void geos::FlowSolverBase::saveAquiferConvergedState ( real64 const &  time, real64 const &  dt, DomainPartition &  domain  )

protectedvirtual

Increment the cumulative flux from each aquifer.

Parameters

[in]	time	the time at the beginning of the time step
[in]	dt	the time step size
[in]	domain	the domain partition

For now this function is here because it can be used for both single-phase flow and multiphase flow This may have to be revisited when aquifer BC is implemented for hybrid FVM

Reimplemented in geos::SinglePhaseHybridFVM, and geos::CompositionalMultiphaseHybridFVM.

saveConvergedState()

virtual void geos::FlowSolverBase::saveConvergedState ( ElementSubRegionBase &  subRegion ) const

protectedvirtual

Utility function to save the converged state.

Parameters

[in]

subRegion

the element subRegion

Reimplemented in geos::CompositionalMultiphaseBase, and geos::SinglePhaseBase.

saveSequentialIterationState()

virtual void geos::FlowSolverBase::saveSequentialIterationState ( DomainPartition &  domain )

overridevirtual

Utility function to save the iteration state (useful for sequential simulations)

Parameters

[in]

domain

the domain partition

Reimplemented from geos::PhysicsSolverBase.

Reimplemented in geos::CompositionalMultiphaseBase.

setConstitutiveNamesCallSuper()

virtual void geos::FlowSolverBase::setConstitutiveNamesCallSuper ( ElementSubRegionBase &  subRegion ) const

overrideprotectedvirtual

This function sets constitutive name fields on an ElementSubRegionBase, and calls the base function it overrides.

Parameters

subRegion

The ElementSubRegionBase that will have constitutive names set.

Reimplemented from geos::PhysicsSolverBase.

Reimplemented in geos::SinglePhaseBase.

setKeepVariablesConstantDuringInitStep()

void geos::FlowSolverBase::setKeepVariablesConstantDuringInitStep ( bool const  keepVariablesConstantDuringInitStep )

inline

Utility function to keep the flow variables during a time step (used in poromechanics simulations)

Parameters

[in]

keepVariablesConstantDuringInitStep

flag to tell the solver to freeze its primary variables during a time step @detail This function is meant to be called by a specific task before/after the initialization step

Definition at line 161 of file FlowSolverBase.hpp.

updateStencilWeights()

void geos::FlowSolverBase::updateStencilWeights

(

DomainPartition &

domain

)

const

Update the stencil weights by adding the contribution of the hydraulic aperture after the aperture is updated.

Parameters

[in]

domain

the domain partition

validatePoreVolumes()

virtual void geos::FlowSolverBase::validatePoreVolumes ( DomainPartition const &  domain ) const

protectedvirtual

Helper function to compute/report the elements with small pore volumes.

Parameters

[in]

domain

the domain partition

---

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

• coreComponents/physicsSolvers/fluidFlow/FlowSolverBase.hpp