www.mooseframework.org
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Types | Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
PorousFlowHeatAdvection Class Reference

Advection of heat via flux of component k in fluid phase alpha. More...

#include <PorousFlowHeatAdvection.h>

Inheritance diagram for PorousFlowHeatAdvection:
[legend]

Public Types

enum  ResidualTagType { ResidualTagType::NonReference, ResidualTagType::Reference }
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 PorousFlowHeatAdvection (const InputParameters &parameters)
 
virtual void computeResidualAndJacobian () override
 
virtual void computeOffDiagJacobianScalar (unsigned int jvar) override
 
virtual const MooseVariablevariable () const override
 
virtual void computeNonlocalJacobian ()
 
virtual void computeNonlocalOffDiagJacobian (unsigned int)
 
const SubProblemsubProblem () const
 
virtual void prepareShapes (unsigned int var_num)
 
virtual bool enabled () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
virtual const std::string & name () const
 
std::string typeAndName () const
 
std::string errorPrefix (const std::string &error_type) const
 
void callMooseError (std::string msg, const bool with_prefix) const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
const InputParametersparameters () const
 
MooseObjectName uniqueName () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &nm) const
 
void paramError (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
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void mooseError (Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &name, const std::string *param=nullptr) const
 
virtual void initialSetup ()
 
virtual void jacobianSetup ()
 
virtual void residualSetup ()
 
virtual void subdomainSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnumgetExecuteOnEnum () 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
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const T & getUserObject (const std::string &param_name, bool is_dependency=true) const
 
const T & getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
bool isImplicit ()
 
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
 
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)
 
virtual void meshChanged ()
 
void useVectorTag (const TagName &tag_name, VectorTagsKey)
 
void useVectorTag (TagID tag_id, VectorTagsKey)
 
void useMatrixTag (const TagName &tag_name, MatrixTagsKey)
 
void useMatrixTag (TagID tag_id, MatrixTagsKey)
 
bool isVectorTagged ()
 
bool isMatrixTagged ()
 
bool hasVectorTags () const
 
const std::set< TagID > & getVectorTags (VectorTagsKey) const
 
const std::set< TagID > & getMatrixTags (MatrixTagsKey) 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 (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
 
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 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={})
 
std::set< MooseVariableFieldBase *> checkVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const MaterialProperty< T > & getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const ADMaterialProperty< T > & getADMaterialPropertyByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOldByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
const MaterialProperty< T > & getMaterialPropertyOlderByName (const MaterialPropertyName &name)
 
std::pair< const MaterialProperty< T > *, std::set< SubdomainID > > getBlockMaterialProperty (const MaterialPropertyName &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialProperty ()
 
const GenericMaterialProperty< T, is_ad > & getGenericZeroMaterialPropertyByName (const std::string &prop_name)
 
const MaterialProperty< T > & getZeroMaterialProperty (Ts... args)
 
std::set< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (const std::string &name)
 
std::vector< BoundaryName > getMaterialPropertyBoundaryNames (const std::string &name)
 
void checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)
 
std::unordered_map< SubdomainID, std::vector< MaterialBase *> > buildRequiredMaterials (bool allow_stateful=true)
 
void statefulPropertiesAllowed (bool)
 
bool getMaterialPropertyCalled () const
 
const std::unordered_set< unsigned int > & getMatPropDependencies () const
 
virtual void resolveOptionalProperties ()
 
const GenericMaterialProperty< T, is_ad > & getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)
 
virtual unsigned int getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const
 
virtual unsigned int getElementIDIndexByName (const std::string &id_name) const
 
virtual const dof_id_typegetElementID (const std::string &id_parameter_name, unsigned int comp=0) const
 
dof_id_type getElementID (const Elem *elem, unsigned int elem_id_index) const
 
virtual const dof_id_typegetElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const
 
virtual const dof_id_typegetElementIDByName (const std::string &id_name) const
 
virtual const dof_id_typegetElementIDNeighborByName (const std::string &id_name) const
 
bool hasElementID (const std::string &id_name) const
 
dof_id_type maxElementID (unsigned int elem_id_index) const
 
dof_id_type minElementID (unsigned int elem_id_index) const
 
bool areElemIDsIdentical (const std::string &id_name1, const std::string &id_name2) const
 
std::unordered_map< dof_id_type, std::set< dof_id_type > > getElemIDMapping (const std::string &id_name1, const std::string &id_name2) const
 
std::set< dof_id_typegetAllElemIDs (unsigned int elem_id_index) const
 
std::set< dof_id_typegetElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) const
 
MooseVariableBasemooseVariableBase () const
 
MooseVariableField< Real > & mooseVariableField ()
 
MooseVariableFE< Real > * mooseVariable () const
 
MooseVariableFV< Real > * mooseVariableFV () const
 
