#include <CompositionalMultiphaseFVM.hpp>

Inheritance diagram for geos::CompositionalMultiphaseFVM:

[legend]

Classes
struct	DBCParameters

struct	valueAndLocation
	Storage for value and element location, used to determine global max + location. More...

struct	viewKeyStruct

Public Types
enum class	ScalingType : integer { Global , Local }
	Solution scaling type. More...

typedef valueAndLocation< real64, globalIndex >	valueAndLocationType

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.

Public Member Functions
	CompositionalMultiphaseFVM (const string &name, Group *const parent)
	main constructor for Group Objects More...

	CompositionalMultiphaseFVM ()=delete
	deleted default constructor

	CompositionalMultiphaseFVM (CompositionalMultiphaseFVM const &)=delete
	deleted copy constructor

	CompositionalMultiphaseFVM (CompositionalMultiphaseFVM &&)=default
	default move constructor

CompositionalMultiphaseFVM &	operator= (CompositionalMultiphaseFVM const &)=delete
	deleted assignment operator

CompositionalMultiphaseFVM &	operator= (CompositionalMultiphaseFVM &&)=delete
	deleted move operator

virtual	~CompositionalMultiphaseFVM () override=default
	default destructor

string	getCatalogName () const override

virtual void	setupDofs (DomainPartition const &domain, DofManager &dofManager) const override
	Populate degree-of-freedom manager with fields relevant to this solver. More...

virtual void	applyBoundaryConditions (real64 const time_n, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) override
	apply boundary condition to system More...

virtual real64	calculateResidualNorm (real64 const &time_n, real64 const &dt, DomainPartition const &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localRhs) override
	calculate the norm of the global system residual More...

virtual real64	scalingForSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution) override
	Function to determine if the solution vector should be scaled back in order to maintain a known constraint. More...

virtual bool	checkSystemSolution (DomainPartition &domain, DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor) override
	Function to check system solution for physical consistency and constraint violation. More...

virtual void	applySystemSolution (DofManager const &dofManager, arrayView1d< real64 const > const &localSolution, real64 const scalingFactor, real64 const dt, DomainPartition &domain) override
	Function to apply the solution vector to the state. More...

