https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Private Types | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
GeochemistryTimeDependentReactor Class Reference

Class that controls the time-dependent (but not space-dependent) geochemistry reactions. More...

#include <GeochemistryTimeDependentReactor.h>

Inheritance diagram for GeochemistryTimeDependentReactor:
[legend]

Public Types

enum  TEST_TYPE
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 GeochemistryTimeDependentReactor (const InputParameters &parameters)
 
virtual void initialize () override
 
virtual void finalize () override
 
virtual void initialSetup () override
 
virtual void execute () override
 
virtual const GeochemicalSystemgetGeochemicalSystem (dof_id_type node_id) const override
 
virtual const std::stringstream & getSolverOutput (dof_id_type node_id) const override
 
virtual unsigned getSolverIterations (dof_id_type node_id) const override
 
virtual Real getSolverResidual (dof_id_type node_id) const override
 
virtual const DenseVector< Real > & getMoleAdditions (dof_id_type node_id) const override
 
virtual Real getMolesDumped (dof_id_type node_id, const std::string &species) const override
 
virtual void threadJoin (const UserObject &) override
 
const PertinentGeochemicalSystemgetPertinentGeochemicalSystem () const
 returns a reference to the PertinentGeochemicalSystem used to creat the ModelGeochemicalDatabase
 
virtual void subdomainSetup () override
 
bool isUniqueNodeExecute ()
 
virtual Real spatialValue (const Point &) const
 
virtual const std::vector< Point > spatialPoints () const
 
void setPrimaryThreadCopy (UserObject *primary)
 
UserObjectprimaryThreadCopy ()
 
SubProblemgetSubProblem () const
 
bool shouldDuplicateInitialExecution () const
 
void gatherSum (T &value)
 
void gatherMax (T &value)
 
void gatherMin (T &value)
 
void gatherProxyValueMax (T1 &proxy, T2 &value)
 
void gatherProxyValueMin (T1 &proxy, T2 &value)
 
std::set< UserObjectName > getDependObjects () const
 
const std::set< std::string > & getRequestedItems () override
 
const std::set< std::string > & getSuppliedItems () override
 
unsigned int systemNumber () const
 
virtual bool needThreadedCopy () const
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
bool isKokkosObject () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
const std::string & name () const
 
std::string typeAndName () const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
MooseObjectName uniqueName () const
 
const InputParametersparameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const TgetParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const TqueryParam (const std::string &name) const
 
const TgetRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &name) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
std::string messagePrefix (const bool hit_prefix=true) const
 
std::string errorPrefix (const std::string &) const
 