MooseLinearVariableFV< Real > * mooseLinearVariableFV () 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 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 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 GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericOptionalMaterialProperty< T, is_ad > & getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalMaterialProperty< T > & getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalADMaterialProperty< T > & getOptionalADMaterialProperty (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOld (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
const OptionalMaterialProperty< T > & getOptionalMaterialPropertyOlder (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialProperty (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialProperty (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasADMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialProperty (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
bool hasGenericMaterialPropertyByName (const std::string &name)
 
PenetrationLocatorgetPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
PenetrationLocatorgetQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)
 
NearestNodeLocatorgetNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
NearestNodeLocatorgetQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)
 
bool requiresGeometricSearch () const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

static InputParameters validParams ()
 

Public Attributes

const ConsoleStream _console
 

Static Public Attributes

static constexpr PropertyValue::id_type default_property_id
 
static constexpr PropertyValue::id_type zero_property_id
 

Protected Types

enum  JacRes { JacRes::CALCULATE_RESIDUAL = 0, JacRes::CALCULATE_JACOBIAN = 1 }
 
enum  FallbackEnum { FallbackEnum::QUICK, FallbackEnum::HARMONIC }
 If full upwinding is failing due to nodes swapping between upwind and downwind in successive nonlinear iterations (within one timestep and element) a fallback scheme is used instead. More...
 

Protected Member Functions

virtual Real mobility (unsigned nodenum, unsigned phase) const override
 The mobility of the fluid. More...
 
virtual Real dmobility (unsigned nodenum, unsigned phase, unsigned pvar) const override
 The derivative of mobility with respect to PorousFlow variable pvar. More...
 
virtual void timestepSetup () override
 
virtual Real computeQpResidual () override
 
virtual void computeResidual () override
 
virtual void computeJacobian () override
 
virtual void computeOffDiagJacobian (unsigned int jvar) override
 
virtual Real darcyQp (unsigned int ph) const
 The Darcy part of the flux (this is the non-upwinded part) More...
 
virtual Real darcyQpJacobian (unsigned int jvar, unsigned int ph) const
 Jacobian of the Darcy part of the flux. More...
 
void computeResidualAndJacobian (JacRes res_or_jac, unsigned int jvar)
 Computation of the residual and Jacobian. More...
 
virtual void computeResidualAndJacobian ()
 
void fullyUpwind (JacRes res_or_jac, unsigned int ph, unsigned int pvar)
 Calculate the residual or Jacobian using full upwinding. More...
 
void quickUpwind (JacRes res_or_jac, unsigned int ph, unsigned int pvar)
 Calculate the residual or Jacobian using the nodal mobilities, but without conserving fluid mass. More...
 
void harmonicMean (JacRes res_or_jac, unsigned int ph, unsigned int pvar)
 Calculate the residual or Jacobian by using the harmonic mean of the nodal mobilities for the entire element. More...
 
virtual Real computeQpJacobian ()
 
virtual Real computeQpOffDiagJacobian (unsigned int)
 
virtual Real computeQpOffDiagJacobianScalar (unsigned int)
 
virtual RealEigenVector computeQpOffDiagJacobianArray (const ArrayMooseVariable &jvar)
 
virtual void precalculateResidual ()
 
virtual void precalculateJacobian ()
 
virtual void precalculateOffDiagJacobian (unsigned int)
 
const MooseVariableFieldBasegetVariable (unsigned int jvar_num) const
 
virtual void addUserObjectDependencyHelper (const UserObject &) const
 
Moose::StateArg determineState () const
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &) const
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const
 
T & declareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< T > declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const T & getRestartableData (const std::string &data_name) const
 
T & declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
T & declareRecoverableData (const std::string &data_name, Args &&... args)
 
T & declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
T & declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTagNeighbor (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTagLower (Assembly &assembly, unsigned int ivar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar, DenseMatrix< Number > &k) const
 
void prepareMatrixTagNonlocal (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, DenseMatrix< Number > &k) const
 
void prepareMatrixTagLower (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)
 
void accumulateTaggedLocalResidual ()
 
void assignTaggedLocalResidual ()
 
void accumulateTaggedLocalMatrix ()
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, const DenseMatrix< Number > &k)
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, const DenseMatrix< Number > &k)
 
void accumulateTaggedNonlocalMatrix ()
 
void assignTaggedLocalMatrix ()
 
void addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const DenseVector< T > &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobian (Assembly &assembly, DenseMatrix< Real > &local_k, const std::vector< dof_id_type > &row_indices, const std::vector< dof_id_type > &column_indices, Real scaling_factor)
 
void addResidualsWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianElement (Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)
 
void setResidual (SystemBase &sys, const T &residual, MooseVariableFE< T > &var)
 
void setResidual (SystemBase &sys, Real residual, dof_id_type dof_index)
 
void setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
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
 
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
 
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
 
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
 
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 T & coupledNodalValue (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 T & coupledNodalValueOld (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
const T & coupledNodalDot (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)
 
const T * getVarHelper (const std::string &var_name, unsigned int comp) const
 
T * getVarHelper (const std::string &var_name, unsigned int comp)
 
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< T > coupledVectorHelper (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
 
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
 
void addMooseVariableDependency (MooseVariableFieldBase *var)
 
void addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
 
virtual void checkMaterialProperty (const std::string &name, const unsigned int state)
 
void markMatPropRequested (const std::string &)
 
MaterialPropertyName getMaterialPropertyName (const std::string &name) const
 
void checkExecutionStage ()
 
virtual const OutputTools< Real >::VariableValuevalue ()
 
virtual const OutputTools< Real >::VariableValuevalueOld ()
 
virtual const OutputTools< Real >::VariableValuevalueOlder ()
 
virtual const OutputTools< Real >::VariableValuedot ()
 
virtual const OutputTools< Real >::VariableValuedotDot ()
 
virtual const OutputTools< Real >::VariableValuedotOld ()
 
virtual const OutputTools< Real >::VariableValuedotDotOld ()
 
virtual const VariableValuedotDu ()
 
virtual const VariableValuedotDotDu ()
 
virtual const OutputTools< Real >::VariableGradientgradient ()
 
virtual const OutputTools< Real >::VariableGradientgradientOld ()
 
virtual const OutputTools< Real >::VariableGradientgradientOlder ()
 
virtual const OutputTools< Real >::VariableSecondsecond ()
 
virtual const OutputTools< Real >::VariableSecondsecondOld ()
 
virtual const OutputTools< Real >::VariableSecondsecondOlder ()
 
virtual const OutputTools< Real >::VariableTestSecondsecondTest ()
 
virtual const OutputTools< Real >::VariableTestSecondsecondTestFace ()
 
virtual const OutputTools< Real >::VariablePhiSecondsecondPhi ()
 
virtual const OutputTools< Real >::VariablePhiSecondsecondPhiFace ()
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > * defaultGenericMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const MaterialProperty< T > * defaultMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > * defaultADMaterialProperty (const std::string &name)
 

Protected Attributes

const MaterialProperty< std::vector< Real > > & _enthalpy
 Enthalpy of each phase. More...
 
const MaterialProperty< std::vector< std::vector< Real > > > & _denthalpy_dvar
 Derivative of the enthalpy wrt PorousFlow variables. More...
 
const MaterialProperty< std::vector< Real > > & _relative_permeability
 Relative permeability of each phase. More...
 
const MaterialProperty< std::vector< std::vector< Real > > > & _drelative_permeability_dvar
 Derivative of relative permeability of each phase wrt PorousFlow variables. More...
 
const MaterialProperty< RealTensorValue > & _permeability
 Permeability of porous material. More...
 
const MaterialProperty< std::vector< RealTensorValue > > & _dpermeability_dvar
 d(permeabiity)/d(PorousFlow variable) More...
 
const MaterialProperty< std::vector< std::vector< RealTensorValue > > > & _dpermeability_dgradvar
 d(permeabiity)/d(grad(PorousFlow variable)) More...
 
const MaterialProperty< std::vector< Real > > & _fluid_density_node
 Fluid density for each phase (at the node) More...
 
const MaterialProperty< std::vector< std::vector< Real > > > & _dfluid_density_node_dvar
 Derivative of the fluid density for each phase wrt PorousFlow variables (at the node) More...
 
const MaterialProperty< std::vector< Real > > & _fluid_density_qp
 Fluid density for each phase (at the qp) More...
 
const MaterialProperty< std::vector< std::vector< Real > > > & _dfluid_density_qp_dvar
 Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp) More...
 
const MaterialProperty< std::vector< Real > > & _fluid_viscosity
 Viscosity of each component in each phase. More...
 
const MaterialProperty< std::vector< std::vector< Real > > > & _dfluid_viscosity_dvar
 Derivative of the fluid viscosity for each phase wrt PorousFlow variables. More...
 
const MaterialProperty< std::vector< Real > > & _pp
 Nodal pore pressure in each phase. More...
 
const MaterialProperty< std::vector< RealGradient > > & _grad_p
 Gradient of the pore pressure in each phase. More...
 
const MaterialProperty< std::vector< std::vector< Real > > > & _dgrad_p_dgrad_var
 Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables) More...
 
const MaterialProperty< std::vector< std::vector< RealGradient > > > & _dgrad_p_dvar
 Derivative of Grad porepressure in each phase wrt PorousFlow variables. More...
 
const PorousFlowDictator_dictator
 PorousFlowDictator UserObject. More...
 
const unsigned int _num_phases
 The number of fluid phases. More...
 
const RealVectorValue _gravity
 Gravity. Defaults to 9.81 m/s^2. More...
 
const bool _perm_derivs
 Flag to check whether permeabiity derivatives are non-zero. More...
 
const unsigned _full_upwind_threshold
 If the number of upwind-downwind swaps is less than this amount then full upwinding is used. More...
 
enum PorousFlowDarcyBase::FallbackEnum _fallback_scheme
 
std::vector< std::vector< Real > > _proto_flux
 The Darcy flux. More...
 
std::vector< std::vector< std::vector< Real > > > _jacobian
 Derivative of _proto_flux with respect to nodal variables. More...
 
std::unordered_map< unsigned, std::vector< std::vector< unsigned > > > _num_upwinds
 Number of nonlinear iterations (in this timestep and this element) that a node is an upwind node for a given fluid phase. More...
 
std::unordered_map< unsigned, std::vector< std::vector< unsigned > > > _num_downwinds
 Number of nonlinear iterations (in this timestep and this element) that a node is an downwind node for a given fluid phase. More...
 
MooseVariable_var
 
const VariableTestValue_test
 
const VariableTestGradient_grad_test
 
const VariablePhiValue_phi
 
const VariablePhiGradient_grad_phi
 
const VariableValue_u
 
const VariableGradient_grad_u
 
const Elem *const & _current_elem
 
const Real_current_elem_volume
 
unsigned int _qp
 
const MooseArray< Point > & _q_point
 
const QBase *const & _qrule
 
const MooseArray< Real > & _JxW
 
const MooseArray< Real > & _coord
 
unsigned int _i
 
unsigned int _j
 
bool _has_save_in
 
std::vector< MooseVariableFEBase *> _save_in
 
std::vector< AuxVariableName > _save_in_strings
 
bool _has_diag_save_in
 
std::vector< MooseVariableFEBase *> _diag_save_in
 
std::vector< AuxVariableName > _diag_save_in_strings
 
std::vector< unsigned int_displacements
 
const bool _use_displaced_mesh
 
SubProblem_subproblem
 
FEProblemBase_fe_problem
 
SystemBase_sys
 
THREAD_ID _tid
 
Assembly_assembly
 
MooseMesh_mesh
 
const bool & _enabled
 
MooseApp_app
 
const std::string _type
 
const std::string _name
 
const InputParameters_pars
 
Factory_factory
 
ActionFactory_action_factory
 
const ExecFlagEnum_execute_enum
 
const ExecFlagType_current_execute_flag
 
const InputParameters_ti_params
 
FEProblemBase_ti_feproblem
 
bool _is_implicit
 
Real_t
 
int_t_step
 
Real_dt
 
Real_dt_old
 
bool _is_transient
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
DenseVector< Number_local_re
 
DenseMatrix< Number_local_ke
 
DenseMatrix< Number_nonlocal_ke
 
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< MooseVariableFV< Real > *> _coupled_standard_fv_moose_vars
 
std::vector< MooseLinearVariableFV< Real > *> _coupled_standard_linear_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< DualReal > > > _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
 
const VariableValue_zero
 
const VariablePhiValue_phi_zero
 
const MooseArray< DualReal > & _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
 
ArrayVariableValue _default_array_value_zero
 
ArrayVariableGradient _default_array_gradient
 
ArrayVariableCurl _default_array_curl
 
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_mi_params
 
const std::string _mi_name
 
const MooseObjectName _mi_moose_object_name
 
FEProblemBase_mi_feproblem
 
SubProblem_mi_subproblem
 
const THREAD_ID _mi_tid
 
const Moose::MaterialDataType _material_data_type
 
MaterialData_material_data
 
bool _stateful_allowed
 
bool _get_material_property_called
 
std::vector< std::unique_ptr< PropertyValue > > _default_properties
 
std::unordered_set< unsigned int_material_property_dependencies
 
const MaterialPropertyName _get_suffix
 
const bool _use_interpolated_state
 
GeometricSearchData_geometric_search_data
 
bool _requires_geometric_search
 
bool _nodal
 
MooseVariableFE< Real > * _variable
 
MooseVariableFV< Real > * _fv_variable
 
MooseLinearVariableFV< Real > * _linear_fv_variable
 
MooseVariableField< Real > * _field_variable
 
Assembly_mvi_assembly
 
const Parallel::Communicator & _communicator
 

Static Protected Attributes

static const std::string _interpolated_old
 
static const std::string _interpolated_older
 

Detailed Description

Advection of heat via flux of component k in fluid phase alpha.

A fully-updwinded version is implemented, where the mobility of the upstream nodes is used.

Definition at line 19 of file PorousFlowHeatAdvection.h.

Member Enumeration Documentation

◆ FallbackEnum

enum PorousFlowDarcyBase::FallbackEnum
strongprotectedinherited

If full upwinding is failing due to nodes swapping between upwind and downwind in successive nonlinear iterations (within one timestep and element) a fallback scheme is used instead.

QUICK: use the nodal mobilities, as in full upwinding, but do not attempt to conserve fluid mass. HARMONIC: assign the harmonic mean of the mobilities to the entire element.

Enumerator
QUICK 
HARMONIC 

Definition at line 150 of file PorousFlowDarcyBase.h.

150 { QUICK, HARMONIC } _fallback_scheme;
enum PorousFlowDarcyBase::FallbackEnum _fallback_scheme

◆ JacRes

enum PorousFlowDarcyBase::JacRes
strongprotectedinherited
Enumerator
CALCULATE_RESIDUAL 
CALCULATE_JACOBIAN 

Definition at line 59 of file PorousFlowDarcyBase.h.

60  {
61  CALCULATE_RESIDUAL = 0,
62  CALCULATE_JACOBIAN = 1
63  };

Constructor & Destructor Documentation

◆ PorousFlowHeatAdvection()

PorousFlowHeatAdvection::PorousFlowHeatAdvection ( const InputParameters parameters)

Definition at line 22 of file PorousFlowHeatAdvection.C.

24  _enthalpy(getMaterialProperty<std::vector<Real>>("PorousFlow_fluid_phase_enthalpy_nodal")),
25  _denthalpy_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
26  "dPorousFlow_fluid_phase_enthalpy_nodal_dvar")),
28  getMaterialProperty<std::vector<Real>>("PorousFlow_relative_permeability_nodal")),
29  _drelative_permeability_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
30  "dPorousFlow_relative_permeability_nodal_dvar"))
31 {
32 }
const MaterialProperty< std::vector< Real > > & _enthalpy
Enthalpy of each phase.
const MaterialProperty< std::vector< std::vector< Real > > > & _drelative_permeability_dvar
Derivative of relative permeability of each phase wrt PorousFlow variables.
const MaterialProperty< std::vector< Real > > & _relative_permeability
Relative permeability of each phase.
PorousFlowDarcyBase(const InputParameters &parameters)
const MaterialProperty< T > & getMaterialProperty(const std::string &name, MaterialData &material_data, const unsigned int state=0)
const InputParameters & parameters() const
const MaterialProperty< std::vector< std::vector< Real > > > & _denthalpy_dvar
Derivative of the enthalpy wrt PorousFlow variables.