virtual void	assembleFluxTerms (real64 const dt, DomainPartition const &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const override
	assembles the flux terms for all cells More...

virtual void	assembleStabilizedFluxTerms (real64 const dt, DomainPartition const &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const override
	assembles the flux terms for all cells with pressure jump stabilization More...

virtual void	updatePhaseMobility (ObjectManagerBase &dataGroup) const override
	Recompute phase mobility from constitutive and primary variables. More...

virtual void	applyAquiferBC (real64 const time, real64 const dt, DofManager const &dofManager, DomainPartition &domain, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const override
	Apply aquifer boundary conditions to the system. More...

Public Member Functions inherited from geos::CompositionalMultiphaseBase
	CompositionalMultiphaseBase (const string &name, Group *const parent)
	main constructor for Group Objects More...

	CompositionalMultiphaseBase ()=delete
	deleted default constructor

	CompositionalMultiphaseBase (CompositionalMultiphaseBase const &)=delete
	deleted copy constructor

	CompositionalMultiphaseBase (CompositionalMultiphaseBase &&)=default
	default move constructor

CompositionalMultiphaseBase &	operator= (CompositionalMultiphaseBase const &)=delete
	deleted assignment operator

CompositionalMultiphaseBase &	operator= (CompositionalMultiphaseBase &&)=delete
	deleted move operator

virtual	~CompositionalMultiphaseBase () override=default
	default destructor

virtual void	registerDataOnMesh (Group &meshBodies) override
	Register wrappers that contain data on the mesh objects. More...

virtual void	implicitStepSetup (real64 const &time_n, real64 const &dt, DomainPartition &domain) override
	function to perform setup for implicit timestep More...

virtual void	assembleSystem (real64 const time_n, real64 const dt, DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) override
	function to assemble the linear system matrix and rhs More...

virtual void	resetStateToBeginningOfStep (DomainPartition &domain) override
	reset state of physics back to the beginning of the step. More...

virtual void	implicitStepComplete (real64 const &time, real64 const &dt, DomainPartition &domain) override
	perform cleanup for implicit timestep More...

void	updateGlobalComponentFraction (ObjectManagerBase &dataGroup) const
	Recompute global component fractions from primary variables (component densities) More...

real64	updatePhaseVolumeFraction (ObjectManagerBase &dataGroup) const
	Recompute phase volume fractions (saturations) from constitutive and primary variables. More...

void	updateFluidModel (ObjectManagerBase &dataGroup) const
	Update all relevant fluid models using current values of pressure and composition. More...

void	updateRelPermModel (ObjectManagerBase &dataGroup) const
	Update all relevant relperm models using current values of phase volume fraction. More...

void	updateCapPressureModel (ObjectManagerBase &dataGroup) const
	Update all relevant capillary pressure models using current values of phase volume fraction. More...

void	updateCompAmount (ElementSubRegionBase &subRegion) const
	Update components mass/moles. More...

void	updateEnergy (ElementSubRegionBase &subRegion) const
	Update energy. More...

void	updateSolidInternalEnergyModel (ObjectManagerBase &dataGroup) const
	Update all relevant solid internal energy models using current values of temperature. More...

real64	updateFluidState (ElementSubRegionBase &subRegion) const

virtual void	saveConvergedState (ElementSubRegionBase &subRegion) const override final
	Utility function to save the converged state. More...

virtual void	saveSequentialIterationState (DomainPartition &domain) override final
	Utility function to save the iteration state (useful for sequential simulations) More...

virtual void	updateState (DomainPartition &domain) override final
	Recompute all dependent quantities from primary variables (including constitutive models) More...

integer	numFluidComponents () const
	Getter for the number of fluid components (species) More...

integer	numFluidPhases () const
	Getter for the number of fluid phases. More...

string	referenceFluidModelName () const
	Getter for the name of the reference fluid model name. More...

virtual units::Unit	getMassUnit () const override

void	assembleAccumulationAndVolumeBalanceTerms (DomainPartition &domain, DofManager const &dofManager, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const
	assembles the accumulation and volume balance terms for all cells More...

virtual void	initializeFluidState (MeshLevel &mesh, arrayView1d< string const > const &regionNames) override
	Initialize all variables from initial conditions. More...

virtual void	initializeThermalState (MeshLevel &mesh, arrayView1d< string const > const &regionNames) override

virtual void	computeHydrostaticEquilibrium (DomainPartition &domain) override
	Compute the hydrostatic equilibrium using the compositions and temperature input tables.

void	applyDirichletBC (real64 const time, real64 const dt, DofManager const &dofManager, DomainPartition &domain, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const
	Function to perform the Application of Dirichlet type BC's. More...

void	applySourceFluxBC (real64 const time, real64 const dt, DofManager const &dofManager, DomainPartition &domain, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const
	Apply source flux boundary conditions to the system. More...

void	keepVariablesConstantDuringInitStep (real64 const time, real64 const dt, DofManager const &dofManager, DomainPartition &domain, CRSMatrixView< real64, globalIndex const > const &localMatrix, arrayView1d< real64 > const &localRhs) const
	Function to fix the initial state during the initialization step in coupled problems. More...

void	chopNegativeDensities (DomainPartition &domain)
	Sets all the negative component densities (if any) to zero. More...

void	chopNegativeDensities (ElementSubRegionBase &subRegion)

virtual real64	setNextDtBasedOnStateChange (real64 const &currentDt, DomainPartition &domain) override
	function to set the next dt based on state change More...

void	computeCFLNumbers (DomainPartition &domain, real64 const &dt, real64 &maxPhaseCFL, real64 &maxCompCFL)

real64	setNextDt (real64 const &currentDt, DomainPartition &domain) override
	function to set the next time step size More...

virtual real64	setNextDtBasedOnCFL (real64 const &currentDt, DomainPartition &domain) override
	function to set the next dt based on state change More...

virtual void	initializePostInitialConditionsPreSubGroups () override
	Called by InitializePostInitialConditions() prior to initializing sub-Groups.

integer	useSimpleAccumulation () const

integer	useTotalMassEquation () const

virtual bool	checkSequentialSolutionIncrements (DomainPartition &domain) const override
	Check if the solution increments are ok to use. More...

Public Member Functions inherited from geos::FlowSolverBase
	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

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

integer &	isThermal ()

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

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)

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

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 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	setNextDtBasedOnNewtonIter (real64 const &currentDt)
	function to set the next time step size based on Newton convergence 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	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...

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 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	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 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	catalogName ()
	name of the solver in the object catalog More...

Static Public Member Functions inherited from geos::FlowSolverBase
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	postInputInitialization () override

virtual void	initializePreSubGroups () override
	Called by Initialize() prior to initializing sub-Groups.

Protected Member Functions inherited from geos::CompositionalMultiphaseBase
void	validateConstitutiveModels (DomainPartition const &domain) const
	Utility function that checks the consistency of the constitutive models. More...

void	initializeAquiferBC (constitutive::ConstitutiveManager const &cm) const
	Initialize the aquifer boundary condition (gravity vector, water phase index) More...

template<typename OBJECT_TYPE >
void	applyFieldValue (real64 const &time_n, real64 const &dt, MeshLevel &mesh, char const logMessage[], string const fieldKey, string const boundaryFieldKey) const
	Utility function that encapsulates the call to FieldSpecificationBase::applyFieldValue in BC application. More...

Protected Member Functions inherited from geos::FlowSolverBase
virtual void	saveAquiferConvergedState (real64 const &time, real64 const &dt, DomainPartition &domain)
	Increment the cumulative flux from each aquifer. 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)

void	initialize (DomainPartition &domain)

void	initializePorosityAndPermeability (MeshLevel &mesh, arrayView1d< string const > const &regionNames)

void	initializeHydraulicAperture (MeshLevel &mesh, const arrayView1d< const string > &regionNames)

void	saveInitialPressureAndTemperature (MeshLevel &mesh, const arrayView1d< const string > &regionNames)

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	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
struct geos::CompositionalMultiphaseFVM::DBCParameters	m_dbcParams

ScalingType	m_scalingType
	Solution scaling type.

Protected Attributes inherited from geos::CompositionalMultiphaseBase
integer	m_numPhases
	the max number of fluid phases

integer	m_numComponents
	the number of fluid components

integer	m_useMass
	flag indicating whether mass or molar formulation should be used

integer	m_hasCapPressure
	flag to determine whether or not to apply capillary pressure

integer	m_hasDiffusion
	flag to determine whether or not to apply diffusion

integer	m_hasDispersion
	flag to determine whether or not to apply dispersion

real64	m_maxCompFracChange
	maximum (absolute) change in a component fraction in a Newton iteration

real64	m_maxRelativePresChange
	maximum (relative) change in pressure in a Newton iteration

real64	m_maxRelativeTempChange
	maximum (relative) change in temperature in a Newton iteration

real64	m_maxRelativeCompDensChange
	maximum (relative) change in component density in a Newton iteration

real64	m_solutionChangeScalingFactor
	damping factor for solution change targets

real64	m_targetRelativePresChange
	target (relative) change in pressure in a time step

real64	m_targetRelativeTempChange
	target (relative) change in temperature in a time step

real64	m_targetPhaseVolFracChange
	target (absolute) change in phase volume fraction in a time step

real64	m_targetRelativeCompDensChange
	target (relative) change in component density in a time step

real64	m_minScalingFactor
	minimum value of the scaling factor obtained by enforcing maxCompFracChange

integer	m_allowCompDensChopping
	flag indicating whether local (cell-wise) chopping of negative compositions is allowed

integer	m_useTotalMassEquation
	flag indicating whether total mass equation is used

integer	m_useSimpleAccumulation
	flag indicating whether simple accumulation form is used

real64	m_minCompDens
	minimum allowed global component density

string	m_referenceFluidModelName
	name of the fluid constitutive model used as a reference for component/phase description

real64	m_sequentialCompDensChange
	maximum (absolute) component density change in a sequential iteration

real64	m_maxSequentialCompDensChange

real64	m_targetFlowCFL
	the targeted CFL for timestep

Protected Attributes inherited from geos::FlowSolverBase
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
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

A compositional multiphase solver using only cell-centered variables works with both TPFA and MPFA

Definition at line 36 of file CompositionalMultiphaseFVM.hpp.

Member Enumeration Documentation

◆ ScalingType

enum geos::CompositionalMultiphaseFVM::ScalingType : integer

strong

Solution scaling type.

Enumerator
Global	Scale the Newton update with a unique scaling factor.
Local	Scale the Newton update locally (modifies the Newton direction)

Definition at line 171 of file CompositionalMultiphaseFVM.hpp.

Constructor & Destructor Documentation

◆ CompositionalMultiphaseFVM()

geos::CompositionalMultiphaseFVM::CompositionalMultiphaseFVM	(	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

◆ applyAquiferBC()

virtual void geos::CompositionalMultiphaseFVM::applyAquiferBC	(	real64 const	time,
		real64 const	dt,
		DofManager const &	dofManager,
		DomainPartition &	domain,
		CRSMatrixView< real64, globalIndex const > const &	localMatrix,
		arrayView1d< real64 > const &	localRhs
	)		const

overridevirtual

Apply aquifer boundary conditions to the system.

Parameters

time	current time
dt	time step
dofManager	degree-of-freedom manager associated with the linear system
domain	the domain
localMatrix	local system matrix
localRhs	local system right-hand side vector

Implements geos::CompositionalMultiphaseBase.

◆ assembleFluxTerms()

virtual void geos::CompositionalMultiphaseFVM::assembleFluxTerms	(	real64 const	dt,
		DomainPartition const &	domain,
		DofManager const &	dofManager,
		CRSMatrixView< real64, globalIndex const > const &	localMatrix,
		arrayView1d< real64 > const &	localRhs
	)		const

overridevirtual

assembles the flux terms for all cells

Parameters

time_n	previous time value
dt	time step
domain	the physical domain object
dofManager	degree-of-freedom manager associated with the linear system
matrix	the system matrix
rhs	the system right-hand side vector

Implements geos::CompositionalMultiphaseBase.

◆ assembleStabilizedFluxTerms()

virtual void geos::CompositionalMultiphaseFVM::assembleStabilizedFluxTerms	(	real64 const	dt,
		DomainPartition const &	domain,
		DofManager const &	dofManager,
		CRSMatrixView< real64, globalIndex const > const &	localMatrix,
		arrayView1d< real64 > const &	localRhs
	)		const

overridevirtual

assembles the flux terms for all cells with pressure jump stabilization

Parameters

dt	time step
domain	the physical domain object
dofManager	degree-of-freedom manager associated with the linear system
matrix	the system matrix
rhs	the system right-hand side vector

Implements geos::CompositionalMultiphaseBase.

◆ catalogName()

static string geos::CompositionalMultiphaseFVM::catalogName ( )

inlinestatic

name of the solver in the object catalog

Returns: string that contains the catalog name to generate a new object through the object catalog.

Definition at line 74 of file CompositionalMultiphaseFVM.hpp.

◆ getCatalogName()

string geos::CompositionalMultiphaseFVM::getCatalogName ( ) const

inlineoverridevirtual

Returns: Get the final class Catalog name

Implements geos::PhysicsSolverBase.

Definition at line 78 of file CompositionalMultiphaseFVM.hpp.

◆ postInputInitialization()

virtual void geos::CompositionalMultiphaseFVM::postInputInitialization ( )

overrideprotectedvirtual

This function provides capability to post process input values prior to any other initialization operations.

Reimplemented from geos::CompositionalMultiphaseBase.

◆ updatePhaseMobility()

virtual void geos::CompositionalMultiphaseFVM::updatePhaseMobility ( ObjectManagerBase & dataGroup ) const

overridevirtual

Recompute phase mobility from constitutive and primary variables.

Parameters

dataGroup the group storing the required field

Implements geos::CompositionalMultiphaseBase.

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

coreComponents/physicsSolvers/fluidFlow/CompositionalMultiphaseFVM.hpp

Classes

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ ScalingType

Constructor & Destructor Documentation

◆ CompositionalMultiphaseFVM()

Member Function Documentation

◆ applyAquiferBC()

◆ assembleFluxTerms()

◆ assembleStabilizedFluxTerms()

◆ catalogName()

◆ getCatalogName()

◆ postInputInitialization()

◆ updatePhaseMobility()