void mooseError (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseDeprecatedNoTrace (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
void callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void timestepSetup ()
 
virtual void jacobianSetup ()
 
virtual void residualSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnumgetExecuteOnEnum () const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const TgetUserObject (const std::string &param_name, bool is_dependency=true) const
 
const TgetUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const std::vector< SubdomainName > & blocks () const
 
unsigned int numBlocks () const
 
virtual const std::set< SubdomainID > & blockIDs () const
 
unsigned int blocksMaxDimension () const
 
bool hasBlocks (const SubdomainName &name) const
 
bool hasBlocks (const std::vector< SubdomainName > &names) const
 
bool hasBlocks (const std::set< SubdomainName > &names) const
 
bool hasBlocks (SubdomainID id) const
 
bool hasBlocks (const std::vector< SubdomainID > &ids) const
 
bool hasBlocks (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::vector< SubdomainID > &ids) const
 
bool hasBlockMaterialProperty (const std::string &prop_name)
 
const std::set< SubdomainID > & meshBlockIDs () const
 
virtual bool blockRestricted () const
 
virtual void checkVariable (const MooseVariableFieldBase &variable) const
 
virtual const std::set< BoundaryID > & boundaryIDs () const
 
const std::vector< BoundaryName > & boundaryNames () const
 
unsigned int numBoundaryIDs () const
 
bool hasBoundary (const BoundaryName &name) const
 
bool hasBoundary (const std::vector< BoundaryName > &names) const
 
bool hasBoundary (const BoundaryID &id) const
 
bool hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const
 
bool hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const
 
bool hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const
 
bool hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const
 
bool isBoundarySubset (const std::set< BoundaryID > &ids) const
 
bool isBoundarySubset (const std::vector< BoundaryID > &ids) const
 
bool hasBoundaryMaterialProperty (const std::string &prop_name) const
 
virtual bool boundaryRestricted () const
 
const std::set< BoundaryID > & meshBoundaryIDs () const
 
virtual bool checkVariableBoundaryIntegrity () const
 
virtual const VariableValuecoupledValueByName (const std::string &var_name)
 
virtual const ArrayVariableValuecoupledArrayValueByName (const std::string &var_name)
 
const std::unordered_map< std::string, std::vector< MooseVariableFieldBase * > > & getCoupledVars () const
 
const std::vector< MooseVariableFieldBase * > & getCoupledMooseVars () const
 
const std::vector< MooseVariable * > & getCoupledStandardMooseVars () const
 
const std::vector< VectorMooseVariable * > & getCoupledVectorMooseVars () const
 
const std::vector< ArrayMooseVariable * > & getCoupledArrayMooseVars () const
 
void addFEVariableCoupleableVectorTag (TagID tag)
 
void addFEVariableCoupleableMatrixTag (TagID tag)
 
std::set< TagID > & getFEVariableCoupleableVectorTags ()
 
const std::set< TagID > & getFEVariableCoupleableVectorTags () const
 
std::set< TagID > & getFEVariableCoupleableMatrixTags ()
 
const std::set< TagID > & getFEVariableCoupleableMatrixTags () const
 
auto & getWritableCoupledVariables () const
 
bool hasWritableCoupledVariables () const
 
const ADVariableValuegetADDefaultValue (const std::string &var_name) const
 
const ADVectorVariableValuegetADDefaultVectorValue (const std::string &var_name) const
 
const ADVariableGradientgetADDefaultGradient () const
 
const ADVectorVariableGradientgetADDefaultVectorGradient () const
 
const ADVariableSecondgetADDefaultSecond () const
 
const ADVectorVariableCurlgetADDefaultCurl () const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorTagVectorGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorTagVectorGradientsByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorDofValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorTagVectorGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorTagVectorGradients (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorNodalValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorTagVectorDofValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledValues (const std::string &var_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradients (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValues (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValues (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientsOld (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableGradient kokkosCoupledGradientsOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDofValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledDots (const std::string &var_name) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosCoupledNodalDots (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDofValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorValues (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorGradients (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalValues (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDofValues (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDofValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorGradientsOld (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDofValuesOld (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDofValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableGradient kokkosCoupledVectorGradientsOlder (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDofValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorDots (const std::string &var_name) const
 
Moose::Kokkos::VectorVariableValue kokkosCoupledVectorNodalDots (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradientByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradientsByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValueByName (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValuesByName (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradient (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledVectorTagGradients (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagNodalValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValue (const std::string &var_name, const std::string &tag_param_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledVectorTagDofValues (const std::string &var_name, const std::string &tag_param_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValues (const std::string &var_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradients (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValues (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValue (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValues (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValuesOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientsOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValuesOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValueOld (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValuesOld (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableGradient kokkosADCoupledGradientsOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValueOlder (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDofValuesOlder (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledDots (const std::string &var_name) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::ADVariableValue kokkosADCoupledNodalDots (const std::string &var_name) const
 
Moose::Kokkos::Scalar< const Real > kokkosCoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
Moose::Kokkos::VariableValue kokkosZeroValue () const
 
Moose::Kokkos::VariableGradient kokkosZeroGradient () const
 
Moose::Kokkos::VariableValue kokkosZeroNodalValue () const
 
const std::vector< MooseVariableScalar * > & getCoupledMooseScalarVars ()
 
const std::set< TagID > & getScalarVariableCoupleableVectorTags () const
 
const std::set< TagID > & getScalarVariableCoupleableMatrixTags () const
 
const std::set< MooseVariableFieldBase * > & getMooseVariableDependencies () const
 
std::set< MooseVariableFieldBase * > checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})
 
virtual std::set< MooseVariableFieldBase * > checkVariables (const libMesh::Node &node, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
std::set< MooseVariableFieldBase * > checkVariables (const libMesh::Elem &element, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
void addMooseVariableDependency (MooseVariableFieldBase *var)
 
void addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
 
bool isImplicit ()
 
Moose::StateArg determineState () const
 
void setRandomResetFrequency (ExecFlagType exec_flag)
 
unsigned long getRandomLong () const
 
Real getRandomReal () const
 
unsigned int getSeed (std::size_t id)
 
unsigned int getMasterSeed () const
 
bool isNodal () const
 
ExecFlagType getResetOnTime () const
 
void setRandomDataPointer (RandomData *random_data)
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
const FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
Moose::Kokkos::Function getKokkosFunction (const std::string &name) const
 
const TgetKokkosFunction (const std::string &name) const
 
Moose::Kokkos::Function getKokkosFunctionByName (const FunctionName &name) const
 
const TgetKokkosFunctionByName (const FunctionName &name) const
 
bool hasKokkosFunction (const std::string &param_name) const
 
bool hasKokkosFunctionByName (const FunctionName &name) const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
 
const VectorPostprocessorName & getVectorPostprocessorName (const std::string &param_name) const
 
TgetSampler (const std::string &name)
 
SamplergetSampler (const std::string &name)
 
TgetSamplerByName (const SamplerName &name)
 
SamplergetSamplerByName (const SamplerName &name)
 
virtual void meshChanged ()
 
virtual void meshDisplaced ()
 
PerfGraphperfGraph ()
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

static InputParameters sharedParams ()
 params that are shared with AddTimeDependentReactionSolverAction
 
static InputParameters validParams ()
 
static void callMooseError (MooseApp *const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node *node, const bool show_trace=true)
 
static void sort (typename std::vector< T > &vector)
 
static void sortDFS (typename std::vector< T > &vector)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header, NameFunc &&name_func)
 
static void cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)
 
static bool restricted (const std::set< BoundaryID > &ids)
 

Public Attributes

 ALL
 
 ANY
 
 usingCombinedWarningSolutionWarnings
 
const ConsoleStream _console
 

Static Public Attributes

static const std::string type_param
 
static const std::string name_param
 
static const std::string unique_name_param
 
static const std::string app_param
 
static const std::string moose_base_param
 
static const std::string kokkos_object_param
 
static constexpr auto SYSTEM
 
static constexpr auto NAME
 

Protected Member Functions

Real newTemperature (const DenseVector< Real > &mole_additions) const
 Based on _temperature[0], mole_additions, the current mass of the system and the current temperature of the system, work out the resulting temperature.
 
void preSolveDump ()
 Activate the special "dump" mode prior to solving the geochemical system:
 
void preSolveFlush ()
 Activate the special "flush" mode prior to solving the geochemical system:
 
void postSolveFlowThrough ()
 Activate the special "flow-through" mode after solving the geochemical system:
 
void removeCurrentSpecies ()
 Alter _mole_additions so that it will represent the situation in which all current species are removed and replaced with source_species.
 
void postSolveExchanger ()
 This is relevant for mode=4 simulations (heat-exchanger simulations).
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &name) const override
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override
 
virtual void addUserObjectDependencyHelper (const UserObjectBase &uo) const override
 
void addReporterDependencyHelper (const ReporterName &reporter_name) override
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
const ReporterContextBasegetReporterContextBaseByName (const ReporterName &reporter_name) const
 
const ReporterNamegetReporterName (const std::string &param_name) const
 
TdeclareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< TdeclareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const TgetRestartableData (const std::string &data_name) const
 
TdeclareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
TdeclareRecoverableData (const std::string &data_name, Args &&... args)
 
TdeclareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
TdeclareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
const TgetMeshProperty (const std::string &data_name, const std::string &prefix)
 
const TgetMeshProperty (const std::string &data_name)
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name) const
 
bool hasMeshProperty (const std::string &data_name) const
 
std::string meshPropertyName (const std::string &data_name) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const
 
std::string timedSectionName (const std::string &section_name) const
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
void initializeKokkosBlockRestrictable ()
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockElements () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockNodes () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBlockSides () const
 
KOKKOS_FUNCTION ContiguousElementID kokkosBlockElementID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION ContiguousElementID kokkosBlockNodeID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION auto kokkosBlockElementSideID (Moose::Kokkos::ThreadID tid) const
 
bool hasBoundaryMaterialPropertyHelper (const std::string &prop_name) const
 
KOKKOS_FUNCTION dof_id_type numKokkosBoundaryNodes () const
 
KOKKOS_FUNCTION dof_id_type numKokkosBoundarySides () const
 
KOKKOS_FUNCTION ContiguousNodeID kokkosBoundaryNodeID (Moose::Kokkos::ThreadID tid) const
 
KOKKOS_FUNCTION auto kokkosBoundaryElementSideID (Moose::Kokkos::ThreadID tid) const
 
virtual void coupledCallback (const std::string &, bool) const
 
virtual bool isCoupled (const std::string &var_name, unsigned int i=0) const
 
virtual bool isCoupledConstant (const std::string &var_name) const
 
unsigned int coupledComponents (const std::string &var_name) const
 
VariableName coupledName (const std::string &var_name, unsigned int comp=0) const
 
std::vector< VariableName > coupledNames (const std::string &var_name) const
 
virtual unsigned int coupled (const std::string &var_name, unsigned int comp=0) const
 
std::vector< unsigned intcoupledIndices (const std::string &var_name) const
 
virtual const VariableValuecoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledValues (const std::string &var_name) const
 
std::vector< const VectorVariableValue * > coupledVectorValues (const std::string &var_name) const
 
const GenericVariableValue< is_ad > & coupledGenericValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericValue (const std::string &var_name, unsigned int comp) const
 
const GenericVectorVariableValue< is_ad > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVectorVariableValue< false > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
 
const GenericVectorVariableValue< true > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
 
std::vector< const GenericVariableValue< is_ad > * > coupledGenericValues (const std::string &var_name) const
 
std::vector< const GenericVariableValue< false > * > coupledGenericValues (const std::string &var_name) const
 
std::vector< const GenericVariableValue< true > * > coupledGenericValues (const std::string &var_name) const
 
const GenericVariableValue< is_ad > & coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< is_ad > & coupledGenericDot (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< is_ad > & coupledGenericDotDot (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
 
virtual const VariableValuecoupledValueLower (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue * > adCoupledValues (const std::string &var_name) const
 
const ADVariableValueadCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVectorVariableValue * > adCoupledVectorValues (const std::string &var_name) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableValue * > coupledVectorTagValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue * > coupledVectorTagValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const ArrayVariableValue * > coupledVectorTagArrayValues (const std::string &var_names, TagID tag) const
 
std::vector< const ArrayVariableValue * > coupledVectorTagArrayValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableGradientcoupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableGradientcoupledVectorTagGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableGradient * > coupledVectorTagGradients (const std::string &var_names, TagID tag) const
 
std::vector< const VariableGradient * > coupledVectorTagGradients (const std::string &var_names, const std::string &tag_name) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const ArrayVariableGradient * > coupledVectorTagArrayGradients (const std::string &var_names, TagID tag) const
 
std::vector< const ArrayVariableGradient * > coupledVectorTagArrayGradients (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const ArrayVariableValuecoupledVectorTagArrayDofValue (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledVectorTagDofValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue * > coupledVectorTagDofValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableValuecoupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledMatrixTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableValue * > coupledMatrixTagValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue * > coupledMatrixTagValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VectorVariableValuecoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ArrayVariableValue * > coupledArrayValues (const std::string &var_name) const
 
MooseWritableVariablewritableVariable (const std::string &var_name, unsigned int comp=0)
 
virtual VariableValuewritableCoupledValue (const std::string &var_name, unsigned int comp=0)
 
void checkWritableVar (MooseWritableVariable *var)
 
virtual const VariableValuecoupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledValuesOld (const std::string &var_name) const
 
std::vector< const VectorVariableValue * > coupledVectorValuesOld (const std::string &var_name) const
 
virtual const VariableValuecoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledValuesOlder (const std::string &var_name) const
 
virtual const VariableValuecoupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient * > coupledGradients (const std::string &var_name) const
 
const ADVariableGradientadCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableGradientadCoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableGradient * > adCoupledGradients (const std::string &var_name) const
 
const GenericVariableGradient< is_ad > & coupledGenericGradient (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableGradient< false > & coupledGenericGradient (const std::string &var_name, unsigned int comp) const
 
const GenericVariableGradient< true > & coupledGenericGradient (const std::string &var_name, unsigned int comp) const
 
std::vector< const GenericVariableGradient< is_ad > * > coupledGenericGradients (const std::string &var_name) const
 
std::vector< const GenericVariableGradient< false > * > coupledGenericGradients (const std::string &var_name) const
 
std::vector< const GenericVariableGradient< true > * > coupledGenericGradients (const std::string &var_name) const
 
const ADVectorVariableGradientadCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableSecondadCoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableSecondadCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient * > coupledGradientsOld (const std::string &var_name) const
 
virtual const VariableGradientcoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableCurladCoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDiv (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledDots (const std::string &var_name) const
 
virtual const VariableValuecoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue * > adCoupledDots (const std::string &var_name) const
 
const ADVariableValueadCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const Moose::ADType< T >::typeadCoupledNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
const TcoupledNodalDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledAllDofValues (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledAllDofValuesOld (const std::string &var_name) const
 
virtual const VariableValuecoupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue * > coupledAllDofValuesOlder (const std::string &var_name) const
 
virtual const ArrayVariableValuecoupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableValueadCoupledDofValues (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadZeroValue () const
 
const ADVariableGradientadZeroGradient () const
 
const ADVariableSecondadZeroSecond () const
 
const GenericVariableValue< is_ad > & genericZeroValue ()
 
const GenericVariableValue< false > & genericZeroValue ()
 
const GenericVariableValue< true > & genericZeroValue ()
 
const GenericVariableGradient< is_ad > & genericZeroGradient ()
 
const GenericVariableGradient< false > & genericZeroGradient ()
 
const GenericVariableGradient< true > & genericZeroGradient ()
 
const GenericVariableSecond< is_ad > & genericZeroSecond ()
 
const GenericVariableSecond< false > & genericZeroSecond ()
 
const GenericVariableSecond< true > & genericZeroSecond ()
 
bool checkVar (const std::string &var_name, unsigned int comp=0, unsigned int comp_bound=0) const
 
const MooseVariableFieldBasegetFEVar (const std::string &var_name, unsigned int comp) const
 
const MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp) const
 
MooseVariableFieldBasegetFieldVar (const std::string &var_name, unsigned int comp)
 
std::vector< const MooseVariableFieldBase * > getFieldVars (const std::string &var_name) const
 
const TgetVarHelper (const std::string &var_name, unsigned int comp) const
 
TgetVarHelper (const std::string &var_name, unsigned int comp)
 
std::vector< const T * > getVarsHelper (const std::string &var_name) const
 
MooseVariablegetVar (const std::string &var_name, unsigned int comp)
 
const MooseVariablegetVar (const std::string &var_name, unsigned int comp) const
 
VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp)
 
const VectorMooseVariablegetVectorVar (const std::string &var_name, unsigned int comp) const
 
ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp)
 
const ArrayMooseVariablegetArrayVar (const std::string &var_name, unsigned int comp) const
 
void validateExecutionerType (const std::string &name, const std::string &fn_name) const
 
std::vector< TcoupledVectorHelper (const std::string &var_name, const Func &func) const
 
bool isCoupledScalar (const std::string &var_name, unsigned int i=0) const
 
unsigned int coupledScalarComponents (const std::string &var_name) const
 
unsigned int coupledScalar (const std::string &var_name, unsigned int comp=0) const
 
libMesh::Order coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< is_ad > & coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (const std::string &var_name, unsigned int comp) const
 
const TgetReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValue (const std::string &param_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 
bool hasReporterValueByName (const ReporterName &reporter_name) const
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 

Protected Attributes

const VariableValue_temperature
 Temperature specified by user.
 
const VariableValue_cold_temperature
 Cold temperature specified by user, which is used only when mode==4.
 
const unsigned _heating_increments
 If mode=4 then temperature is ramped from _cold_temperature to _temperature in heating_increments increments.
 
Real _new_temperature
 Temperature at which the solution is required.
 
Real _previous_temperature
 Temperature at which the _egs was last made consistent.
 
GeochemicalSystem _egs
 The equilibrium geochemical system that holds all the molalities, activities, etc.
 
GeochemicalSolver _solver
 The solver.
 
const unsigned _num_kin
 Number of kinetic species.
 
const Real _close_system_at_time
 Defines the time at which to close the system.
 
bool _closed_system
 Whether the system has been closed.
 
const std::vector< std::string > _source_species_names
 Names of the source species.
 
const unsigned _num_source_species
 Number of source species.
 
const std::vector< const VariableValue * > _source_species_rates
 Rates of the source species.
 
const std::vector< std::string > _remove_fixed_activity_name
 Names of species to remove the fixed activity or fugacity constraint from.
 
const std::vector< Real > _remove_fixed_activity_time
 Times at which to remove the fixed activity or fugacity from the species in _remove_fixed_activity_name.
 
const unsigned _num_removed_fixed
 Number of elements in the vector _remove_fixed_activity_name;.
 
std::vector< bool > _removed_fixed_activity
 Whether the activity or activity constraint has been removfed.
 
const std::vector< std::string > _controlled_activity_species_names
 Names of the species with controlled activity or fugacity.
 
const unsigned _num_controlled_activity
 Number of species with controlled activity or fugacity.
 
const std::vector< const VariableValue * > _controlled_activity_species_values
 Activity or fugacity of the species with controlled activity or fugacity.
 
DenseVector< Real > _mole_additions
 Moles of each basis species added at the current timestep, along with kinetic rates.
 
DenseMatrix< Real > _dmole_additions
 Derivative of moles_added.
 
const VariableValue_mode
 Mode of the system (flush, flow-through, etc)
 
std::unordered_map< std::string, Real > _minerals_dumped
 Moles of mineral removed by dump and flow-through.
 
const unsigned _ramp_subsequent
 the ramp_max_ionic_strength to use during time-stepping
 
const unsigned _num_my_nodes
 Number of nodes handled by this processor (will need to be made un-const when mesh adaptivity is handled)
 
ModelGeochemicalDatabase _mgd
 my copy of the underlying ModelGeochemicalDatabase
 
const PertinentGeochemicalSystem_pgs
 Reference to the original PertinentGeochemicalSystem used to create the ModelGeochemicalDatabase.
 
const unsigned _num_basis
 number of basis species
 
const unsigned _num_eqm
 number of equilibrium species
 
const Real _initial_max_ionic_str
 Initial value of maximum ionic strength.
 
GeochemistryIonicStrength _is
 The ionic strength calculator.
 
GeochemistryActivityCoefficientsDebyeHuckel _gac
 The activity calculator.
 
const unsigned _max_swaps_allowed
 Maximum number of swaps allowed during a single solve.
 
GeochemistrySpeciesSwapper _swapper
 The species swapper.
 
const Real _small_molality
 A small value of molality.
 
std::vector< std::stringstream > _solver_output
 The solver output at each node.
 
std::vector< unsigned > _tot_iter
 Number of iterations used by the solver at each node.
 
std::vector< Real > _abs_residual
 L1norm of the solver residual at each node.
 
MooseMesh_mesh
 
const unsigned int _qp
 
const Node *const & _current_node
 
const bool & _unique_node_execute
 
const Moose::CoordinateSystemType_coord_sys
 
const THREAD_ID _tid
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
SystemBase_sys
 
Assembly_assembly
 
const bool _duplicate_initial_execution
 
std::set< std::string > _depend_uo
 
const bool & _enabled
 
MooseApp_app
 
Factory_factory
 
ActionFactory_action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters_pars
 
const ExecFlagEnum_execute_enum
 
const ExecFlagType_current_execute_flag
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
FEProblemBase_mdi_feproblem
 
MooseApp_pg_moose_app
 
const std::string _prefix
 
const MaterialData_blk_material_data
 
const InputParameters_c_parameters
 
const std::string & _c_name
 
const std::string & _c_type
 
FEProblemBase_c_fe_problem
 
const SystemBase *const _c_sys
 
std::unordered_map< std::string, std::vector< MooseVariableFieldBase * > > _coupled_vars
 
std::vector< MooseVariableFieldBase * > _coupled_moose_vars
 
std::vector< MooseVariable * > _coupled_standard_moose_vars
 
std::vector< VectorMooseVariable * > _coupled_vector_moose_vars
 
std::vector< ArrayMooseVariable * > _coupled_array_moose_vars
 
std::vector< MooseVariableField< Real > * > _coupled_fv_moose_vars
 
const std::unordered_map< std::string, std::string > & _new_to_deprecated_coupled_vars
 
bool _c_nodal
 
bool _c_is_implicit
 
const bool _c_allow_element_to_nodal_coupling
 
THREAD_ID _c_tid
 
std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > > _default_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADReal > > > _ad_default_value
 
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > > _default_vector_value
 
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > > _default_array_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > > _ad_default_vector_value
 
VariableValue _default_value_zero
 
VariableGradient _default_gradient
 
MooseArray< ADRealVectorValue_ad_default_gradient
 
MooseArray< ADRealTensorValue_ad_default_vector_gradient
 
VariableSecond _default_second
 
MooseArray< ADRealTensorValue_ad_default_second
 
MooseArray< ADRealVectorValue_ad_default_curl
 
const VariableValue_zero
 
const VariablePhiValue_phi_zero
 
const MooseArray< ADReal > & _ad_zero
 
const VariableGradient_grad_zero
 
const MooseArray< ADRealVectorValue > & _ad_grad_zero
 
const VariablePhiGradient_grad_phi_zero
 
const VariableSecond_second_zero
 
const MooseArray< ADRealTensorValue > & _ad_second_zero
 
const VariablePhiSecond_second_phi_zero
 
const VectorVariableValue_vector_zero
 
const VectorVariableCurl_vector_curl_zero
 
VectorVariableValue _default_vector_value_zero
 
VectorVariableGradient _default_vector_gradient
 
VectorVariableCurl _default_vector_curl
 
VectorVariableDivergence _default_div
 
ArrayVariableValue _default_array_value_zero
 
ArrayVariableGradient _default_array_gradient
 
bool _coupleable_neighbor
 
FEProblemBase_sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real & _real_zero
 
const VariableValue_scalar_zero
 
const Point & _point_zero
 
const InputParameters_ti_params
 
FEProblemBase_ti_feproblem
 
bool _is_implicit
 
Real & _t
 
const Real & _t_old
 
int_t_step
 
Real & _dt
 
Real & _dt_old
 
bool _is_transient
 
const Parallel::Communicator & _communicator
 

Private Types

enum  FuncAge
 
enum  VarType
 

Private Member Functions

const Moose::FunctionBasegetKokkosFunctionByNameHelper (const FunctionName &name) const
 
const UserObjectBasegetUserObjectFromFEProblem (const UserObjectName &object_name, const THREAD_ID tid=0) const
 
const TcastUserObject (const UserObjectBase &uo_base, const std::string &param_name="") const
 
void mooseObjectError (const std::string &param_name, std::stringstream &oss) const
 
const std::string & userObjectType (const UserObjectBase &uo) const
 
const std::string & userObjectName (const UserObjectBase &uo) const
 
const PostprocessorName & getPostprocessorNameInternal (const std::string &param_name, const unsigned int index, const bool allow_default_value=true) const
 
bool isDefaultPostprocessorValueByName (const PostprocessorName &name) const
 
PostprocessorValue getDefaultPostprocessorValueByName (const PostprocessorName &name) const
 
void checkParam (const std::string &param_name, const unsigned int index=std::numeric_limits< unsigned int >::max()) const
 
bool postprocessorsAdded () const
 
const VectorPostprocessorValuegetVectorPostprocessorByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name, bool broadcast, std::size_t t_index) const
 
const VectorPostprocessorContext< VectorPostprocessorValue > & getVectorPostprocessorContextByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool vectorPostprocessorsAdded () const
 
bool reportersAdded () const
 
void possiblyCheckHasReporter (const ReporterName &reporter_name, const std::string &param_name="") const
 
RestartableDataValueregisterRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const
 
void registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)
 
RestartableData< T > & declareRestartableDataHelper (const std::string &data_name, void *context, Args &&... args) const
 
virtual std::string meshPropertyPrefix (const std::string &data_name) const
 
const RestartableDataValuegetMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const
 
void mooseErrorInternal (Args &&... args) const
 
void initializeBoundaryRestrictable ()
 
void initializeKokkosBoundaryRestrictable ()
 
const OutputTools< T >::VariableValuevectorTagValueHelper (const std::string &var_names, TagID tag, unsigned int index=0) const
 
const OutputTools< T >::VariableValuevectorTagValueHelper (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const OutputTools< T >::VariableValuevectorTagDofValueHelper (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const OutputTools< T >::VariableValuevectorTagDofValueHelper (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
void requestStates (const std::string &var_name, const TagName &tag_name, const unsigned int comp)
 
void checkFuncType (const std::string var_name, VarType t, FuncAge age) const
 
const VariableValuegetDefaultValue (const std::string &var_name, unsigned int comp) const
 
const VectorVariableValuegetDefaultVectorValue (const std::string &var_name) const
 
const ArrayVariableValuegetDefaultArrayValue (const std::string &var_name) const
 
const TgetDefaultNodalValue (const std::string &var_name, unsigned int comp=0) const
 
const RealVectorValue & getDefaultNodalValue (const std::string &var_name, unsigned int) const
 
const RealEigenVector & getDefaultNodalValue (const std::string &var_name, unsigned int) const
 
const Moose::Functor< T > & getDefaultFunctor (const std::string &var_name) const
 
Moose::Kokkos::Variable kokkosCoupledVectorTagVariable (const std::string &var_name, const std::string &tag_name, unsigned int comp) const
 
Moose::Kokkos::Variable kokkosCoupledVectorTagVariables (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::Variable kokkosCoupledVectorTagVectorVariable (const std::string &var_name, const std::string &tag_name, unsigned int comp) const
 
Moose::Kokkos::Variable kokkosCoupledVectorTagVectorVariables (const std::string &var_name, const std::string &tag_name) const
 
Moose::Kokkos::Variable kokkosZeroVariable () const
 
std::set< MooseVariableFieldBase * > checkVariablesHelper (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
const PostprocessorValuegetPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const
 
const PostprocessorValuegetPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const
 
const PostprocessorValuegetPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const
 
const PostprocessorValuegetPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const
 
void possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 

Static Private Member Functions

static const hit::Node * getHitNode (const InputParameters &params)
 
static std::string messagePrefix (const InputParameters &params, const bool hit_prefix)
 

Private Attributes

UserObject_primary_thread_copy
 
std::set< std::string > _supplied_uo
 
const ParallelParamObject_parent
 
const MooseBase_si_moose_base
 
const FEProblemBase_si_problem
 
ExecFlagEnum _empty_execute_enum
 
const MooseObject_fni_object
 
const InputParameters_fni_params
 
FEProblemBase_fni_feproblem
 
const THREAD_ID _fni_tid
 
const MooseObject_uoi_moose_object
 
const FEProblemBase_uoi_feproblem
 
const THREAD_ID _uoi_tid
 
const MooseObject_ppi_moose_object
 
const InputParameters_ppi_params
 
const FEProblemBase_ppi_feproblem
 
std::map< PostprocessorName, std::unique_ptr< PostprocessorValue > > _default_values
 
const bool _broadcast_by_default
 
const MooseObject_vpi_moose_object
 
const FEProblemBase_vpi_feproblem
 
const THREAD_ID _vpi_tid
 
const InputParameters_ri_params
 
FEProblemBase_ri_fe_problem_base
 
const ReporterData_ri_reporter_data
 
const MooseObject_ri_moose_object
 
const InputParameters_dni_params
 
FEProblemBase_dni_feproblem
 
const MooseObject *const _dni_moose_object_ptr
 
const InputParameters_si_params
 
FEProblemBase_si_feproblem
 
THREAD_ID _si_tid
 
const RestartableDataMapName _metaname
 
std::string _restartable_name
 
MooseApp_meta_data_app
 
const MooseObject *const _meta_data_object
 
std::set< SubdomainID_blk_ids
 
std::vector< SubdomainID_vec_ids
 
std::vector< SubdomainName > _blocks
 
const bool _blk_dual_restrictable
 
FEProblemBase_blk_feproblem
 
const MooseMesh_blk_mesh
 
const std::set< BoundaryID_empty_boundary_ids
 
const std::set< BoundaryID > & _boundary_ids
 
THREAD_ID _blk_tid
 
const std::string & _blk_name
 
unsigned int _blk_dim
 
const MooseObject_moose_object
 
Moose::Kokkos::Array< ContiguousElementID_kokkos_element_ids
 
Moose::Kokkos::Array< ContiguousNodeID_kokkos_node_ids
 
Moose::Kokkos::Array< Moose::Kokkos::Pair< ContiguousElementID, unsigned int > > _kokkos_element_side_ids
 
FEProblemBase_bnd_feproblem
 
MooseMesh_bnd_mesh
 
std::set< BoundaryID_bnd_ids
 
std::vector< BoundaryName > _boundary_names
 
const bool _bnd_dual_restrictable
 
const std::set< SubdomainID_empty_block_ids
 
const std::set< SubdomainID > & _block_ids
 
THREAD_ID _bnd_tid
 
const MaterialData_bnd_material_data
 
bool _bnd_nodal
 
unsigned int _coupleable_max_qps
 
std::unordered_map< std::string, std::vector< unsigned int > > _optional_var_index
 
std::unordered_map< std::string, std::vector< MooseVariableScalar * > > _c_coupled_scalar_vars
 
std::set< TagID_fe_coupleable_vector_tags
 
std::set< TagID_fe_coupleable_matrix_tags
 
const bool _is_fv
 
const MooseObject *const _obj
 
const std::set< std::string > _older_state_tags
 
std::vector< std::set< MooseWritableVariable * > > _writable_coupled_variables
 
const Real _kokkos_default_value_zero
 
const InputParameters_sc_parameters
 
const std::string & _sc_name
 
const bool _sc_is_implicit
 
std::unordered_map< std::string, std::vector< MooseVariableScalar * > > _coupled_scalar_vars
 
std::unordered_map< std::string, std::unique_ptr< ADVariableValue > > _dual_default_value
 
std::vector< MooseVariableScalar * > _coupled_moose_scalar_vars
 
std::unordered_map< std::string, std::vector< MooseVariableFieldBase * > > _sc_coupled_vars
 
std::set< TagID_sc_coupleable_vector_tags
 
std::set< TagID_sc_coupleable_matrix_tags
 
std::set< MooseVariableFieldBase * > _moose_variable_dependencies
 
std::vector< libMesh::dof_id_type_dof_indices
 
const std::string _ti_name
 
RandomData_random_data
 
MooseRandom_generator
 
FEProblemBase_ri_problem
 
const std::string _ri_name
 
unsigned int _master_seed
 
bool _is_nodal
 
ExecFlagType _reset_on
 
const Node *const & _curr_node
 
const Elem *const & _curr_element
 

Detailed Description

Class that controls the time-dependent (but not space-dependent) geochemistry reactions.

Definition at line 18 of file GeochemistryTimeDependentReactor.h.

Constructor & Destructor Documentation

◆ GeochemistryTimeDependentReactor()

GeochemistryTimeDependentReactor::GeochemistryTimeDependentReactor ( const InputParameters parameters)

Definition at line 139 of file GeochemistryTimeDependentReactor.C.

142 _temperature(coupledValue("temperature")),
143 _cold_temperature(coupledValue("cold_temperature")),
144 _heating_increments(getParam<unsigned>("heating_increments")),
145 _new_temperature(getParam<Real>("initial_temperature")),
146 _previous_temperature(getParam<Real>("initial_temperature")),
147 _egs(_mgd,
148 _gac,
149 _is,
150 _swapper,
151 getParam<std::vector<std::string>>("swap_out_of_basis"),
152 getParam<std::vector<std::string>>("swap_into_basis"),
153 getParam<std::string>("charge_balance_species"),
154 getParam<std::vector<std::string>>("constraint_species"),
155 getParam<std::vector<Real>>("constraint_value"),
156 getParam<MultiMooseEnum>("constraint_unit"),
157 getParam<MultiMooseEnum>("constraint_meaning"),
159 getParam<unsigned>("extra_iterations_to_make_consistent"),
160 getParam<Real>("min_initial_molality"),
161 getParam<std::vector<std::string>>("kinetic_species_name"),
162 getParam<std::vector<Real>>("kinetic_species_initial_value"),
163 getParam<MultiMooseEnum>("kinetic_species_unit")),
165 _mgd.kin_species_name.size(),
166 _is,
167 getParam<Real>("abs_tol"),
168 getParam<Real>("rel_tol"),
169 getParam<unsigned>("max_iter"),
170 getParam<Real>("max_initial_residual"),
173 getParam<std::vector<std::string>>("prevent_precipitation"),
174 getParam<Real>("max_ionic_strength"),
175 getParam<unsigned>("ramp_max_ionic_strength_initial"),
176 getParam<bool>("evaluate_kinetic_rates_always")),
178 _close_system_at_time(getParam<Real>("close_system_at_time")),
179 _closed_system(false),
180 _source_species_names(getParam<std::vector<std::string>>("source_species_names")),
182 _source_species_rates(coupledValues("source_species_rates")),
183 _remove_fixed_activity_name(getParam<std::vector<std::string>>("remove_fixed_activity_name")),
184 _remove_fixed_activity_time(getParam<std::vector<Real>>("remove_fixed_activity_time")),
188 getParam<std::vector<std::string>>("controlled_activity_name")),
190 _controlled_activity_species_values(coupledValues("controlled_activity_value")),
193 _mode(coupledValue("mode")),
195 _ramp_subsequent(getParam<unsigned>("ramp_max_ionic_strength_subsequent"))
196{
197 // check sources and set the rates
198 if (coupledComponents("source_species_rates") != _num_source_species)
199 paramError("source_species_names", "must have the same size as source_species_rates");
200 for (unsigned i = 0; i < _num_source_species; ++i)
201 if (!(_mgd.basis_species_index.count(_source_species_names[i]) == 1 ||
204 paramError("source_species_names",
205 "The name " + _source_species_names[i] +
206 " does not appear in the basis, equilibrium or kinetic species list");
207
208 // check and set activity controllers
209 for (unsigned i = 0; i < _num_removed_fixed; ++i)
210 {
213 "remove_fixed_activity_name",
214 "The species ",
216 " is not in the basis, so cannot have its activity or fugacity constraint removed");
217 else
218 {
219 const unsigned basis_ind = _mgd.basis_species_index.at(_remove_fixed_activity_name[i]);
223 paramError("remove_fixed_activity_name",
224 "The species ",
226 " is does not have an activity or fugacity constraint so cannot have such a "
227 "constraint removed");
228 }
229 }
231 paramError("remove_fixed_activity_name",
232 "must be of the same size as remove_fixed_activity_time");
233 if (coupledComponents("controlled_activity_value") != _num_controlled_activity)
234 paramError("controlled_activity_name", "must have the same size as controlled_activity_value");
235
236 // record coupled variables
237 const std::vector<MooseVariableFEBase *> & coupled_vars = getCoupledMooseVars();
238 for (unsigned int i = 0; i < coupled_vars.size(); i++)
239 addMooseVariableDependency(coupled_vars[i]);
240
241 // setup minerals_dumped
242 for (unsigned i = 0; i < _mgd.basis_species_name.size(); ++i)
245 for (unsigned j = 0; j < _mgd.eqm_species_name.size(); ++j)
248 for (unsigned k = 0; k < _mgd.kin_species_name.size(); ++k)
251}
virtual const VariableValue & coupledValue(const std::string &var_name, unsigned int comp=0) const
const std::vector< MooseVariableFieldBase * > & getCoupledMooseVars() const
std::vector< const VariableValue * > coupledValues(const std::string &var_name) const
unsigned int coupledComponents(const std::string &var_name) const
const std::vector< GeochemicalSystem::ConstraintMeaningEnum > & getConstraintMeaning() const
unsigned getNumKinetic() const
returns the number of kinetic species
ConstraintMeaningEnum
Each basis species has one of the following constraints.
Base class that controls the spatio-temporal solution of geochemistry reactions.
GeochemistryIonicStrength _is
The ionic strength calculator.
const unsigned _num_basis
number of basis species
const unsigned _max_swaps_allowed
Maximum number of swaps allowed during a single solve.
GeochemistryActivityCoefficientsDebyeHuckel _gac
The activity calculator.
ModelGeochemicalDatabase _mgd
my copy of the underlying ModelGeochemicalDatabase
GeochemistrySpeciesSwapper _swapper
The species swapper.
const Real _small_molality
A small value of molality.
const unsigned _ramp_subsequent
the ramp_max_ionic_strength to use during time-stepping
const VariableValue & _cold_temperature
Cold temperature specified by user, which is used only when mode==4.
const unsigned _heating_increments
If mode=4 then temperature is ramped from _cold_temperature to _temperature in heating_increments inc...
Real _previous_temperature
Temperature at which the _egs was last made consistent.
std::unordered_map< std::string, Real > _minerals_dumped
Moles of mineral removed by dump and flow-through.
const unsigned _num_controlled_activity
Number of species with controlled activity or fugacity.
const std::vector< const VariableValue * > _source_species_rates
Rates of the source species.
const unsigned _num_source_species
Number of source species.
DenseMatrix< Real > _dmole_additions
Derivative of moles_added.
const std::vector< Real > _remove_fixed_activity_time
Times at which to remove the fixed activity or fugacity from the species in _remove_fixed_activity_na...
const VariableValue & _temperature
Temperature specified by user.
const std::vector< std::string > _controlled_activity_species_names
Names of the species with controlled activity or fugacity.
const std::vector< std::string > _remove_fixed_activity_name
Names of species to remove the fixed activity or fugacity constraint from.
const unsigned _num_kin
Number of kinetic species.
const unsigned _num_removed_fixed
Number of elements in the vector _remove_fixed_activity_name;.
bool _closed_system
Whether the system has been closed.
DenseVector< Real > _mole_additions
Moles of each basis species added at the current timestep, along with kinetic rates.
const Real _close_system_at_time
Defines the time at which to close the system.
const std::vector< std::string > _source_species_names
Names of the source species.
std::vector< bool > _removed_fixed_activity
Whether the activity or activity constraint has been removfed.
const VariableValue & _mode
Mode of the system (flush, flow-through, etc)
const std::vector< const VariableValue * > _controlled_activity_species_values
Activity or fugacity of the species with controlled activity or fugacity.
Real _new_temperature
Temperature at which the solution is required.
GeochemicalSystem _egs
The equilibrium geochemical system that holds all the molalities, activities, etc.
const InputParameters & parameters() const
void paramError(const std::string &param, Args... args) const
const T & getParam(const std::string &name) const
void addMooseVariableDependency(MooseVariableFieldBase *var)
std::vector< bool > eqm_species_mineral
eqm_species_mineral[i] = true iff the i^th equilibrium species is a mineral
std::vector< bool > kin_species_mineral
kin_species_mineral[j] = true iff the j^th kinetic species is a mineral
std::vector< std::string > kin_species_name
kin_species_name[j] = name of the j^th kinetic species
std::vector< bool > basis_species_mineral
basis_species_mineral[j] = true iff the j^th basis species is a mineral
std::unordered_map< std::string, unsigned > eqm_species_index
eqm_species_index[name] = index of the equilibrium species (secondary aqueous species,...
std::vector< std::string > basis_species_name
basis_species_name[j] = name of the j^th basis species
std::unordered_map< std::string, unsigned > basis_species_index
basis_species_index[name] = index of the basis species, within all ModelGeochemicalDatabase internal ...
std::unordered_map< std::string, unsigned > kin_species_index
kin_species_index[name] = index of the kinetic species, within all ModelGeochemicalDatabase internal ...
std::vector< std::string > eqm_species_name
eqm_species_name[i] = name of the i^th eqm species

Member Function Documentation

◆ execute()

void GeochemistryTimeDependentReactor::execute ( )
overridevirtual

Reimplemented from GeochemistryReactorBase.

Definition at line 282 of file GeochemistryTimeDependentReactor.C.

283{
284 if (_current_node->id() != 0)
285 return;
286
288
289 _mole_additions.zero();
291
292 // remove appropriate constraints
294 {
296 _closed_system = true;
297 }
298 for (unsigned i = 0; i < _num_removed_fixed; ++i)
299 {
301 {
304 _removed_fixed_activity[i] = true;
305 }
306 }
307
308 // control activity
309 for (unsigned ca = 0; ca < _num_controlled_activity; ++ca)
310 {
311 const std::vector<GeochemicalSystem::ConstraintMeaningEnum> & cm = _egs.getConstraintMeaning();
313 {
314 const unsigned basis_ind =
319 }
320 }
321
322 // compute moles added
323 for (unsigned i = 0; i < _num_source_species; ++i)
324 {
325 const Real this_rate = (*_source_species_rates[i])[0];
327 {
328 const unsigned basis_ind = _mgd.basis_species_index.at(_source_species_names[i]);
329 _mole_additions(basis_ind) += this_rate;
330 }
332 {
333 const unsigned eqm_j = _mgd.eqm_species_index.at(_source_species_names[i]);
334 for (unsigned basis_ind = 0; basis_ind < _num_basis; ++basis_ind)
335 _mole_additions(basis_ind) += _mgd.eqm_stoichiometry(eqm_j, basis_ind) * this_rate;
336 }
337 else
338 {
339 const unsigned kin_ind = _mgd.kin_species_index.at(_source_species_names[i]);
340 _mole_additions(_num_basis + kin_ind) += this_rate;
341 }
342 }
343 for (unsigned i = 0; i < _num_basis + _num_kin; ++i)
344 _mole_additions(i) *= _dt;
345
346 // activate special modes
347 if (_mode[0] == 1.0)
348 preSolveDump();
349 else if (_mode[0] == 3.0)
351 else if (_mode[0] == 4.0)
352 {
355 }
356 else // nothing special: simply compute the desired temperature
358
359 // set temperature, if necessary
361 {
364 }
366
367 // solve the geochemical system
370 _tot_iter[0],
371 _abs_residual[0],
372 _dt,
375
376 // activate special modes
377 if (_mode[0] == 2.0)
379 else if (_mode[0] == 4.0)
381}
void solveSystem(GeochemicalSystem &egs, std::stringstream &ss, unsigned &tot_iter, Real &abs_residual, Real dt, DenseVector< Real > &mole_additions, DenseMatrix< Real > &dmole_additions)
Solve the system.
void setRampMaxIonicStrength(unsigned ramp_max_ionic_strength)
Sets the value of _ramp_max_ionic_strength.
void computeConsistentConfiguration()
Compute a reasonably consistent configuration that can be used as an initial condition in a Newton pr...
void changeConstraintToBulk(unsigned basis_ind)
Changes the constraint to MOLES_BULK_SPECIES (or MOLES_BULK_WATER if basis_ind = 0) for the basis ind...
void setConstraintValue(unsigned basis_ind, Real value)
Set the constraint value for the basis species.
void closeSystem()
Changes a KG_SOLVENT_WATER constraint to MOLES_BULK_WATER (if there is such a constraint) and all FRE...
void setTemperature(Real temperature)
Sets the temperature to the given quantity.
std::vector< unsigned > _tot_iter
Number of iterations used by the solver at each node.
std::vector< Real > _abs_residual
L1norm of the solver residual at each node.
std::vector< std::stringstream > _solver_output
The solver output at each node.
void postSolveExchanger()
This is relevant for mode=4 simulations (heat-exchanger simulations).
Real newTemperature(const DenseVector< Real > &mole_additions) const
Based on _temperature[0], mole_additions, the current mass of the system and the current temperature ...
void preSolveFlush()
Activate the special "flush" mode prior to solving the geochemical system:
void postSolveFlowThrough()
Activate the special "flow-through" mode after solving the geochemical system:
void removeCurrentSpecies()
Alter _mole_additions so that it will represent the situation in which all current species are remove...
void preSolveDump()
Activate the special "dump" mode prior to solving the geochemical system:
const Node *const & _current_node
virtual void zero() override final
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
DenseMatrix< Real > eqm_stoichiometry
eqm_stoichiometry(i, j) = stoichiometric coefficient for equilibrium species "i" in terms of the basi...

◆ finalize()

void GeochemistryTimeDependentReactor::finalize ( )
overridevirtual

Reimplemented from GeochemistryReactorBase.

Definition at line 259 of file GeochemistryTimeDependentReactor.C.

260{
262}
virtual void finalize() override

◆ getGeochemicalSystem()

const GeochemicalSystem & GeochemistryTimeDependentReactor::getGeochemicalSystem ( dof_id_type  node_id) const
overridevirtual
Returns
a reference to the equilibrium geochemical system at the given node
Parameters
node_idthe ID of the node

Implements GeochemistryReactorBase.

Definition at line 384 of file GeochemistryTimeDependentReactor.C.

385{
386 return _egs;
387}

◆ getMoleAdditions()

const DenseVector< Real > & GeochemistryTimeDependentReactor::getMoleAdditions ( dof_id_type  node_id) const
overridevirtual
Returns
the mole additions (the first num_basis of these are additions to the bulk composition of the equilibrium system, the last num_kin of these are -dt*kinetic_reaction_rate, ie dt*dissolution_rate) at the specified node_id
Parameters
node_idthe node ID

Implements GeochemistryReactorBase.

Definition at line 390 of file GeochemistryTimeDependentReactor.C.

391{
392 return _mole_additions;
393}

◆ getMolesDumped()

Real GeochemistryTimeDependentReactor::getMolesDumped ( dof_id_type  node_id,
const std::string &  species 
) const
overridevirtual
Returns
the moles dumped of the given species at the specified node_id
Parameters
speciesthe name of the species
node_idthe node ID

Implements GeochemistryReactorBase.

Definition at line 414 of file GeochemistryTimeDependentReactor.C.

416{
417 if (_minerals_dumped.count(species) == 1)
418 return _minerals_dumped.at(species);
419 return 0.0;
420}

◆ getPertinentGeochemicalSystem()

const PertinentGeochemicalSystem & GeochemistryReactorBase::getPertinentGeochemicalSystem ( ) const
inlineinherited

returns a reference to the PertinentGeochemicalSystem used to creat the ModelGeochemicalDatabase

Definition at line 73 of file GeochemistryReactorBase.h.

73{ return _pgs; };
const PertinentGeochemicalSystem & _pgs
Reference to the original PertinentGeochemicalSystem used to create the ModelGeochemicalDatabase.

Referenced by GeochemistryQuantityAux::computeValue(), GeochemistryQuantityAux::GeochemistryQuantityAux(), and GeochemistryConsoleOutput::output().

◆ getSolverIterations()

unsigned GeochemistryTimeDependentReactor::getSolverIterations ( dof_id_type  node_id) const
overridevirtual
Returns
the total number of iterations used by the most recent solve at the specified node_id
Parameters
node_idthe node ID

Implements GeochemistryReactorBase.

Definition at line 402 of file GeochemistryTimeDependentReactor.C.

403{
404 return _tot_iter[0];
405}

◆ getSolverOutput()

const std::stringstream & GeochemistryTimeDependentReactor::getSolverOutput ( dof_id_type  node_id) const
overridevirtual
Returns
a reference to the most recent solver output (containing iteration info, swap info, residuals, etc) at the specified node_id
Parameters
node_idthe node ID

Implements GeochemistryReactorBase.

Definition at line 396 of file GeochemistryTimeDependentReactor.C.

397{
398 return _solver_output[0];
399}

◆ getSolverResidual()

Real GeochemistryTimeDependentReactor::getSolverResidual ( dof_id_type  node_id) const
overridevirtual
Returns
the L1norm of the residual at the end of the most recent solve at the specified node_id
Parameters
node_idthe node ID

Implements GeochemistryReactorBase.

Definition at line 408 of file GeochemistryTimeDependentReactor.C.

409{
410 return _abs_residual[0];
411}

◆ initialize()

void GeochemistryTimeDependentReactor::initialize ( )
overridevirtual

Reimplemented from GeochemistryReactorBase.

Definition at line 254 of file GeochemistryTimeDependentReactor.C.

◆ initialSetup()

void GeochemistryTimeDependentReactor::initialSetup ( )
overridevirtual

Reimplemented from NodalUserObject.

Definition at line 265 of file GeochemistryTimeDependentReactor.C.

266{
267 // solve the geochemical system with its initial composition and with dt=0 so no kinetic additions
268 if (_num_my_nodes == 0)
269 return; // rather peculiar case where user has used many processors
270 _mole_additions.zero();
274 _tot_iter[0],
275 _abs_residual[0],
276 0.0,
279}
const unsigned _num_my_nodes
Number of nodes handled by this processor (will need to be made un-const when mesh adaptivity is hand...

◆ newTemperature()

Real GeochemistryTimeDependentReactor::newTemperature ( const DenseVector< Real > &  mole_additions) const
protected

Based on _temperature[0], mole_additions, the current mass of the system and the current temperature of the system, work out the resulting temperature.

This is either:

  • _temperature[0] if no mole_additions are positive (nothing is being added to the system)
  • a weighted sum of _temperature[0] (the temperature of the reactants being added) and the current system temperature, assuming the heat capacities of the inputs equal the heat capacities of the existing constituents in the system
    Parameters
    mole_additionsNumber of moles of each species being added to the system

Definition at line 423 of file GeochemistryTimeDependentReactor.C.

424{
426 return _temperature[0];
427
428 // if no reactants are being added, the system temperature will be _temperature[0]
429 bool any_additions = false;
430 for (unsigned i = 0; i < _num_basis + _num_kin; ++i)
431 if (mole_additions(i) > 0)
432 {
433 any_additions = true;
434 break;
435 }
436 if (!any_additions)
437 return _temperature[0];
438
439 // assume heat capacities of inputs and outputs are the same, so final temperature is dictated
440 // by masses, also assume that the input happens first, then temperature equilibration, then
441 // the outputs occur
442 Real new_temperature = _temperature[0];
443 const std::vector<Real> & current_bulk = _egs.getBulkMolesOld();
444 Real current_kg = current_bulk[0] / GeochemistryConstants::MOLES_PER_KG_WATER;
445 Real input_kg = std::max(mole_additions(0), 0.0) / GeochemistryConstants::MOLES_PER_KG_WATER;
446 for (unsigned i = 1; i < _num_basis; ++i)
447 {
448 current_kg += current_bulk[i] * _mgd.basis_species_molecular_weight[i] / 1000.0;
449 input_kg += std::max(mole_additions(i), 0.0) * _mgd.basis_species_molecular_weight[i] / 1000.0;
450 }
451 for (unsigned k = 0; k < _num_kin; ++k)
452 input_kg += std::max(mole_additions(k + _num_basis), 0.0) *
454 new_temperature =
455 (_previous_temperature * current_kg + _temperature[0] * input_kg) / (current_kg + input_kg);
456 return new_temperature;
457}
const std::vector< Real > & getBulkMolesOld() const
std::vector< Real > basis_species_molecular_weight
all quantities have a molecular weight (g)
std::vector< Real > kin_species_molecular_weight
all quantities have a molecular weight (g/mol)

Referenced by execute(), preSolveDump(), and preSolveFlush().

◆ postSolveExchanger()

void GeochemistryTimeDependentReactor::postSolveExchanger ( )
protected

This is relevant for mode=4 simulations (heat-exchanger simulations).

Any precipitates are removed from the system, the system is then heated to temperature and re-solved, and any precipitates are removed from the system

Definition at line 604 of file GeochemistryTimeDependentReactor.C.

605{
606 // remove precipitates
608
609 DenseVector<Real> zero_mole_additions(_num_basis + _num_kin);
610 DenseMatrix<Real> zero_dmole_additions(_num_basis + _num_kin, _num_basis + _num_kin);
611
612 const Real del_temperature = (_temperature[0] - _previous_temperature) / _heating_increments;
613
614 for (unsigned incr = 0; incr < _heating_increments; ++incr)
615 {
616 // set temperature
617 _new_temperature = _previous_temperature + del_temperature;
621
622 zero_mole_additions.zero();
623 zero_dmole_additions.zero();
624
625 // solve the geochemical system
628 _tot_iter[0],
629 _abs_residual[0],
630 _dt,
631 zero_mole_additions,
632 zero_dmole_additions);
633
634 // remove precipitates
636 }
637}

Referenced by execute().

◆ postSolveFlowThrough()

void GeochemistryTimeDependentReactor::postSolveFlowThrough ( )
protected

Activate the special "flow-through" mode after solving the geochemical system:

  • precipitated minerals are removed by setting the free mole number of minerals (basis and kinetic) is set to a small number
  • _minerals_dumped is updated accordingly

Definition at line 560 of file GeochemistryTimeDependentReactor.C.

561{
562 // copy the current_molal values into a new vector
563 const std::vector<Real> current_molal = _egs.getSolventMassAndFreeMolalityAndMineralMoles();
564 // remove minerals
565 for (unsigned i = 1; i < _num_basis; ++i)
567 {
568 const Real to_remove = current_molal[i] - _small_molality * 10.0;
569 _egs.addToBulkMoles(i, -to_remove);
571 }
572
573 // copy the current kinetic moles into a new vector
574 const std::vector<Real> kin_moles = _egs.getKineticMoles();
575 // remove minerals
576 for (unsigned k = 0; k < _num_kin; ++k)
578 {
579 const Real to_remove = kin_moles[k] - _small_molality;
580 for (unsigned i = 0; i < _num_basis; ++i)
581 if (_mgd.kin_stoichiometry(k, i) != 0)
583 -_mgd.kin_stoichiometry(k, i) *
584 to_remove); // remember bulk moles contains kinetic contributions
586 _minerals_dumped[_mgd.kin_species_name[k]] += to_remove;
587 }
588
590}
void setMineralRelatedFreeMoles(Real value)
Set the free mole number of mineral-related species to the value provided.
void setKineticMoles(unsigned kin, Real moles)
Sets the current AND old mole number for a kinetic species.
Real getKineticMoles(unsigned kin) const
void addToBulkMoles(unsigned basis_ind, Real value)
Add to the MOLES_BULK_SPECIES (or MOLES_BULK_WATER if basis_ind = 0) for the basis species.
const std::vector< Real > & getSolventMassAndFreeMolalityAndMineralMoles() const
DenseMatrix< Real > kin_stoichiometry
kin_stoichiometry(i, j) = stoichiometric coefficient for kinetic species "i" in terms of the basis sp...

Referenced by execute(), and postSolveExchanger().

◆ preSolveDump()

void GeochemistryTimeDependentReactor::preSolveDump ( )
protected

Activate the special "dump" mode prior to solving the geochemical system:

  • _mole_additions for basis mineral moles are set to the free mineral mole numbers (potentially overwriting a user's choice)
  • _minerals_dumped is updated
  • the new temperature is computed
  • the free mineral moles and bulk mole numbers in the GeochemicalSystem are updated accordingly, and _mole_additions is zero
  • all minerals are swapped out of the basis

Definition at line 460 of file GeochemistryTimeDependentReactor.C.

461{
462 // remove basis mineral moles
463 const std::vector<Real> & current_molal = _egs.getSolventMassAndFreeMolalityAndMineralMoles();
464 for (unsigned i = 1; i < _num_basis; ++i)
466 {
467 _mole_additions(i) = -current_molal[i]; // might overwrite the rates set above, which is good
468 _minerals_dumped[_mgd.basis_species_name[i]] += current_molal[i];
469 }
470
472
473 // add the chemicals immediately instead of during the solve (as occurs for other modes)
474 for (unsigned basis_ind = 0; basis_ind < _num_basis; ++basis_ind)
475 {
476 _egs.addToBulkMoles(basis_ind, _mole_additions(basis_ind));
477 _mole_additions(basis_ind) = 0.0;
478 }
479 // dump needs free mineral moles to be exactly zero and the above addToBulkMoles will have set
480 // this for standard minerals, but not things related to sorption or kinetic minerals, so:
482
483 // Now need to swap all minerals out of the basis
484 const std::vector<Real> & eqm_molality = _egs.getEquilibriumMolality();
485 unsigned swap_into_basis = 0;
486 for (unsigned i = 0; i < _num_basis; ++i)
488 {
489 const bool legitimate_swap_found = _egs.getSwapper().findBestEqmSwap(
490 i, _mgd, eqm_molality, false, false, false, swap_into_basis);
491 if (legitimate_swap_found)
492 {
493 try
494 {
495 _egs.performSwap(i, swap_into_basis);
496 }
497 catch (const MooseException & e)
498 {
499 mooseError(e.what());
500 }
501 }
502 }
503}
const GeochemistrySpeciesSwapper & getSwapper() const
void performSwap(unsigned swap_out_of_basis, unsigned swap_into_basis)
Perform the basis swap, and ensure that the resulting system is consistent.
Real getEquilibriumMolality(unsigned j) const
bool findBestEqmSwap(unsigned basis_ind, const ModelGeochemicalDatabase &mgd, const std::vector< Real > &eqm_molality, bool minerals_allowed, bool gas_allowed, bool sorption_allowed, unsigned &best_eqm_species) const
For the the given basis index, find the equilibrium index that it should be swapped with.
void mooseError(Args &&... args) const
virtual const char * what() const

Referenced by execute().

◆ preSolveFlush()

void GeochemistryTimeDependentReactor::preSolveFlush ( )
protected

Activate the special "flush" mode prior to solving the geochemical system:

  • the new temperature is computed
  • the total mass of aqueous solution entering the system (ie, without free mineral moles and kinetic-mineral moles) is computed
  • the current mass of the aqueous solution is computed
  • _mole_additions is updated by removing species from the current aqueous solution so that the final mass (after _mole_additions has been added) will be the same as the current mass of aqueous solution

Definition at line 506 of file GeochemistryTimeDependentReactor.C.

507{
509
510 // Here we conserve mass, so compute the mass of the solution, without the free mineral moles.
511 // We don't include the free mineral moles because users of GeochemistWorkbench will want
512 // "flush" to operate like Bethke Eqn(13.14)
513 // I assume we also don't include kinetic-mineral moles
515 for (unsigned i = 1; i < _num_basis; ++i)
517 kg_in += _mole_additions(i) * _mgd.basis_species_molecular_weight[i] / 1000.0;
518 for (unsigned kin = 0; kin < _num_kin; ++kin)
519 if (!_mgd.kin_species_mineral[kin])
520 kg_in += _mole_additions(kin + _num_basis) * _mgd.kin_species_molecular_weight[kin] / 1000.0;
521
522 const std::vector<Real> & current_bulk = _egs.getBulkMolesOld();
523 const std::vector<Real> & current_molal = _egs.getSolventMassAndFreeMolalityAndMineralMoles();
524 const std::vector<Real> & kin_moles = _egs.getKineticMoles();
525
526 // compute the current mass, without moles from free minerals and without kinetic minerals
527 Real current_kg = current_bulk[0] / GeochemistryConstants::MOLES_PER_KG_WATER;
528 for (unsigned i = 1; i < _num_basis; ++i)
529 {
530 Real kinetic_contribution = 0.0;
531 for (unsigned k = 0; k < _num_kin; ++k)
533 kinetic_contribution += kin_moles[k] * _mgd.kin_stoichiometry(k, i);
535 current_kg += (current_bulk[i] - current_molal[i] - kinetic_contribution) *
537 else
538 current_kg += (current_bulk[i] - kinetic_contribution) *
540 }
541
542 const Real fraction_to_remove = kg_in / current_kg;
543 for (unsigned i = 0; i < _num_basis; ++i)
544 {
545 Real all_kinetic_contribution = 0.0;
546 for (unsigned k = 0; k < _num_kin; ++k)
547 all_kinetic_contribution += kin_moles[k] * _mgd.kin_stoichiometry(k, i);
549 _mole_additions(i) -=
550 fraction_to_remove * (current_bulk[i] - current_molal[i] - all_kinetic_contribution);
551 else
552 _mole_additions(i) -= fraction_to_remove * (current_bulk[i] - all_kinetic_contribution);
553 }
554 for (unsigned k = 0; k < _num_kin; ++k)
556 _mole_additions(k + _num_basis) -= fraction_to_remove * kin_moles[k];
557}

Referenced by execute().

◆ removeCurrentSpecies()

void GeochemistryTimeDependentReactor::removeCurrentSpecies ( )
protected

Alter _mole_additions so that it will represent the situation in which all current species are removed and replaced with source_species.

Definition at line 593 of file GeochemistryTimeDependentReactor.C.

594{
595 const std::vector<Real> & current_bulk = _egs.getBulkMolesOld();
596 for (unsigned i = 0; i < _num_basis; ++i)
597 _mole_additions(i) -= current_bulk[i];
598 const std::vector<Real> & kin_moles = _egs.getKineticMoles();
599 for (unsigned k = 0; k < _num_kin; ++k)
600 _mole_additions(k + _num_basis) -= kin_moles[k];
601}

Referenced by execute().

◆ sharedParams()

InputParameters GeochemistryTimeDependentReactor::sharedParams ( )
static

params that are shared with AddTimeDependentReactionSolverAction

Definition at line 15 of file GeochemistryTimeDependentReactor.C.

16{
18 params.addParam<unsigned>(
19 "ramp_max_ionic_strength_subsequent",
20 0,
21 "The number of iterations over which to progressively increase the maximum ionic strength "
22 "(from zero to max_ionic_strength) during time-stepping. Unless a great deal occurs in each "
23 "time step, this parameter can be set quite small");
24 params.addCoupledVar(
25 "mode",
26 0.0,
27 "This may vary temporally. If mode=1 then 'dump' mode is used, which means all non-kinetic "
28 "mineral masses are removed from the system before the equilibrium solution is sought (ie, "
29 "removal occurs at the beginning of the time step). If mode=2 then 'flow-through' mode is "
30 "used, which means all mineral masses are removed from the system after it the "
31 "equilbrium solution has been found (ie, at the end of a time step). If mode=3 then 'flush' "
32 "mode is used, then before the equilibrium solution is sought (ie, at the start of a time "
33 "step) water+species is removed from the system at the same rate as pure water + non-mineral "
34 "solutes are entering the system (specified in source_species_rates). If mode=4 then "
35 "'heat-exchanger' mode is used, which means the entire current aqueous solution is removed, "
36 "then the source_species are added, then the temperature is set to 'cold_temperature', the "
37 "system is solved and any precipitated minerals are removed, then the temperature is set to "
38 "'temperature', the system re-solved and any precipitated minerals are removed. If mode is "
39 "any other number, no special mode is active (the system simply responds to the "
40 "source_species_rates, controlled_activity_value, etc).");
41 params.addCoupledVar(
42 "temperature",
43 25,
44 "Temperature. This has two different meanings if mode!=4. (1) If no species are being "
45 "added to the solution (no source_species_rates are positive) then this is the temperature "
46 "of the aqueous solution. (2) If species are being added, this is the temperature of the "
47 "species being added. In case (2), the final aqueous-solution temperature is computed "
48 "assuming the species are added, temperature is equilibrated and then, if species are also "
49 "being removed, they are removed. If you wish to add species and simultaneously alter the "
50 "temperature, you will have to use a sequence of heat-add-heat-add, etc steps. In the case "
51 "that mode=4, temperature is the final temperature of the aqueous solution");
52 params.addCoupledVar(
53 "cold_temperature",
54 25,
55 "This is only used if mode=4, where it is the cold temperature of the heat exchanger.");
56 params.addRangeCheckedParam<unsigned>(
57 "heating_increments",
58 1,
59 "heating_increments > 0",
60 "This is only used if mode=4. Internal to this object, the temperature is ramped from "
61 "cold_temperature to temperature in heating_increments increments. This helps difficult "
62 "problems converge");
63 params.addParam<Real>("initial_temperature",
64 25.0,
65 "The initial aqueous solution is equilibrated at this system before adding "
66 "reactants, changing temperature, etc.");
67 params.addParam<Real>("close_system_at_time",
68 0.0,
69 "Time at which to 'close' the system, that is, change a kg_solvent_water "
70 "constraint to moles_bulk_water, and all free_molality and "
71 "free_moles_mineral_species to moles_bulk_species");
72 params.addParam<std::vector<std::string>>(
73 "remove_fixed_activity_name",
74 {},
75 "The name of the species that should have their activity or fugacity constraint removed at "
76 "time given in remove_fixed_activity_time. There should be an equal number of these names "
77 "as times given in remove_fixed_activity_time. Each of these must be in the basis and have "
78 "an activity or fugacity constraint");
79 params.addParam<std::vector<Real>>("remove_fixed_activity_time",
80 {},
81 "The times at which the species in remove_fixed_activity_name "
82 "should have their activity or fugacity constraint removed.");
83 params.addParam<std::vector<std::string>>(
84 "source_species_names",
85 {},
86 "The name of the species that are added at rates given in source_species_rates. There must "
87 "be an equal number of these as source_species_rates.");
88 params.addCoupledVar("source_species_rates",
89 "Rates, in mols/time_unit, of addition of the species with names given in "
90 "source_species_names. A negative value corresponds to removing a species: "
91 "be careful that you don't cause negative mass problems!");
92 params.addParam<std::vector<std::string>>(
93 "controlled_activity_name",
94 {},
95 "The names of the species that have their activity or fugacity constrained. There should be "
96 "an equal number of these names as values given in controlled_activity_value. NOTE: if "
97 "these species are not in the basis, or they do not have an activity (or fugacity) "
98 "constraint then their activity cannot be controlled: in this case MOOSE will ignore the "
99 "value you prescribe in controlled_activity_value.");
100 params.addCoupledVar("controlled_activity_value",
101 "Values of the activity or fugacity of the species in "
102 "controlled_activity_name list. These should always be positive");
103 params.addParam<bool>(
104 "evaluate_kinetic_rates_always",
105 true,
106 "If true, then, evaluate the kinetic rates at every Newton step during the solve using the "
107 "current values of molality, activity, etc (ie, implement an implicit solve). If false, "
108 "then evaluate the kinetic rates using the values of molality, activity, etc, at the start "
109 "of the current time step (ie, implement an explicit solve)");
110 params.addParam<std::vector<std::string>>(
111 "kinetic_species_name",
112 {},
113 "Names of the kinetic species given initial values in kinetic_species_initial_value");
114 params.addParam<std::vector<Real>>(
115 "kinetic_species_initial_value",
116 {},
117 "Initial number of moles, mass or volume (depending on kinetic_species_unit) for each of the "
118 "species named in kinetic_species_name");
119 MultiMooseEnum kin_species_unit("dimensionless moles molal kg g mg ug kg_per_kg_solvent "
120 "g_per_kg_solvent mg_per_kg_solvent ug_per_kg_solvent cm3");
121 params.addParam<MultiMooseEnum>(
122 "kinetic_species_unit",
123 kin_species_unit,
124 "Units of the numerical values given in kinetic_species_initial_value. Moles: mole number. "
125 "kg: kilograms. g: grams. mg: milligrams. ug: micrograms. cm3: cubic centimeters");
126 return params;
127}
InputParameters emptyInputParameters()
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void addCoupledVar(const std::string &name, const std::string &doc_string)
void addRangeCheckedParam(const std::string &name, const T &value, const std::string &parsed_function, const std::string &doc_string)

Referenced by AddTimeDependentReactionSolverAction::validParams(), and validParams().

◆ threadJoin()

virtual void GeochemistryReactorBase::threadJoin ( const UserObject )
inlineoverridevirtualinherited

Implements NodalUserObject.

Reimplemented in GeochemistrySpatialReactor.

Definition at line 28 of file GeochemistryReactorBase.h.

28{}

◆ validParams()

InputParameters GeochemistryTimeDependentReactor::validParams ( )
static

Definition at line 130 of file GeochemistryTimeDependentReactor.C.

131{
134 params.addClassDescription("UserObject that controls the time-dependent geochemistry reaction "
135 "processes. Spatial dependence is not possible using this class");
136 return params;
137}
static InputParameters validParams()
static InputParameters sharedParams()
params that are shared with AddTimeDependentReactionSolverAction
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _abs_residual

std::vector<Real> GeochemistryReactorBase::_abs_residual
protectedinherited

◆ _close_system_at_time

const Real GeochemistryTimeDependentReactor::_close_system_at_time
protected

Defines the time at which to close the system.

Definition at line 58 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _closed_system

bool GeochemistryTimeDependentReactor::_closed_system
protected

Whether the system has been closed.

Definition at line 60 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _cold_temperature

const VariableValue& GeochemistryTimeDependentReactor::_cold_temperature
protected

Cold temperature specified by user, which is used only when mode==4.

Definition at line 44 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _controlled_activity_species_names

const std::vector<std::string> GeochemistryTimeDependentReactor::_controlled_activity_species_names
protected

Names of the species with controlled activity or fugacity.

Definition at line 76 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _controlled_activity_species_values

const std::vector<const VariableValue *> GeochemistryTimeDependentReactor::_controlled_activity_species_values
protected

Activity or fugacity of the species with controlled activity or fugacity.

Definition at line 80 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _dmole_additions

DenseMatrix<Real> GeochemistryTimeDependentReactor::_dmole_additions
protected

Derivative of moles_added.

Definition at line 84 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and initialSetup().

◆ _egs

GeochemicalSystem GeochemistryTimeDependentReactor::_egs
protected

The equilibrium geochemical system that holds all the molalities, activities, etc.

Definition at line 52 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), GeochemistryTimeDependentReactor(), getGeochemicalSystem(), initialSetup(), newTemperature(), postSolveExchanger(), postSolveFlowThrough(), preSolveDump(), preSolveFlush(), and removeCurrentSpecies().

◆ _gac

GeochemistryActivityCoefficientsDebyeHuckel GeochemistryReactorBase::_gac
protectedinherited

The activity calculator.

Definition at line 91 of file GeochemistryReactorBase.h.

Referenced by GeochemistrySpatialReactor::GeochemistrySpatialReactor().

◆ _heating_increments

const unsigned GeochemistryTimeDependentReactor::_heating_increments
protected

If mode=4 then temperature is ramped from _cold_temperature to _temperature in heating_increments increments.

Definition at line 46 of file GeochemistryTimeDependentReactor.h.

Referenced by postSolveExchanger().

◆ _initial_max_ionic_str

const Real GeochemistryReactorBase::_initial_max_ionic_str
protectedinherited

Initial value of maximum ionic strength.

Definition at line 87 of file GeochemistryReactorBase.h.

◆ _is

GeochemistryIonicStrength GeochemistryReactorBase::_is
protectedinherited

The ionic strength calculator.

Definition at line 89 of file GeochemistryReactorBase.h.

Referenced by GeochemistrySpatialReactor::GeochemistrySpatialReactor().

◆ _max_swaps_allowed

const unsigned GeochemistryReactorBase::_max_swaps_allowed
protectedinherited

Maximum number of swaps allowed during a single solve.

Definition at line 93 of file GeochemistryReactorBase.h.

◆ _mgd

ModelGeochemicalDatabase GeochemistryReactorBase::_mgd
protectedinherited

◆ _minerals_dumped

std::unordered_map<std::string, Real> GeochemistryTimeDependentReactor::_minerals_dumped
protected

Moles of mineral removed by dump and flow-through.

Definition at line 88 of file GeochemistryTimeDependentReactor.h.

Referenced by GeochemistryTimeDependentReactor(), getMolesDumped(), postSolveFlowThrough(), and preSolveDump().

◆ _mode

const VariableValue& GeochemistryTimeDependentReactor::_mode
protected

Mode of the system (flush, flow-through, etc)

Definition at line 86 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _mole_additions

DenseVector<Real> GeochemistryTimeDependentReactor::_mole_additions
protected

Moles of each basis species added at the current timestep, along with kinetic rates.

Definition at line 82 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), getMoleAdditions(), initialSetup(), preSolveDump(), preSolveFlush(), and removeCurrentSpecies().

◆ _new_temperature

Real GeochemistryTimeDependentReactor::_new_temperature
protected

Temperature at which the solution is required.

Definition at line 48 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), postSolveExchanger(), preSolveDump(), and preSolveFlush().

◆ _num_basis

const unsigned GeochemistryReactorBase::_num_basis
protectedinherited

◆ _num_controlled_activity

const unsigned GeochemistryTimeDependentReactor::_num_controlled_activity
protected

Number of species with controlled activity or fugacity.

Definition at line 78 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and GeochemistryTimeDependentReactor().

◆ _num_eqm

const unsigned GeochemistryReactorBase::_num_eqm
protectedinherited

number of equilibrium species

Definition at line 85 of file GeochemistryReactorBase.h.

◆ _num_kin

const unsigned GeochemistryTimeDependentReactor::_num_kin
protected

◆ _num_my_nodes

const unsigned GeochemistryReactorBase::_num_my_nodes
protectedinherited

◆ _num_removed_fixed

const unsigned GeochemistryTimeDependentReactor::_num_removed_fixed
protected

Number of elements in the vector _remove_fixed_activity_name;.

Definition at line 72 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and GeochemistryTimeDependentReactor().

◆ _num_source_species

const unsigned GeochemistryTimeDependentReactor::_num_source_species
protected

Number of source species.

Definition at line 64 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and GeochemistryTimeDependentReactor().

◆ _pgs

const PertinentGeochemicalSystem& GeochemistryReactorBase::_pgs
protectedinherited

Reference to the original PertinentGeochemicalSystem used to create the ModelGeochemicalDatabase.

Definition at line 81 of file GeochemistryReactorBase.h.

Referenced by GeochemistryReactorBase::getPertinentGeochemicalSystem().

◆ _previous_temperature

Real GeochemistryTimeDependentReactor::_previous_temperature
protected

Temperature at which the _egs was last made consistent.

Definition at line 50 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), newTemperature(), and postSolveExchanger().

◆ _ramp_subsequent

const unsigned GeochemistryTimeDependentReactor::_ramp_subsequent
protected

the ramp_max_ionic_strength to use during time-stepping

Definition at line 90 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _remove_fixed_activity_name

const std::vector<std::string> GeochemistryTimeDependentReactor::_remove_fixed_activity_name
protected

Names of species to remove the fixed activity or fugacity constraint from.

Definition at line 68 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and GeochemistryTimeDependentReactor().

◆ _remove_fixed_activity_time

const std::vector<Real> GeochemistryTimeDependentReactor::_remove_fixed_activity_time
protected

Times at which to remove the fixed activity or fugacity from the species in _remove_fixed_activity_name.

Definition at line 70 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and GeochemistryTimeDependentReactor().

◆ _removed_fixed_activity

std::vector<bool> GeochemistryTimeDependentReactor::_removed_fixed_activity
protected

Whether the activity or activity constraint has been removfed.

Definition at line 74 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _small_molality

const Real GeochemistryReactorBase::_small_molality
protectedinherited

A small value of molality.

Definition at line 97 of file GeochemistryReactorBase.h.

Referenced by postSolveFlowThrough().

◆ _solver

GeochemicalSolver GeochemistryTimeDependentReactor::_solver
protected

The solver.

Definition at line 54 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), initialSetup(), and postSolveExchanger().

◆ _solver_output

std::vector<std::stringstream> GeochemistryReactorBase::_solver_output
protectedinherited

◆ _source_species_names

const std::vector<std::string> GeochemistryTimeDependentReactor::_source_species_names
protected

Names of the source species.

Definition at line 62 of file GeochemistryTimeDependentReactor.h.

Referenced by execute(), and GeochemistryTimeDependentReactor().

◆ _source_species_rates

const std::vector<const VariableValue *> GeochemistryTimeDependentReactor::_source_species_rates
protected

Rates of the source species.

Definition at line 66 of file GeochemistryTimeDependentReactor.h.

Referenced by execute().

◆ _swapper

GeochemistrySpeciesSwapper GeochemistryReactorBase::_swapper
protectedinherited

The species swapper.

Definition at line 95 of file GeochemistryReactorBase.h.

Referenced by GeochemistrySpatialReactor::GeochemistrySpatialReactor().

◆ _temperature

const VariableValue& GeochemistryTimeDependentReactor::_temperature
protected

Temperature specified by user.

Definition at line 42 of file GeochemistryTimeDependentReactor.h.

Referenced by newTemperature(), and postSolveExchanger().

◆ _tot_iter

std::vector<unsigned> GeochemistryReactorBase::_tot_iter
protectedinherited

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