Member Function Documentation

◆ computeJacobian()

void PorousFlowDarcyBase::computeJacobian ( )
overrideprotectedvirtualinherited

Reimplemented from Kernel.

Definition at line 140 of file PorousFlowDarcyBase.C.

141 {
143 }
MooseVariable & _var
unsigned int number() const
virtual void computeResidualAndJacobian() override

◆ computeOffDiagJacobian()

void PorousFlowDarcyBase::computeOffDiagJacobian ( unsigned int  jvar)
overrideprotectedvirtualinherited

Reimplemented from Kernel.

Definition at line 146 of file PorousFlowDarcyBase.C.

◆ computeQpResidual()

Real PorousFlowDarcyBase::computeQpResidual ( )
overrideprotectedvirtualinherited

Implements Kernel.

Definition at line 127 of file PorousFlowDarcyBase.C.

128 {
129  mooseError("PorousFlowDarcyBase: computeQpResidual called");
130  return 0.0;
131 }
void mooseError(Args &&... args) const

◆ computeResidual()

void PorousFlowDarcyBase::computeResidual ( )
overrideprotectedvirtualinherited

Reimplemented from Kernel.

Definition at line 134 of file PorousFlowDarcyBase.C.

◆ computeResidualAndJacobian() [1/2]

virtual void ResidualObject::computeResidualAndJacobian
protectedinherited

◆ computeResidualAndJacobian() [2/2]

void PorousFlowDarcyBase::computeResidualAndJacobian ( JacRes  res_or_jac,
unsigned int  jvar 
)
protectedinherited

Computation of the residual and Jacobian.

If res_or_jac=CALCULATE_JACOBIAN then the residual gets calculated anyway (becuase we have to know which nodes are upwind and which are downwind) except we don't actually need to store the residual in _assembly.residualBlock and we don't need to save in save_in.

If res_or_jac=CALCULATE_RESIDUAL then we don't have to do all the Jacobian calculations.

Parameters
res_or_jacWhether to calculate the residual or the jacobian.
jvarPorousFlow variable to take derivatives wrt to

Definition at line 152 of file PorousFlowDarcyBase.C.

153 {
154  if ((res_or_jac == JacRes::CALCULATE_JACOBIAN) && _dictator.notPorousFlowVariable(jvar))
155  return;
156 
157  // The PorousFlow variable index corresponding to the variable number jvar
158  const unsigned int pvar =
159  ((res_or_jac == JacRes::CALCULATE_JACOBIAN) ? _dictator.porousFlowVariableNum(jvar) : 0);
160 
162  if ((_local_ke.n() == 0) && (res_or_jac == JacRes::CALCULATE_JACOBIAN)) // this removes a problem
163  // encountered in the
164  // initial timestep when
165  // use_displaced_mesh=true
166  return;
167 
168  // The number of nodes in the element
169  const unsigned int num_nodes = _test.size();
170 
171  // Compute the residual and jacobian without the mobility terms. Even if we are computing the
172  // Jacobian we still need this in order to see which nodes are upwind and which are downwind.
173  for (unsigned ph = 0; ph < _num_phases; ++ph)
174  {
175  _proto_flux[ph].assign(num_nodes, 0);
176  for (_qp = 0; _qp < _qrule->n_points(); _qp++)
177  {
178  for (_i = 0; _i < num_nodes; ++_i)
179  _proto_flux[ph][_i] += _JxW[_qp] * _coord[_qp] * darcyQp(ph);
180  }
181  }
182 
183  // for this element, record whether each node is "upwind" or "downwind" (or neither)
184  const unsigned elem = _current_elem->id();
185  if (_num_upwinds.find(elem) == _num_upwinds.end())
186  {
187  _num_upwinds[elem] = std::vector<std::vector<unsigned>>(_num_phases);
188  _num_downwinds[elem] = std::vector<std::vector<unsigned>>(_num_phases);
189  for (unsigned ph = 0; ph < _num_phases; ++ph)
190  {
191  _num_upwinds[elem][ph].assign(num_nodes, 0);
192  _num_downwinds[elem][ph].assign(num_nodes, 0);
193  }
194  }
195  // record the information once per nonlinear iteration
196  if (res_or_jac == JacRes::CALCULATE_JACOBIAN && jvar == _var.number())
197  {
198  for (unsigned ph = 0; ph < _num_phases; ++ph)
199  {
200  for (unsigned nod = 0; nod < num_nodes; ++nod)
201  {
202  if (_proto_flux[ph][nod] > 0)
203  _num_upwinds[elem][ph][nod]++;
204  else if (_proto_flux[ph][nod] < 0)
205  _num_downwinds[elem][ph][nod]++;
206  }
207  }
208  }
209 
210  // based on _num_upwinds and _num_downwinds, calculate the maximum number
211  // of upwind-downwind swaps that have been encountered in this timestep
212  // for this element
213  std::vector<unsigned> max_swaps(_num_phases, 0);
214  for (unsigned ph = 0; ph < _num_phases; ++ph)
215  {
216  for (unsigned nod = 0; nod < num_nodes; ++nod)
217  max_swaps[ph] = std::max(
218  max_swaps[ph], std::min(_num_upwinds[elem][ph][nod], _num_downwinds[elem][ph][nod]));
219  }
220 
221  // size the _jacobian correctly and calculate it for the case residual = _proto_flux
222  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
223  {
224  for (unsigned ph = 0; ph < _num_phases; ++ph)
225  {
226  _jacobian[ph].resize(_local_ke.m());
227  for (_i = 0; _i < _test.size(); _i++)
228  {
229  _jacobian[ph][_i].assign(_local_ke.n(), 0.0);
230  for (_j = 0; _j < _phi.size(); _j++)
231  for (_qp = 0; _qp < _qrule->n_points(); _qp++)
232  _jacobian[ph][_i][_j] += _JxW[_qp] * _coord[_qp] * darcyQpJacobian(jvar, ph);
233  }
234  }
235  }
236 
237  // Loop over all the phases, computing the mass flux, which
238  // gets placed into _proto_flux, and the derivative of this
239  // which gets placed into _jacobian
240  for (unsigned int ph = 0; ph < _num_phases; ++ph)
241  {
242  if (max_swaps[ph] < _full_upwind_threshold)
243  fullyUpwind(res_or_jac, ph, pvar);
244  else
245  {
246  switch (_fallback_scheme)
247  {
248  case FallbackEnum::QUICK:
249  quickUpwind(res_or_jac, ph, pvar);
250  break;
252  harmonicMean(res_or_jac, ph, pvar);
253  break;
254  }
255  }
256  }
257 
258  // Add results to the Residual or Jacobian
259  if (res_or_jac == JacRes::CALCULATE_RESIDUAL)
260  {
262  for (_i = 0; _i < _test.size(); _i++)
263  for (unsigned int ph = 0; ph < _num_phases; ++ph)
264  _local_re(_i) += _proto_flux[ph][_i];
266 
267  if (_has_save_in)
268  for (unsigned int i = 0; i < _save_in.size(); i++)
269  _save_in[i]->sys().solution().add_vector(_local_re, _save_in[i]->dofIndices());
270  }
271 
272  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
273  {
274  for (_i = 0; _i < _test.size(); _i++)
275  for (_j = 0; _j < _phi.size(); _j++)
276  for (unsigned int ph = 0; ph < _num_phases; ++ph)
277  _local_ke(_i, _j) += _jacobian[ph][_i][_j];
278 
280 
281  if (_has_diag_save_in && jvar == _var.number())
282  {
283  unsigned int rows = _local_ke.m();
284  DenseVector<Number> diag(rows);
285  for (unsigned int i = 0; i < rows; i++)
286  diag(i) = _local_ke(i, i);
287 
288  for (unsigned int i = 0; i < _diag_save_in.size(); i++)
289  _diag_save_in[i]->sys().solution().add_vector(diag, _diag_save_in[i]->dofIndices());
290  }
291  }
292 }
bool notPorousFlowVariable(unsigned int moose_var_num) const
Returns true if moose_var_num is not a porous flow variabe.
void accumulateTaggedLocalResidual()
MooseVariable & _var
std::vector< MooseVariableFEBase *> _save_in
const MooseArray< Real > & _JxW
enum PorousFlowDarcyBase::FallbackEnum _fallback_scheme
unsigned int number() const
std::unordered_map< unsigned, std::vector< std::vector< unsigned > > > _num_downwinds
Number of nonlinear iterations (in this timestep and this element) that a node is an downwind node fo...
std::vector< std::vector< Real > > _proto_flux
The Darcy flux.
std::vector< std::vector< std::vector< Real > > > _jacobian
Derivative of _proto_flux with respect to nodal variables.
virtual Real darcyQpJacobian(unsigned int jvar, unsigned int ph) const
Jacobian of the Darcy part of the flux.
unsigned int m() const
bool _has_diag_save_in
DenseMatrix< Number > _local_ke
const VariableTestValue & _test
std::vector< MooseVariableFEBase *> _diag_save_in
const PorousFlowDictator & _dictator
PorousFlowDictator UserObject.
const unsigned int _num_phases
The number of fluid phases.
const QBase *const & _qrule
bool _has_save_in
unsigned int _i
std::unordered_map< unsigned, std::vector< std::vector< unsigned > > > _num_upwinds
Number of nonlinear iterations (in this timestep and this element) that a node is an upwind node for ...
void accumulateTaggedLocalMatrix()
void harmonicMean(JacRes res_or_jac, unsigned int ph, unsigned int pvar)
Calculate the residual or Jacobian by using the harmonic mean of the nodal mobilities for the entire ...
const MooseArray< Real > & _coord
Assembly & _assembly
unsigned int _j
void quickUpwind(JacRes res_or_jac, unsigned int ph, unsigned int pvar)
Calculate the residual or Jacobian using the nodal mobilities, but without conserving fluid mass...
const unsigned _full_upwind_threshold
If the number of upwind-downwind swaps is less than this amount then full upwinding is used...
DenseVector< Number > _local_re
const Elem *const & _current_elem
unsigned int porousFlowVariableNum(unsigned int moose_var_num) const
The PorousFlow variable number.
unsigned int n() const
void prepareVectorTag(Assembly &assembly, unsigned int ivar)
void prepareMatrixTag(Assembly &assembly, unsigned int ivar, unsigned int jvar)
virtual Real darcyQp(unsigned int ph) const
The Darcy part of the flux (this is the non-upwinded part)
const VariablePhiValue & _phi
unsigned int _qp
void fullyUpwind(JacRes res_or_jac, unsigned int ph, unsigned int pvar)
Calculate the residual or Jacobian using full upwinding.

◆ darcyQp()

Real PorousFlowDarcyBase::darcyQp ( unsigned int  ph) const
protectedvirtualinherited

The Darcy part of the flux (this is the non-upwinded part)

Definition at line 94 of file PorousFlowDarcyBase.C.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

95 {
96  return _grad_test[_i][_qp] *
98 }
const MaterialProperty< RealTensorValue > & _permeability
Permeability of porous material.
const MaterialProperty< std::vector< Real > > & _fluid_density_qp
Fluid density for each phase (at the qp)
const MaterialProperty< std::vector< RealGradient > > & _grad_p
Gradient of the pore pressure in each phase.
unsigned int _i
const RealVectorValue _gravity
Gravity. Defaults to 9.81 m/s^2.
const VariableTestGradient & _grad_test
unsigned int _qp

◆ darcyQpJacobian()

Real PorousFlowDarcyBase::darcyQpJacobian ( unsigned int  jvar,
unsigned int  ph 
) const
protectedvirtualinherited

Jacobian of the Darcy part of the flux.

Definition at line 101 of file PorousFlowDarcyBase.C.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

102 {
104  return 0.0;
105 
106  const unsigned int pvar = _dictator.porousFlowVariableNum(jvar);
107 
109  _permeability[_qp] * (_grad_phi[_j][_qp] * _dgrad_p_dgrad_var[_qp][ph][pvar] -
110  _phi[_j][_qp] * _dfluid_density_qp_dvar[_qp][ph][pvar] * _gravity);
111 
112  deriv += _permeability[_qp] * (_dgrad_p_dvar[_qp][ph][pvar] * _phi[_j][_qp]);
113 
114  if (_perm_derivs)
115  {
116  deriv += _dpermeability_dvar[_qp][pvar] * _phi[_j][_qp] *
117  (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity);
118  for (const auto i : make_range(Moose::dim))
119  deriv += _dpermeability_dgradvar[_qp][i][pvar] * _grad_phi[_j][_qp](i) *
120  (_grad_p[_qp][ph] - _fluid_density_qp[_qp][ph] * _gravity);
121  }
122 
123  return _grad_test[_i][_qp] * deriv;
124 }
const MaterialProperty< std::vector< std::vector< RealTensorValue > > > & _dpermeability_dgradvar
d(permeabiity)/d(grad(PorousFlow variable))
const MaterialProperty< std::vector< std::vector< Real > > > & _dfluid_density_qp_dvar
Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp) ...
const MaterialProperty< std::vector< RealTensorValue > > & _dpermeability_dvar
d(permeabiity)/d(PorousFlow variable)
bool notPorousFlowVariable(unsigned int moose_var_num) const
Returns true if moose_var_num is not a porous flow variabe.
const VariablePhiGradient & _grad_phi
const MaterialProperty< RealTensorValue > & _permeability
Permeability of porous material.
static constexpr std::size_t dim
const MaterialProperty< std::vector< Real > > & _fluid_density_qp
Fluid density for each phase (at the qp)
const MaterialProperty< std::vector< RealGradient > > & _grad_p
Gradient of the pore pressure in each phase.
const MaterialProperty< std::vector< std::vector< RealGradient > > > & _dgrad_p_dvar
Derivative of Grad porepressure in each phase wrt PorousFlow variables.
Real deriv(unsigned n, unsigned alpha, unsigned beta, Real x)
const PorousFlowDictator & _dictator
PorousFlowDictator UserObject.
unsigned int _i
const RealVectorValue _gravity
Gravity. Defaults to 9.81 m/s^2.
unsigned int _j
const MaterialProperty< std::vector< std::vector< Real > > > & _dgrad_p_dgrad_var
Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables)
const VariableTestGradient & _grad_test
IntRange< T > make_range(T beg, T end)
unsigned int porousFlowVariableNum(unsigned int moose_var_num) const
The PorousFlow variable number.
const VariablePhiValue & _phi
const bool _perm_derivs
Flag to check whether permeabiity derivatives are non-zero.
unsigned int _qp

◆ dmobility()

Real PorousFlowHeatAdvection::dmobility ( unsigned  nodenum,
unsigned  phase,
unsigned  pvar 
) const
overrideprotectedvirtual

The derivative of mobility with respect to PorousFlow variable pvar.

Parameters
nodenumThe node-number to evaluate the mobility for
phasethe fluid phase number
pvarthe PorousFlow variable pvar

Reimplemented from PorousFlowDarcyBase.

Definition at line 42 of file PorousFlowHeatAdvection.C.

43 {
44  Real dm = _denthalpy_dvar[nodenum][phase][pvar] * _fluid_density_node[nodenum][phase] *
45  _relative_permeability[nodenum][phase] / _fluid_viscosity[nodenum][phase];
46  dm += _enthalpy[nodenum][phase] * _dfluid_density_node_dvar[nodenum][phase][pvar] *
47  _relative_permeability[nodenum][phase] / _fluid_viscosity[nodenum][phase];
48  dm += _enthalpy[nodenum][phase] * _fluid_density_node[nodenum][phase] *
49  _drelative_permeability_dvar[nodenum][phase][pvar] / _fluid_viscosity[nodenum][phase];
50  dm -= _enthalpy[nodenum][phase] * _fluid_density_node[nodenum][phase] *
51  _relative_permeability[nodenum][phase] * _dfluid_viscosity_dvar[nodenum][phase][pvar] /
52  std::pow(_fluid_viscosity[nodenum][phase], 2);
53  return dm;
54 }
const MaterialProperty< std::vector< std::vector< Real > > > & _dfluid_density_node_dvar
Derivative of the fluid density for each phase wrt PorousFlow variables (at the node) ...
const MaterialProperty< std::vector< Real > > & _fluid_viscosity
Viscosity of each component in each phase.
const MaterialProperty< std::vector< Real > > & _enthalpy
Enthalpy of each phase.
const MaterialProperty< std::vector< std::vector< Real > > > & _drelative_permeability_dvar
Derivative of relative permeability of each phase wrt PorousFlow variables.
const MaterialProperty< std::vector< std::vector< Real > > > & _dfluid_viscosity_dvar
Derivative of the fluid viscosity for each phase wrt PorousFlow variables.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const MaterialProperty< std::vector< Real > > & _relative_permeability
Relative permeability of each phase.
MooseUnits pow(const MooseUnits &, int)
const MaterialProperty< std::vector< Real > > & _fluid_density_node
Fluid density for each phase (at the node)
const MaterialProperty< std::vector< std::vector< Real > > > & _denthalpy_dvar
Derivative of the enthalpy wrt PorousFlow variables.

◆ fullyUpwind()

void PorousFlowDarcyBase::fullyUpwind ( JacRes  res_or_jac,
unsigned int  ph,
unsigned int  pvar 
)
protectedinherited

Calculate the residual or Jacobian using full upwinding.

Parameters
res_or_jacwhether to compute the residual or jacobian
phfluid phase number
pvardifferentiate wrt to this PorousFlow variable (when computing the jacobian)

Perform the full upwinding by multiplying the residuals at the upstream nodes by their mobilities. Mobility is different for each phase, and in each situation: mobility = density / viscosity for single-component Darcy flow mobility = mass_fraction * density * relative_perm / viscosity for multi-component, multiphase flow mobility = enthalpy * density * relative_perm / viscosity for heat convection

The residual for the kernel is the sum over Darcy fluxes for each phase. The Darcy flux for a particular phase is R_i = int{mobility*flux_no_mob} = int{mobility*grad(pot)*permeability*grad(test_i)} for node i. where int is the integral over the element. However, in fully-upwind, the first step is to take the mobility outside the integral, which was done in the _proto_flux calculation above.

NOTE: Physically _proto_flux[i][ph] is a measure of fluid of phase ph flowing out of node i. If we had left in mobility, it would be exactly the component mass flux flowing out of node i.

This leads to the definition of upwinding:

If _proto_flux(i)[ph] is positive then we use mobility_i. That is we use the upwind value of mobility.

The final subtle thing is we must also conserve fluid mass: the total component mass flowing out of node i must be the sum of the masses flowing into the other nodes.

note the -= means the result is positive

Definition at line 295 of file PorousFlowDarcyBase.C.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

296 {
326  // The number of nodes in the element
327  const unsigned int num_nodes = _test.size();
328 
329  Real mob;
330  Real dmob;
331  // Define variables used to ensure mass conservation
332  Real total_mass_out = 0.0;
333  Real total_in = 0.0;
334 
335  // The following holds derivatives of these
336  std::vector<Real> dtotal_mass_out;
337  std::vector<Real> dtotal_in;
338  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
339  {
340  dtotal_mass_out.assign(num_nodes, 0.0);
341  dtotal_in.assign(num_nodes, 0.0);
342  }
343 
344  // Perform the upwinding using the mobility
345  std::vector<bool> upwind_node(num_nodes);
346  for (unsigned int n = 0; n < num_nodes; ++n)
347  {
348  if (_proto_flux[ph][n] >= 0.0) // upstream node
349  {
350  upwind_node[n] = true;
351  // The mobility at the upstream node
352  mob = mobility(n, ph);
353  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
354  {
355  // The derivative of the mobility wrt the PorousFlow variable
356  dmob = dmobility(n, ph, pvar);
357 
358  for (_j = 0; _j < _phi.size(); _j++)
359  _jacobian[ph][n][_j] *= mob;
360 
361  if (_test.size() == _phi.size())
362  /* mobility at node=n depends only on the variables at node=n, by construction. For
363  * linear-lagrange variables, this means that Jacobian entries involving the derivative
364  * of mobility will only be nonzero for derivatives wrt variables at node=n. Hence the
365  * [n][n] in the line below. However, for other variable types (eg constant monomials)
366  * I cannot tell what variable number contributes to the derivative. However, in all
367  * cases I can possibly imagine, the derivative is zero anyway, since in the full
368  * upwinding scheme, mobility shouldn't depend on these other sorts of variables.
369  */
370  _jacobian[ph][n][n] += dmob * _proto_flux[ph][n];
371 
372  for (_j = 0; _j < _phi.size(); _j++)
373  dtotal_mass_out[_j] += _jacobian[ph][n][_j];
374  }
375  _proto_flux[ph][n] *= mob;
376  total_mass_out += _proto_flux[ph][n];
377  }
378  else
379  {
380  upwind_node[n] = false;
381  total_in -= _proto_flux[ph][n];
382  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
383  for (_j = 0; _j < _phi.size(); _j++)
384  dtotal_in[_j] -= _jacobian[ph][n][_j];
385  }
386  }
387 
388  // Conserve mass over all phases by proportioning the total_mass_out mass to the inflow nodes,
389  // weighted by their proto_flux values
390  for (unsigned int n = 0; n < num_nodes; ++n)
391  {
392  if (!upwind_node[n]) // downstream node
393  {
394  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
395  for (_j = 0; _j < _phi.size(); _j++)
396  {
397  _jacobian[ph][n][_j] *= total_mass_out / total_in;
398  _jacobian[ph][n][_j] +=
399  _proto_flux[ph][n] * (dtotal_mass_out[_j] / total_in -
400  dtotal_in[_j] * total_mass_out / total_in / total_in);
401  }
402  _proto_flux[ph][n] *= total_mass_out / total_in;
403  }
404  }
405 }
std::vector< std::vector< Real > > _proto_flux
The Darcy flux.
std::vector< std::vector< std::vector< Real > > > _jacobian
Derivative of _proto_flux with respect to nodal variables.
const VariableTestValue & _test
unsigned int _j
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual Real mobility(unsigned nodenum, unsigned phase) const
The mobility of the fluid.
virtual Real dmobility(unsigned nodenum, unsigned phase, unsigned pvar) const
The derivative of mobility with respect to PorousFlow variable pvar.
const VariablePhiValue & _phi

◆ harmonicMean()

void PorousFlowDarcyBase::harmonicMean ( JacRes  res_or_jac,
unsigned int  ph,
unsigned int  pvar 
)
protectedinherited

Calculate the residual or Jacobian by using the harmonic mean of the nodal mobilities for the entire element.

Parameters
res_or_jacwhether to compute the residual or jacobian
phfluid phase number
pvardifferentiate wrt to this PorousFlow variable (when computing the jacobian)

Definition at line 445 of file PorousFlowDarcyBase.C.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

446 {
447  // The number of nodes in the element
448  const unsigned int num_nodes = _test.size();
449 
450  std::vector<Real> mob(num_nodes);
451  unsigned num_zero = 0;
452  unsigned zero_mobility_node = std::numeric_limits<unsigned>::max();
453  Real harmonic_mob = 0;
454  for (unsigned n = 0; n < num_nodes; ++n)
455  {
456  mob[n] = mobility(n, ph);
457  if (mob[n] == 0.0)
458  {
459  zero_mobility_node = n;
460  num_zero++;
461  }
462  else
463  harmonic_mob += 1.0 / mob[n];
464  }
465  if (num_zero > 0)
466  harmonic_mob = 0.0;
467  else
468  harmonic_mob = (1.0 * num_nodes) / harmonic_mob;
469 
470  // d(harmonic_mob)/d(PorousFlow variable at node n)
471  std::vector<Real> dharmonic_mob(num_nodes, 0.0);
472  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
473  {
474  const Real harm2 = std::pow(harmonic_mob, 2) / (1.0 * num_nodes);
475  if (num_zero == 0)
476  for (unsigned n = 0; n < num_nodes; ++n)
477  dharmonic_mob[n] = dmobility(n, ph, pvar) * harm2 / std::pow(mob[n], 2);
478  else if (num_zero == 1)
479  dharmonic_mob[zero_mobility_node] =
480  num_nodes * dmobility(zero_mobility_node, ph, pvar); // other derivs are zero
481  // if num_zero > 1 then all dharmonic_mob = 0.0
482  }
483 
484  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
485  for (unsigned n = 0; n < num_nodes; ++n)
486  for (_j = 0; _j < _phi.size(); _j++)
487  {
488  _jacobian[ph][n][_j] *= harmonic_mob;
489  if (_test.size() == _phi.size())
490  _jacobian[ph][n][_j] += dharmonic_mob[_j] * _proto_flux[ph][n];
491  }
492 
493  if (res_or_jac == JacRes::CALCULATE_RESIDUAL)
494  for (unsigned n = 0; n < num_nodes; ++n)
495  _proto_flux[ph][n] *= harmonic_mob;
496 }
std::vector< std::vector< Real > > _proto_flux
The Darcy flux.
std::vector< std::vector< std::vector< Real > > > _jacobian
Derivative of _proto_flux with respect to nodal variables.
const VariableTestValue & _test
unsigned int _j
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual Real mobility(unsigned nodenum, unsigned phase) const
The mobility of the fluid.
virtual Real dmobility(unsigned nodenum, unsigned phase, unsigned pvar) const
The derivative of mobility with respect to PorousFlow variable pvar.
const VariablePhiValue & _phi
MooseUnits pow(const MooseUnits &, int)

◆ mobility()

Real PorousFlowHeatAdvection::mobility ( unsigned  nodenum,
unsigned  phase 
) const
overrideprotectedvirtual

The mobility of the fluid.

For multi-component Darcy flow this is mass_fraction * fluid_density * relative_permeability / fluid_viscosity

Parameters
nodenumThe node-number to evaluate the mobility for
phasethe fluid phase number

Reimplemented from PorousFlowDarcyBase.

Definition at line 35 of file PorousFlowHeatAdvection.C.

36 {
37  return _enthalpy[nodenum][phase] * _fluid_density_node[nodenum][phase] *
38  _relative_permeability[nodenum][phase] / _fluid_viscosity[nodenum][phase];
39 }
const MaterialProperty< std::vector< Real > > & _fluid_viscosity
Viscosity of each component in each phase.
const MaterialProperty< std::vector< Real > > & _enthalpy
Enthalpy of each phase.
const MaterialProperty< std::vector< Real > > & _relative_permeability
Relative permeability of each phase.
const MaterialProperty< std::vector< Real > > & _fluid_density_node
Fluid density for each phase (at the node)

◆ quickUpwind()

void PorousFlowDarcyBase::quickUpwind ( JacRes  res_or_jac,
unsigned int  ph,
unsigned int  pvar 
)
protectedinherited

Calculate the residual or Jacobian using the nodal mobilities, but without conserving fluid mass.

Parameters
res_or_jacwhether to compute the residual or jacobian
phfluid phase number
pvardifferentiate wrt to this PorousFlow variable (when computing the jacobian)

Definition at line 408 of file PorousFlowDarcyBase.C.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

409 {
410  // The number of nodes in the element
411  const unsigned int num_nodes = _test.size();
412 
413  Real mob;
414  Real dmob;
415 
416  // Use the raw nodal mobility
417  for (unsigned int n = 0; n < num_nodes; ++n)
418  {
419  // The mobility at the node
420  mob = mobility(n, ph);
421  if (res_or_jac == JacRes::CALCULATE_JACOBIAN)
422  {
423  // The derivative of the mobility wrt the PorousFlow variable
424  dmob = dmobility(n, ph, pvar);
425 
426  for (_j = 0; _j < _phi.size(); _j++)
427  _jacobian[ph][n][_j] *= mob;
428 
429  if (_test.size() == _phi.size())
430  /* mobility at node=n depends only on the variables at node=n, by construction. For
431  * linear-lagrange variables, this means that Jacobian entries involving the derivative
432  * of mobility will only be nonzero for derivatives wrt variables at node=n. Hence the
433  * [n][n] in the line below. However, for other variable types (eg constant monomials)
434  * I cannot tell what variable number contributes to the derivative. However, in all
435  * cases I can possibly imagine, the derivative is zero anyway, since in the full
436  * upwinding scheme, mobility shouldn't depend on these other sorts of variables.
437  */
438  _jacobian[ph][n][n] += dmob * _proto_flux[ph][n];
439  }
440  _proto_flux[ph][n] *= mob;
441  }
442 }
std::vector< std::vector< Real > > _proto_flux
The Darcy flux.
std::vector< std::vector< std::vector< Real > > > _jacobian
Derivative of _proto_flux with respect to nodal variables.
const VariableTestValue & _test
unsigned int _j
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual Real mobility(unsigned nodenum, unsigned phase) const
The mobility of the fluid.
virtual Real dmobility(unsigned nodenum, unsigned phase, unsigned pvar) const
The derivative of mobility with respect to PorousFlow variable pvar.
const VariablePhiValue & _phi

◆ timestepSetup()

void PorousFlowDarcyBase::timestepSetup ( )
overrideprotectedvirtualinherited

Reimplemented from Kernel.

Definition at line 86 of file PorousFlowDarcyBase.C.

87 {
89  _num_upwinds = std::unordered_map<unsigned, std::vector<std::vector<unsigned>>>();
90  _num_downwinds = std::unordered_map<unsigned, std::vector<std::vector<unsigned>>>();
91 }
std::unordered_map< unsigned, std::vector< std::vector< unsigned > > > _num_downwinds
Number of nonlinear iterations (in this timestep and this element) that a node is an downwind node fo...
virtual void timestepSetup()
std::unordered_map< unsigned, std::vector< std::vector< unsigned > > > _num_upwinds
Number of nonlinear iterations (in this timestep and this element) that a node is an upwind node for ...

◆ validParams()

InputParameters PorousFlowHeatAdvection::validParams ( )
static

Definition at line 15 of file PorousFlowHeatAdvection.C.

16 {
18  params.addClassDescription("Fully-upwinded heat flux, advected by the fluid");
19  return params;
20 }
static InputParameters validParams()
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _denthalpy_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowHeatAdvection::_denthalpy_dvar
protected

Derivative of the enthalpy wrt PorousFlow variables.

Definition at line 34 of file PorousFlowHeatAdvection.h.

Referenced by dmobility().

◆ _dfluid_density_node_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDarcyBase::_dfluid_density_node_dvar
protectedinherited

Derivative of the fluid density for each phase wrt PorousFlow variables (at the node)

Definition at line 97 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowAdvectiveFlux::dmobility(), dmobility(), PorousFlowFullySaturatedUpwindHeatAdvection::dmobility(), PorousFlowFullySaturatedAdvectiveFlux::dmobility(), and PorousFlowDarcyBase::dmobility().

◆ _dfluid_density_qp_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDarcyBase::_dfluid_density_qp_dvar
protectedinherited

Derivative of the fluid density for each phase wrt PorousFlow variables (at the qp)

Definition at line 103 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQpJacobian().

◆ _dfluid_viscosity_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDarcyBase::_dfluid_viscosity_dvar
protectedinherited

◆ _dgrad_p_dgrad_var

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowDarcyBase::_dgrad_p_dgrad_var
protectedinherited

Derivative of Grad porepressure in each phase wrt grad(PorousFlow variables)

Definition at line 118 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQpJacobian().

◆ _dgrad_p_dvar

const MaterialProperty<std::vector<std::vector<RealGradient> > >& PorousFlowDarcyBase::_dgrad_p_dvar
protectedinherited

Derivative of Grad porepressure in each phase wrt PorousFlow variables.

Definition at line 121 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQpJacobian().

◆ _dictator

const PorousFlowDictator& PorousFlowDarcyBase::_dictator
protectedinherited

◆ _dpermeability_dgradvar

const MaterialProperty<std::vector<std::vector<RealTensorValue> > >& PorousFlowDarcyBase::_dpermeability_dgradvar
protectedinherited

d(permeabiity)/d(grad(PorousFlow variable))

Definition at line 91 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQpJacobian().

◆ _dpermeability_dvar

const MaterialProperty<std::vector<RealTensorValue> >& PorousFlowDarcyBase::_dpermeability_dvar
protectedinherited

d(permeabiity)/d(PorousFlow variable)

Definition at line 88 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQpJacobian().

◆ _drelative_permeability_dvar

const MaterialProperty<std::vector<std::vector<Real> > >& PorousFlowHeatAdvection::_drelative_permeability_dvar
protected

Derivative of relative permeability of each phase wrt PorousFlow variables.

Definition at line 40 of file PorousFlowHeatAdvection.h.

Referenced by dmobility().

◆ _enthalpy

const MaterialProperty<std::vector<Real> >& PorousFlowHeatAdvection::_enthalpy
protected

Enthalpy of each phase.

Definition at line 31 of file PorousFlowHeatAdvection.h.

Referenced by dmobility(), and mobility().

◆ _fallback_scheme

enum PorousFlowDarcyBase::FallbackEnum PorousFlowDarcyBase::_fallback_scheme
protectedinherited

◆ _fluid_density_node

const MaterialProperty<std::vector<Real> >& PorousFlowDarcyBase::_fluid_density_node
protectedinherited

◆ _fluid_density_qp

const MaterialProperty<std::vector<Real> >& PorousFlowDarcyBase::_fluid_density_qp
protectedinherited

Fluid density for each phase (at the qp)

Definition at line 100 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQp(), and PorousFlowDarcyBase::darcyQpJacobian().

◆ _fluid_viscosity

const MaterialProperty<std::vector<Real> >& PorousFlowDarcyBase::_fluid_viscosity
protectedinherited

◆ _full_upwind_threshold

const unsigned PorousFlowDarcyBase::_full_upwind_threshold
protectedinherited

If the number of upwind-downwind swaps is less than this amount then full upwinding is used.

Otherwise the fallback scheme is employed

Definition at line 139 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

◆ _grad_p

const MaterialProperty<std::vector<RealGradient> >& PorousFlowDarcyBase::_grad_p
protectedinherited

Gradient of the pore pressure in each phase.

Definition at line 115 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQp(), and PorousFlowDarcyBase::darcyQpJacobian().

◆ _gravity

const RealVectorValue PorousFlowDarcyBase::_gravity
protectedinherited

Gravity. Defaults to 9.81 m/s^2.

Definition at line 130 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQp(), and PorousFlowDarcyBase::darcyQpJacobian().

◆ _jacobian

std::vector<std::vector<std::vector<Real> > > PorousFlowDarcyBase::_jacobian
protectedinherited

Derivative of _proto_flux with respect to nodal variables.

This gets modified with the mobility and its derivative during computation of MOOSE's Jacobian

Definition at line 165 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian(), PorousFlowDarcyBase::fullyUpwind(), PorousFlowDarcyBase::harmonicMean(), and PorousFlowDarcyBase::quickUpwind().

◆ _num_downwinds

std::unordered_map<unsigned, std::vector<std::vector<unsigned> > > PorousFlowDarcyBase::_num_downwinds
protectedinherited

Number of nonlinear iterations (in this timestep and this element) that a node is an downwind node for a given fluid phase.

_num_upwinds[element_number][phase][node_number_in_element]

Definition at line 179 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian(), and PorousFlowDarcyBase::timestepSetup().

◆ _num_phases

const unsigned int PorousFlowDarcyBase::_num_phases
protectedinherited

The number of fluid phases.

Definition at line 127 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian().

◆ _num_upwinds

std::unordered_map<unsigned, std::vector<std::vector<unsigned> > > PorousFlowDarcyBase::_num_upwinds
protectedinherited

Number of nonlinear iterations (in this timestep and this element) that a node is an upwind node for a given fluid phase.

_num_upwinds[element_number][phase][node_number_in_element]

Definition at line 172 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian(), and PorousFlowDarcyBase::timestepSetup().

◆ _perm_derivs

const bool PorousFlowDarcyBase::_perm_derivs
protectedinherited

Flag to check whether permeabiity derivatives are non-zero.

Definition at line 133 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQpJacobian().

◆ _permeability

const MaterialProperty<RealTensorValue>& PorousFlowDarcyBase::_permeability
protectedinherited

Permeability of porous material.

Definition at line 85 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::darcyQp(), and PorousFlowDarcyBase::darcyQpJacobian().

◆ _pp

const MaterialProperty<std::vector<Real> >& PorousFlowDarcyBase::_pp
protectedinherited

Nodal pore pressure in each phase.

Definition at line 112 of file PorousFlowDarcyBase.h.

◆ _proto_flux

std::vector<std::vector<Real> > PorousFlowDarcyBase::_proto_flux
protectedinherited

The Darcy flux.

When multiplied by the mobility, this is the mass flux of fluid moving out of a node. This multiplication occurs during the formation of the residual and Jacobian. _proto_flux[num_phases][num_nodes]

Definition at line 158 of file PorousFlowDarcyBase.h.

Referenced by PorousFlowDarcyBase::computeResidualAndJacobian(), PorousFlowDarcyBase::fullyUpwind(), PorousFlowDarcyBase::harmonicMean(), and PorousFlowDarcyBase::quickUpwind().

◆ _relative_permeability

const MaterialProperty<std::vector<Real> >& PorousFlowHeatAdvection::_relative_permeability
protected

Relative permeability of each phase.

Definition at line 37 of file PorousFlowHeatAdvection.h.

Referenced by dmobility(), and mobility().


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