Loading [MathJax]/extensions/tex2jax.js
https://mooseframework.inl.gov
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
INSFVMomentumAdvection Class Reference

An advection kernel that implements interpolation schemes specific to Navier-Stokes flow physics. More...

#include <INSFVMomentumAdvection.h>

Inheritance diagram for INSFVMomentumAdvection:
[legend]

Public Types

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

Public Member Functions

 INSFVMomentumAdvection (const InputParameters &params)
 
void gatherRCData (const Elem &) override final
 Should be a non-empty implementation if the residual object is a FVElementalKernel and introduces residuals that are a function of the velocity, e.g. More...
 
void gatherRCData (const FaceInfo &fi) override final
 Should be a non-empty implementation if the residual object is a FVFluxKernel and introduces residuals that are a function of the velocity, e.g. More...
 
void initialSetup () override
 
void computeResidual (const FaceInfo &fi) override final
 
void computeJacobian (const FaceInfo &fi) override final
 
void computeResidual () override
 
void computeJacobian () override
 
void computeResidualAndJacobian () override
 
void computeResidualAndJacobian (const FaceInfo &fi) override
 
const MooseVariableFV< Real > & variable () const override
 
bool hasFaceSide (const FaceInfo &fi, const bool fi_elem_side) const override
 
virtual void computeOffDiagJacobian (unsigned int)
 
virtual void computeOffDiagJacobianScalar (unsigned int)
 
virtual void computeNonlocalJacobian ()
 
virtual void computeNonlocalOffDiagJacobian (unsigned int)
 
const SubProblemsubProblem () const
 
virtual void prepareShapes (unsigned int var_num)
 
virtual std::set< std::string > additionalROVariables ()
 
virtual bool enabled () const
 
std::shared_ptr< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () 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 * queryParam (const std::string &name) 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 &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void timestepSetup ()
 
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 ()
 
Moose::StateArg determineState () const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
 
const VectorPostprocessorName & getVectorPostprocessorName (const std::string &param_name) const
 
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 (const std::set< SubdomainName > &names) const
 
bool hasBlocks (SubdomainID id) const
 
bool hasBlocks (const std::vector< SubdomainID > &ids) const
 
bool hasBlocks (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::vector< SubdomainID > &ids) const
 
bool hasBlockMaterialProperty (const std::string &prop_name)
 
const std::set< SubdomainID > & meshBlockIDs () const
 
virtual bool blockRestricted () const
 
virtual void checkVariable (const MooseVariableFieldBase &variable) const
 
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
 
virtual 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)
 
MooseVariableBasemooseVariableBase () const
 
MooseVariableField< Real > & mooseVariableField ()
 
MooseVariableFE< Real > * mooseVariable () const
 
MooseVariableFV< Real > * mooseVariableFV () const
 
MooseLinearVariableFV< Real > * mooseLinearVariableFV () const
 
virtual const VariableValuecoupledNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledNeighborValues (const std::string &var_name) const
 
std::vector< const VariableValue *> coupledNeighborValuesOld (const std::string &var_name) const
 
std::vector< const VariableValue *> coupledNeighborValuesOlder (const std::string &var_name) const
 
virtual const ADVariableValueadCoupledNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
const auto & coupledGenericNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
const auto & coupledGenericNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableValueadCoupledNeighborValueDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledNeighborValues (const std::string &var_name) const
 
virtual const ADVectorVariableValueadCoupledVectorNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual std::vector< const VariableGradient *> coupledNeighborGradients (const std::string &var_name) const
 
virtual const VariableGradientcoupledNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ADVariableGradientadCoupledNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayNeighborValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayNeighborGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayNeighborGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayNeighborGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledNeighborSecond (const std::string &var_name, unsigned int i=0) const
 
virtual const VariableValuecoupledNeighborDofValues (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborDofValuesOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledNeighborDofValuesOlder (const std::string &var_name, unsigned int comp=0) 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 ADVectorVariableCurlgetADDefaultCurl () 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)
 
void addMooseVariableDependency (MooseVariableFieldBase *var)
 
void addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
 
Moose::FaceArg makeFace (const FaceInfo &fi, const Moose::FV::LimiterType limiter_type, const bool elem_is_upwind, const bool correct_skewness=false, const Moose::StateArg *state_limiter=nullptr) const
 
Moose::FaceArg makeCDFace (const FaceInfo &fi, const bool correct_skewness=false) 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 GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialProperty (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getNeighborADMaterialProperty (const std::string &name)
 
const ADMaterialProperty< T > & getNeighborADMaterialProperty (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOld (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOlder (const std::string &name)
 
const MaterialProperty< T > & getNeighborMaterialPropertyOlder (const std::string &name)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const GenericMaterialProperty< T, is_ad > & getGenericNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const MaterialProperty< T > & getNeighborMaterialPropertyByName (const std::string &name, const unsigned int state=0)
 
const ADMaterialProperty< T > & getNeighborADMaterialPropertyByName (const std::string &name)
 
const ADMaterialProperty< T > & getNeighborADMaterialPropertyByName (const std::string &name)
 
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)
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

static InputParameters validParams ()
 
static InputParameters uniqueParams ()
 Parameters of this object that should be added to the NSFV action that are unique to this object. More...
 
static std::vector< std::string > listOfCommonParams ()
 
static std::string deduceFunctorName (const std::string &name, const InputParameters &params)
 
static void setRMParamsAdvection (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short conditional_extended_layers)
 
static void setRMParamsDiffusion (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short conditional_extended_layers)
 

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 Member Functions

virtual ADReal computeQpResidual () override final
 
virtual bool hasMaterialTimeDerivative () const override
 
virtual const Moose::FunctorBase< ADReal > & epsilon () const
 A virtual method that allows us to reuse all the code from free-flow for porous. More...
 
virtual void computeResidualsAndAData (const FaceInfo &fi)
 Helper method that computes the 'a' coefficients and AD residuals. More...
 
template<typename T >
void initialSetup (T &insfv_fk)
 setup all the boundary condition member information More...
 
bool skipForBoundary (const FaceInfo &fi) const override
 
ADRealVectorValue velocity () const
 
std::string deduceFunctorName (const std::string &name) const
 
virtual ADReal gradUDotNormal (const Moose::StateArg &time, const bool correct_skewness) const
 
const RealVectorValuenormal () const
 
bool onBoundary (const FaceInfo &fi) const
 
Moose::ElemArg elemArg (bool correct_skewness=false) const
 
Moose::ElemArg neighborArg (bool correct_skewness=false) const
 
Moose::FaceArg singleSidedFaceArg (const FaceInfo *fi=nullptr, Moose::FV::LimiterType limiter_type=Moose::FV::LimiterType::CentralDifference, bool correct_skewness=false, const Moose::StateArg *state_limiter=nullptr) const
 
bool avoidBoundary (const FaceInfo &fi) const
 
virtual void precalculateResidual ()
 
virtual void precalculateJacobian ()
 
virtual void precalculateOffDiagJacobian (unsigned int)
 
const MooseVariableFieldBasegetVariable (unsigned int jvar_num) const
 
virtual void addUserObjectDependencyHelper (const UserObject &) 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 addResiduals (Assembly &assembly, const ADResidualsPacket &packet)
 
void addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
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 ()
 
const Moose::Functor< T > & getFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctor (const std::string &name, THREAD_ID tid)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
bool isFunctor (const std::string &name) const
 
bool isFunctor (const std::string &name, const SubProblem &subproblem) const
 
Moose::ElemArg makeElemArg (const Elem *elem, bool correct_skewnewss=false) const
 
void checkFunctorSupportsSideIntegration (const std::string &name, bool qp_integration)
 
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 >::VariableValueneighborValue ()
 
const VectorVariableValueneighborValue ()
 
virtual const OutputTools< Real >::VariableValueneighborValueOld ()
 
const VectorVariableValueneighborValueOld ()
 
virtual const OutputTools< Real >::VariableValueneighborValueOlder ()
 
const VectorVariableValueneighborValueOlder ()
 
virtual const OutputTools< Real >::VariableGradientneighborGradient ()
 
virtual const OutputTools< Real >::VariableGradientneighborGradientOld ()
 
virtual const OutputTools< Real >::VariableGradientneighborGradientOlder ()
 
virtual const OutputTools< Real >::VariableSecondneighborSecond ()
 
virtual const OutputTools< Real >::VariableSecondneighborSecondOld ()
 
virtual const OutputTools< Real >::VariableSecondneighborSecondOlder ()
 
virtual const OutputTools< Real >::VariableTestSecondneighborSecondTest ()
 
virtual const OutputTools< Real >::VariablePhiSecondneighborSecondPhi ()
 
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 ()
 
virtual void coupledCallback (const std::string &, bool) const
 
virtual bool isCoupled (const std::string &var_name, unsigned int i=0) const
 
virtual bool isCoupledConstant (const std::string &var_name) const
 
unsigned int coupledComponents (const std::string &var_name) const
 
VariableName coupledName (const std::string &var_name, unsigned int comp=0) const
 
std::vector< VariableName > coupledNames (const std::string &var_name) const
 
virtual unsigned int coupled (const std::string &var_name, unsigned int comp=0) const
 
std::vector< unsigned intcoupledIndices (const std::string &var_name) const
 
virtual const VariableValuecoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValues (const std::string &var_name) const
 
std::vector< const VectorVariableValue *> coupledVectorValues (const std::string &var_name) const
 
const GenericVariableValue< is_ad > & coupledGenericValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericValue (const std::string &var_name, unsigned int comp) const
 
const GenericVectorVariableValue< is_ad > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVectorVariableValue< false > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
 
const GenericVectorVariableValue< true > & coupledGenericVectorValue (const std::string &var_name, unsigned int comp) const
 
std::vector< const GenericVariableValue< is_ad > *> coupledGenericValues (const std::string &var_name) const
 
std::vector< const GenericVariableValue< false > *> coupledGenericValues (const std::string &var_name) const
 
std::vector< const GenericVariableValue< true > *> coupledGenericValues (const std::string &var_name) const
 
const GenericVariableValue< is_ad > & coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDofValue (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< is_ad > & coupledGenericDot (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< is_ad > & coupledGenericDotDot (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericDotDot (const std::string &var_name, unsigned int comp) const
 
virtual const VariableValuecoupledValueLower (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledValues (const std::string &var_name) const
 
const ADVariableValueadCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVectorVariableValue *> adCoupledVectorValues (const std::string &var_name) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableValue *> coupledVectorTagValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue *> coupledVectorTagValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableValuecoupledVectorTagArrayValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const ArrayVariableValue *> coupledVectorTagArrayValues (const std::string &var_names, TagID tag) const
 
std::vector< const ArrayVariableValue *> coupledVectorTagArrayValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableGradientcoupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableGradientcoupledVectorTagGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableGradient *> coupledVectorTagGradients (const std::string &var_names, TagID tag) const
 
std::vector< const VariableGradient *> coupledVectorTagGradients (const std::string &var_names, const std::string &tag_name) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const ArrayVariableGradientcoupledVectorTagArrayGradient (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const ArrayVariableGradient *> coupledVectorTagArrayGradients (const std::string &var_names, TagID tag) const
 
std::vector< const ArrayVariableGradient *> coupledVectorTagArrayGradients (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
const ArrayVariableValuecoupledVectorTagArrayDofValue (const std::string &var_name, const std::string &tag_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledVectorTagDofValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue *> coupledVectorTagDofValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VariableValuecoupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const
 
virtual const VariableValuecoupledMatrixTagValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const
 
std::vector< const VariableValue *> coupledMatrixTagValues (const std::string &var_names, TagID tag) const
 
std::vector< const VariableValue *> coupledMatrixTagValues (const std::string &var_names, const std::string &tag_name) const
 
virtual const VectorVariableValuecoupledVectorValue (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValue (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ArrayVariableValue *> coupledArrayValues (const std::string &var_name) const
 
MooseWritableVariablewritableVariable (const std::string &var_name, unsigned int comp=0)
 
virtual VariableValuewritableCoupledValue (const std::string &var_name, unsigned int comp=0)
 
void checkWritableVar (MooseWritableVariable *var)
 
virtual const VariableValuecoupledValueOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOld (const std::string &var_name) const
 
std::vector< const VectorVariableValue *> coupledVectorValuesOld (const std::string &var_name) const
 
virtual const VariableValuecoupledValueOlder (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledValuesOlder (const std::string &var_name) const
 
virtual const VariableValuecoupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradients (const std::string &var_name) const
 
const ADVariableGradientadCoupledGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableGradientadCoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableGradient *> adCoupledGradients (const std::string &var_name) const
 
const GenericVariableGradient< is_ad > & coupledGenericGradient (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableGradient< false > & coupledGenericGradient (const std::string &var_name, unsigned int comp) const
 
const GenericVariableGradient< true > & coupledGenericGradient (const std::string &var_name, unsigned int comp) const
 
std::vector< const GenericVariableGradient< is_ad > *> coupledGenericGradients (const std::string &var_name) const
 
std::vector< const GenericVariableGradient< false > *> coupledGenericGradients (const std::string &var_name) const
 
std::vector< const GenericVariableGradient< true > *> coupledGenericGradients (const std::string &var_name) const
 
const ADVectorVariableGradientadCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableSecondadCoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableSecondadCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientOld (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableGradient *> coupledGradientsOld (const std::string &var_name) const
 
virtual const VariableGradientcoupledGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableGradientcoupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableGradientcoupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradient (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableGradientcoupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableCurlcoupledCurlOlder (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableCurladCoupledCurl (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDiv (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableDivergencecoupledDivOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecond (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondOlder (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableSecondcoupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const VariableValue *> coupledDots (const std::string &var_name) const
 
virtual const VariableValuecoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledDot (const std::string &var_name, unsigned int comp=0) const
 
std::vector< const ADVariableValue *> adCoupledDots (const std::string &var_name) const
 
const ADVariableValueadCoupledDotDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVectorVariableValueadCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VectorVariableValuecoupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const ArrayVariableValuecoupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDu (const std::string &var_name, unsigned int comp=0) const
 
virtual const VariableValuecoupledDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledArrayDotDu (const std::string &var_name, unsigned int comp=0) const
 
const 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
 
libMesh::Order coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< is_ad > & coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const GenericVariableValue< false > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const GenericVariableValue< true > & coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const
 
const VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (const std::string &var_name, unsigned int comp) const
 
const 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 Moose::Functor< ADReal > & _rho
 Density. More...
 
const bool _approximate_as
 Whether to approximately calculate the 'a' coefficients. More...
 
const Real _cs
 Characteristic speed. More...
 
std::unique_ptr< PiecewiseByBlockLambdaFunctor< ADReal > > _rho_u
 Our local momentum functor. More...
 
ADReal _ae = 0
 The a coefficient for the element. More...
 
ADReal _an = 0
 The a coefficient for the neighbor. More...
 
ADReal _elem_residual = 0
 The element residual. More...
 
ADReal _neighbor_residual = 0
 The neighbor residual. More...
 
const Moose::ConstantFunctor< ADReal_unity_functor {1}
 A unity functor used in the epsilon virtual method. More...
 
Moose::FV::InterpMethod _advected_interp_method
 The interpolation method to use for the advected quantity. More...
 
Moose::FV::InterpMethod _velocity_interp_method
 The interpolation method to use for the velocity. More...
 
const RhieChowInterpolatorBase_rc_vel_provider
 The Rhie-Chow user object that provides us with the velocity. More...
 
MooseVariableFV< Real > & _var
 
const unsigned int _qp
 
const ADVariableValue_u_elem
 
const ADVariableValue_u_neighbor
 
RealVectorValue _normal
 
const FaceInfo_face_info
 
FaceInfo::VarFaceNeighbors _face_type
 
std::unordered_set< BoundaryID_boundaries_to_force
 
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
 
const Real_t_old
 
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
 
MaterialData_neighbor_material_data
 
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
 
bool _nodal
 
MooseVariableFE< Real > * _variable
 
MooseVariableFV< Real > * _fv_variable
 
MooseLinearVariableFV< Real > * _linear_fv_variable
 
MooseVariableField< Real > * _field_variable
 
Assembly_mvi_assembly
 
bool _neighbor_nodal
 
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< ADReal > > > _ad_default_value
 
std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > > _default_vector_value
 
std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > > _default_array_value
 
std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > > _ad_default_vector_value
 
VariableValue _default_value_zero
 
VariableGradient _default_gradient
 
MooseArray< ADRealVectorValue_ad_default_gradient
 
MooseArray< ADRealTensorValue_ad_default_vector_gradient
 
VariableSecond _default_second
 
MooseArray< ADRealTensorValue_ad_default_second
 
MooseArray< ADRealVectorValue_ad_default_curl
 
const VariableValue_zero
 
const VariablePhiValue_phi_zero
 
const MooseArray< ADReal > & _ad_zero
 
const VariableGradient_grad_zero
 
const MooseArray< ADRealVectorValue > & _ad_grad_zero
 
const VariablePhiGradient_grad_phi_zero
 
const VariableSecond_second_zero
 
const MooseArray< ADRealTensorValue > & _ad_second_zero
 
const VariablePhiSecond_second_phi_zero
 
const VectorVariableValue_vector_zero
 
const VectorVariableCurl_vector_curl_zero
 
VectorVariableValue _default_vector_value_zero
 
VectorVariableGradient _default_vector_gradient
 
VectorVariableCurl _default_vector_curl
 
VectorVariableDivergence _default_div
 
ArrayVariableValue _default_array_value_zero
 
ArrayVariableGradient _default_array_gradient
 
bool _coupleable_neighbor
 
FEProblemBase_sc_fe_problem
 
const THREAD_ID _sc_tid
 
const Real_real_zero
 
const VariableValue_scalar_zero
 
const Point & _point_zero
 
const Parallel::Communicator & _communicator
 
std::set< BoundaryID_no_slip_wall_boundaries
 Boundary IDs with no slip walls. More...
 
std::set< BoundaryID_slip_wall_boundaries
 Boundary IDs with slip walls. More...
 
std::set< BoundaryID_free_surface_boundaries
 Boundary IDs for free surfaces. More...
 
std::set< BoundaryID_flow_boundaries
 Flow Boundary IDs. More...
 
std::set< BoundaryID_fully_developed_flow_boundaries
 Fully Developed Flow Boundary IDs. This is a subset of _flow_boundaries. More...
 
std::set< BoundaryID_symmetry_boundaries
 Symmetry Boundary IDs. More...
 
std::set< BoundaryID_all_boundaries
 All the BoundaryIDs covered by our different types of INSFVBCs. More...
 
RhieChowInterpolatorBase_rc_uo
 The Rhie Chow user object that is responsible for generating face velocities for advection terms. More...
 
const unsigned int _index
 index x|y|z More...
 

Static Protected Attributes

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

Detailed Description

An advection kernel that implements interpolation schemes specific to Navier-Stokes flow physics.

Definition at line 20 of file INSFVMomentumAdvection.h.

Constructor & Destructor Documentation

◆ INSFVMomentumAdvection()

INSFVMomentumAdvection::INSFVMomentumAdvection ( const InputParameters params)

Definition at line 47 of file INSFVMomentumAdvection.C.

48  : INSFVAdvectionKernel(params),
50  _rho(getFunctor<ADReal>(NS::density)),
51  _approximate_as(isParamValid("characteristic_speed")),
52  _cs(_approximate_as ? getParam<Real>("characteristic_speed") : 0)
53 {
54 }
static const std::string density
Definition: NS.h:33
bool isParamValid(const std::string &name) const
INSFVAdvectionKernel(const InputParameters &params)
const Real _cs
Characteristic speed.
const bool _approximate_as
Whether to approximately calculate the &#39;a&#39; coefficients.
const Moose::Functor< ADReal > & _rho
Density.

Member Function Documentation

◆ computeJacobian()

void INSFVMomentumAdvection::computeJacobian ( const FaceInfo fi)
finaloverridevirtual

Reimplemented from FVFluxKernel.

Definition at line 238 of file INSFVMomentumAdvection.C.

239 {
240  if (skipForBoundary(fi))
241  return;
242 
244 
247  {
248  mooseAssert(_var.dofIndices().size() == 1, "We're currently built to use CONSTANT MONOMIALS");
249 
251  std::array<ADReal, 1>{{_elem_residual}},
252  _var.dofIndices(),
253  _var.scalingFactor());
254  }
255 
258  {
260  (_var.dofIndices().size() == 0),
261  "If the variable is only defined on the neighbor hand side of the face, then that "
262  "means it should have no dof indices on the elem element. Conversely if "
263  "the variable is defined on both sides of the face, then it should have a non-zero "
264  "number of degrees of freedom on the elem element");
265 
266  mooseAssert(_var.dofIndicesNeighbor().size() == 1,
267  "We're currently built to use CONSTANT MONOMIALS");
268 
270  std::array<ADReal, 1>{{_neighbor_residual}},
272  _var.scalingFactor());
273  }
274 }
virtual const std::vector< dof_id_type > & dofIndicesNeighbor() const final
ADReal _neighbor_residual
The neighbor residual.
void addResidualsAndJacobian(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
bool skipForBoundary(const FaceInfo &fi) const override
FaceInfo::VarFaceNeighbors _face_type
Assembly & _assembly
ADReal _elem_residual
The element residual.
MooseVariableFV< Real > & _var
virtual void computeResidualsAndAData(const FaceInfo &fi)
Helper method that computes the &#39;a&#39; coefficients and AD residuals.
void scalingFactor(const std::vector< Real > &factor)
virtual const std::vector< dof_id_type > & dofIndices() const final

◆ computeQpResidual()

ADReal INSFVMomentumAdvection::computeQpResidual ( )
finaloverrideprotectedvirtual

Implements FVFluxKernel.

Definition at line 71 of file INSFVMomentumAdvection.C.

72 {
73  mooseError("Should never be called");
74 }
void mooseError(Args &&... args) const

◆ computeResidual()

void INSFVMomentumAdvection::computeResidual ( const FaceInfo fi)
finaloverridevirtual

Reimplemented from FVFluxKernel.

Definition at line 211 of file INSFVMomentumAdvection.C.

212 {
213  if (skipForBoundary(fi))
214  return;
215 
217 
220  {
221  // residual contribution of this kernel to the elem element
225  }
226 
229  {
230  // residual contribution of this kernel to the neighbor element
234  }
235 }
ADReal _neighbor_residual
The neighbor residual.
void accumulateTaggedLocalResidual()
unsigned int number() const
auto raw_value(const Eigen::Map< T > &in)
bool skipForBoundary(const FaceInfo &fi) const override
FaceInfo::VarFaceNeighbors _face_type
Assembly & _assembly
void prepareVectorTagNeighbor(Assembly &assembly, unsigned int ivar)
DenseVector< Number > _local_re
ADReal _elem_residual
The element residual.
MooseVariableFV< Real > & _var
void prepareVectorTag(Assembly &assembly, unsigned int ivar)
virtual void computeResidualsAndAData(const FaceInfo &fi)
Helper method that computes the &#39;a&#39; coefficients and AD residuals.

◆ computeResidualsAndAData()

void INSFVMomentumAdvection::computeResidualsAndAData ( const FaceInfo fi)
protectedvirtual

Helper method that computes the 'a' coefficients and AD residuals.

Reimplemented in WCNSFV2PMomentumAdvectionSlip.

Definition at line 77 of file INSFVMomentumAdvection.C.

Referenced by computeJacobian(), computeResidual(), and gatherRCData().

78 {
79  mooseAssert(!skipForBoundary(fi),
80  "We shouldn't get in here if we're supposed to skip for a boundary");
81 
82  _face_info = &fi;
83  _normal = fi.normal();
84  _face_type = fi.faceType(std::make_pair(_var.number(), _var.sys().number()));
85  const auto state = determineState();
86 
87  using namespace Moose::FV;
88 
89  const bool correct_skewness = _advected_interp_method == InterpMethod::SkewCorrectedAverage;
90  const bool on_boundary = onBoundary(fi);
91 
92  _elem_residual = 0, _neighbor_residual = 0, _ae = 0, _an = 0;
93 
94  const auto v_face = velocity();
95  const auto vdotn = _normal * v_face;
96 
97  if (on_boundary)
98  {
99  const auto ssf = singleSidedFaceArg();
100  const Elem * const sided_elem = ssf.face_side;
101  const auto dof_number = sided_elem->dof_number(_sys.number(), _var.number(), 0);
102  const auto rho_face = _rho(ssf, state);
103  const auto eps_face = epsilon()(ssf, state);
104  const auto u_face = _var(ssf, state);
105  const Real d_u_face_d_dof = u_face.derivatives()[dof_number];
106  const auto coeff = vdotn * rho_face / eps_face;
107 
108  if (sided_elem == fi.elemPtr())
109  {
110  _ae = coeff * d_u_face_d_dof;
111  _elem_residual = coeff * u_face;
112  if (_approximate_as)
113  _ae = _cs / fi.elem().n_sides() * rho_face / eps_face;
114  }
115  else
116  {
117  _an = -coeff * d_u_face_d_dof;
118  _neighbor_residual = -coeff * u_face;
119  if (_approximate_as)
120  _an = _cs / fi.neighborPtr()->n_sides() * rho_face / eps_face;
121  }
122  }
123  else // we are an internal fluid flow face
124  {
125  const bool elem_is_upwind = MetaPhysicL::raw_value(v_face) * _normal > 0;
126  const auto & limiter_time = _subproblem.isTransient()
128  : Moose::StateArg(1, Moose::SolutionIterationType::Nonlinear);
129  const Moose::FaceArg advected_face_arg{&fi,
131  elem_is_upwind,
132  correct_skewness,
133  nullptr,
134  &limiter_time};
135  if (const auto [is_jump, eps_elem_face, eps_neighbor_face] =
136  NS::isPorosityJumpFace(epsilon(), fi, state);
137  is_jump)
138  {
139  // For a weakly compressible formulation, the density should not depend on pressure and
140  // consequently the density should not be impacted by the pressure jump that occurs at a
141  // porosity jump. Consequently we will allow evaluation of the density using both upstream and
142  // downstream information
143  const auto rho_face = _rho(advected_face_arg, state);
144 
145  // We set the + and - sides of the superficial velocity equal to the interpolated value
146  const auto & var_elem_face = v_face(_index);
147  const auto & var_neighbor_face = v_face(_index);
148 
149  const auto elem_dof_number = fi.elem().dof_number(_sys.number(), _var.number(), 0);
150  const auto neighbor_dof_number = fi.neighbor().dof_number(_sys.number(), _var.number(), 0);
151 
152  const auto d_var_elem_face_d_elem_dof = var_elem_face.derivatives()[elem_dof_number];
153  const auto d_var_neighbor_face_d_neighbor_dof =
154  var_neighbor_face.derivatives()[neighbor_dof_number];
155 
156  // We allow a discontintuity in the advective momentum flux at the jump face. Mainly the
157  // advective flux is:
158  // elem side:
159  // rho * v_superficial / eps_elem * v_superficial = rho * v_interstitial_elem * v_superficial
160  // neighbor side:
161  // rho * v_superficial / eps_neigh * v_superficial = rho * v_interstitial_neigh *
162  // v_superficial
163  const auto elem_coeff = vdotn * rho_face / eps_elem_face;
164  const auto neighbor_coeff = vdotn * rho_face / eps_neighbor_face;
165  _ae = elem_coeff * d_var_elem_face_d_elem_dof;
166  _elem_residual = elem_coeff * var_elem_face;
167  _an = -neighbor_coeff * d_var_neighbor_face_d_neighbor_dof;
168  _neighbor_residual = -neighbor_coeff * var_neighbor_face;
169  if (_approximate_as)
170  {
171  _ae = _cs / fi.elem().n_sides() * rho_face / eps_elem_face;
172  _an = _cs / fi.neighborPtr()->n_sides() * rho_face / eps_elem_face;
173  }
174  }
175  else
176  {
177  const auto [interp_coeffs, advected] =
178  interpCoeffsAndAdvected(*_rho_u, advected_face_arg, state);
179 
180  const auto elem_arg = elemArg();
181  const auto neighbor_arg = neighborArg();
182 
183  const auto rho_elem = _rho(elem_arg, state), rho_neighbor = _rho(neighbor_arg, state);
184  const auto eps_elem = epsilon()(elem_arg, state),
185  eps_neighbor = epsilon()(neighbor_arg, state);
186  const auto var_elem = advected.first / rho_elem * eps_elem,
187  var_neighbor = advected.second / rho_neighbor * eps_neighbor;
188 
189  _ae = vdotn * rho_elem / eps_elem * interp_coeffs.first;
190  // Minus sign because we apply a minus sign to the residual in computeResidual
191  _an = -vdotn * rho_neighbor / eps_neighbor * interp_coeffs.second;
192 
193  _elem_residual = _ae * var_elem - _an * var_neighbor;
195 
196  if (_approximate_as)
197  {
198  _ae = _cs / fi.elem().n_sides() * rho_elem / eps_elem;
199  _an = _cs / fi.neighborPtr()->n_sides() * rho_neighbor / eps_neighbor;
200  }
201  }
202  }
203 
204  _ae *= fi.faceArea() * fi.faceCoord();
205  _an *= fi.faceArea() * fi.faceCoord();
206  _elem_residual *= fi.faceArea() * fi.faceCoord();
207  _neighbor_residual *= fi.faceArea() * fi.faceCoord();
208 }
ADReal _neighbor_residual
The neighbor residual.
Moose::FV::InterpMethod _advected_interp_method
The interpolation method to use for the advected quantity.
Moose::ElemArg elemArg(bool correct_skewness=false) const
unsigned int number() const
std::unique_ptr< PiecewiseByBlockLambdaFunctor< ADReal > > _rho_u
Our local momentum functor.
const Elem & elem() const
const FaceInfo * _face_info
Moose::StateArg determineState() const
std::tuple< bool, ADReal, ADReal > isPorosityJumpFace(const Moose::Functor< ADReal > &porosity, const FaceInfo &fi, const Moose::StateArg &time)
Checks to see whether the porosity value jumps from one side to the other of the provided face...
Definition: NSFVUtils.C:58
const unsigned int _index
index x|y|z
auto raw_value(const Eigen::Map< T > &in)
Real & faceCoord()
RealVectorValue _normal
Real faceArea() const
ADReal _ae
The a coefficient for the element.
SystemBase & _sys
bool onBoundary(const FaceInfo &fi) const
LimiterType limiterType(InterpMethod interp_method)
const Elem * neighborPtr() const
const Real _cs
Characteristic speed.
const bool _approximate_as
Whether to approximately calculate the &#39;a&#39; coefficients.
bool skipForBoundary(const FaceInfo &fi) const override
const Moose::Functor< ADReal > & _rho
Density.
virtual const Moose::FunctorBase< ADReal > & epsilon() const
A virtual method that allows us to reuse all the code from free-flow for porous.
SubProblem & _subproblem
virtual bool isTransient() const=0
const Elem & neighbor() const
const Point & normal() const
unsigned int number() const
FaceInfo::VarFaceNeighbors _face_type
const Elem * elemPtr() const
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Moose::ElemArg neighborArg(bool correct_skewness=false) const
ADReal _an
The a coefficient for the neighbor.
ADReal _elem_residual
The element residual.
Moose::FaceArg singleSidedFaceArg(const FaceInfo *fi=nullptr, Moose::FV::LimiterType limiter_type=Moose::FV::LimiterType::CentralDifference, bool correct_skewness=false, const Moose::StateArg *state_limiter=nullptr) const
ADRealVectorValue velocity() const
MooseVariableFV< Real > & _var
std::pair< std::pair< T, T >, std::pair< T, T > > interpCoeffsAndAdvected(const FunctorBase< T > &functor, const FaceArg &face, const StateArg &time)
SolutionIterationType
VarFaceNeighbors faceType(const std::pair< unsigned int, unsigned int > &var_sys) const

◆ epsilon()

virtual const Moose::FunctorBase<ADReal>& INSFVMomentumAdvection::epsilon ( ) const
inlineprotectedvirtual

A virtual method that allows us to reuse all the code from free-flow for porous.

Reimplemented in PINSFVMomentumAdvection.

Definition at line 52 of file INSFVMomentumAdvection.h.

Referenced by WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), computeResidualsAndAData(), and initialSetup().

52 { return _unity_functor; }
const Moose::ConstantFunctor< ADReal > _unity_functor
A unity functor used in the epsilon virtual method.

◆ gatherRCData() [1/2]

void INSFVMomentumAdvection::gatherRCData ( const Elem &  elem)
inlinefinaloverridevirtual

Should be a non-empty implementation if the residual object is a FVElementalKernel and introduces residuals that are a function of the velocity, e.g.

friction and time-derivative terms.

Implements INSFVMomentumResidualObject.

Definition at line 37 of file INSFVMomentumAdvection.h.

37 {}

◆ gatherRCData() [2/2]

void INSFVMomentumAdvection::gatherRCData ( const FaceInfo fi)
finaloverridevirtual

Should be a non-empty implementation if the residual object is a FVFluxKernel and introduces residuals that are a function of the velocity, e.g.

advection, viscosity/diffusion, symmetry boundary conditions, etc.

Implements INSFVMomentumResidualObject.

Definition at line 277 of file INSFVMomentumAdvection.C.

278 {
279  if (skipForBoundary(fi))
280  return;
281 
282  const auto saved_velocity_interp_method = _velocity_interp_method;
284  // Fill-in the coefficients _ae and _an
286  _velocity_interp_method = saved_velocity_interp_method;
287 
290  _rc_uo.addToA(&fi.elem(), _index, _ae);
294 }
Moose::FV::InterpMethod _velocity_interp_method
The interpolation method to use for the velocity.
const Elem & elem() const
const unsigned int _index
index x|y|z
ADReal _ae
The a coefficient for the element.
RhieChowInterpolatorBase & _rc_uo
The Rhie Chow user object that is responsible for generating face velocities for advection terms...
const Elem * neighborPtr() const
bool skipForBoundary(const FaceInfo &fi) const override
FaceInfo::VarFaceNeighbors _face_type
virtual void addToA(const libMesh::Elem *elem, unsigned int component, const ADReal &value)=0
API for momentum residual objects that have on-diagonals for velocity call.
ADReal _an
The a coefficient for the neighbor.
virtual void computeResidualsAndAData(const FaceInfo &fi)
Helper method that computes the &#39;a&#39; coefficients and AD residuals.

◆ hasMaterialTimeDerivative()

virtual bool INSFVMomentumAdvection::hasMaterialTimeDerivative ( ) const
inlineoverrideprotectedvirtual
Returns
whether the equation this object adds to has a material time derivative. An instance when this will be false is the mass equation for incompressible flows

Implements INSFVAdvectionKernel.

Definition at line 47 of file INSFVMomentumAdvection.h.

47 { return true; }

◆ initialSetup() [1/2]

template<typename T >
void INSFVBCInterface::initialSetup ( T &  insfv_fk)
protectedinherited

setup all the boundary condition member information

Definition at line 79 of file INSFVBCInterface.h.

Referenced by INSFVAdvectionKernel::initialSetup().

80 {
81  const auto & mesh = insfv_fk.subProblem().mesh();
82 
83  std::set<BoundaryID> all_connected_boundaries;
84  const auto & blk_ids = insfv_fk.blockRestricted() ? insfv_fk.blockIDs() : mesh.meshSubdomains();
85  for (const auto blk_id : blk_ids)
86  {
87  const auto & connected_boundaries = mesh.getSubdomainInterfaceBoundaryIds(blk_id);
88  for (const auto bnd_id : connected_boundaries)
89  all_connected_boundaries.insert(bnd_id);
90  }
91 
92  for (const auto bnd_id : all_connected_boundaries)
93  {
94  setupFlowBoundaries(insfv_fk, bnd_id);
95  setupBoundaries<INSFVNoSlipWallBC>(
97  setupBoundaries<INSFVSlipWallBC>(
99  setupBoundaries<INSFVFreeSurfaceBC>(
101  setupBoundaries<INSFVSymmetryBC>(
103  }
104 }
std::set< BoundaryID > _slip_wall_boundaries
Boundary IDs with slip walls.
MeshBase & mesh
void setupFlowBoundaries(T &insfv_fk, BoundaryID bnd_id)
Query for INSFVBCs::INSFVFlowBC on bc_id and add if query successful.
std::set< BoundaryID > _symmetry_boundaries
Symmetry Boundary IDs.
std::set< BoundaryID > _no_slip_wall_boundaries
Boundary IDs with no slip walls.
std::set< BoundaryID > _free_surface_boundaries
Boundary IDs for free surfaces.

◆ initialSetup() [2/2]

void INSFVMomentumAdvection::initialSetup ( )
overridevirtual

Reimplemented from INSFVAdvectionKernel.

Definition at line 57 of file INSFVMomentumAdvection.C.

58 {
60 
61  _rho_u = std::make_unique<PiecewiseByBlockLambdaFunctor<ADReal>>(
62  "rho_u",
63  [this](const auto & r, const auto & t) -> ADReal
64  { return _rho(r, t) * _var(r, t) / epsilon()(r, t); },
65  std::set<ExecFlagType>({EXEC_ALWAYS}),
66  _mesh,
67  this->blockIDs());
68 }
MooseMesh & _mesh
std::unique_ptr< PiecewiseByBlockLambdaFunctor< ADReal > > _rho_u
Our local momentum functor.
virtual const std::set< SubdomainID > & blockIDs() const
const ExecFlagType EXEC_ALWAYS
const Moose::Functor< ADReal > & _rho
Density.
virtual const Moose::FunctorBase< ADReal > & epsilon() const
A virtual method that allows us to reuse all the code from free-flow for porous.
void initialSetup() override
MooseVariableFV< Real > & _var

◆ listOfCommonParams()

std::vector< std::string > INSFVMomentumAdvection::listOfCommonParams ( )
static
Returns
A list of the parameters that are common between this object and the NSFV action

Definition at line 31 of file INSFVMomentumAdvection.C.

Referenced by WCNSFVFlowPhysics::addMomentumAdvectionKernels().

32 {
33  return {"characteristic_speed"};
34 }

◆ skipForBoundary()

bool INSFVAdvectionKernel::skipForBoundary ( const FaceInfo fi) const
overrideprotectedvirtualinherited

Reimplemented from FVFluxKernel.

Definition at line 98 of file INSFVAdvectionKernel.C.

Referenced by INSFVTurbulentAdvection::computeJacobian(), computeJacobian(), INSFVTurbulentAdvection::computeResidual(), computeResidual(), WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), computeResidualsAndAData(), and gatherRCData().

99 {
100  // Boundaries to avoid come first, since they are always obeyed
101  if (avoidBoundary(fi))
102  return true;
103 
104  // We get this to check if we are on a kernel boundary or not
105  const bool on_boundary = onBoundary(fi);
106 
107  // We are either on a kernel boundary or on an internal sideset
108  // which is handled as a boundary
109  if (on_boundary || !fi.boundaryIDs().empty())
110  {
111  // Selected boundaries to force
112  for (const auto bnd_to_force : _boundaries_to_force)
113  if (fi.boundaryIDs().count(bnd_to_force))
114  return false;
115 
116  // If we have flux bcs then we do skip
117  const auto & [have_flux_bcs, flux_bcs] = _var.getFluxBCs(fi);
118  libmesh_ignore(have_flux_bcs);
119  for (const auto * const flux_bc : flux_bcs)
120  // If we have something like an average-value pressure constraint on a flow boundary, then we
121  // still want to execute this advection kernel on the boundary to ensure we're enforcing local
122  // conservation (mass in this example)
123  if (!dynamic_cast<const FVBoundaryScalarLagrangeMultiplierConstraint *>(flux_bc))
124  return true;
125 
126  // If we have a flow boundary without a replacement flux BC, then we must not skip. Mass and
127  // momentum are transported via advection across boundaries
128  for (const auto bc_id : fi.boundaryIDs())
129  if (_flow_boundaries.find(bc_id) != _flow_boundaries.end())
130  return false;
131  }
132 
133  // If not a flow boundary, then there should be no advection/flow in the normal direction, e.g. we
134  // should not contribute any advective flux. If we are on an internal face though, we still
135  // execute.
136  return on_boundary;
137 }
const std::set< BoundaryID > & boundaryIDs() const
std::unordered_set< BoundaryID > _boundaries_to_force
bool avoidBoundary(const FaceInfo &fi) const
std::pair< bool, std::vector< const FVFluxBC * > > getFluxBCs(const FaceInfo &fi) const
bool onBoundary(const FaceInfo &fi) const
void libmesh_ignore(const Args &...)
std::set< BoundaryID > _flow_boundaries
Flow Boundary IDs.
MooseVariableFV< Real > & _var

◆ uniqueParams()

InputParameters INSFVMomentumAdvection::uniqueParams ( )
static

Parameters of this object that should be added to the NSFV action that are unique to this object.

Definition at line 20 of file INSFVMomentumAdvection.C.

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

21 {
22  auto params = emptyInputParameters();
23  params.addParam<Real>(
24  "characteristic_speed",
25  "The characteristic speed. For porous medium simulations, this characteristic speed should "
26  "correspond to the superficial velocity, not the interstitial.");
27  return params;
28 }
InputParameters emptyInputParameters()
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ validParams()

InputParameters INSFVMomentumAdvection::validParams ( )
static

Definition at line 37 of file INSFVMomentumAdvection.C.

Referenced by PINSFVMomentumAdvection::validParams(), WCNSFVFlowPhysics::validParams(), and WCNSFV2PMomentumAdvectionSlip::validParams().

38 {
42  params.addRequiredParam<MooseFunctorName>(NS::density, "Density functor");
43  params.addClassDescription("Object for advecting momentum, e.g. rho*u");
44  return params;
45 }
static InputParameters uniqueParams()
Parameters of this object that should be added to the NSFV action that are unique to this object...
static const std::string density
Definition: NS.h:33
void addRequiredParam(const std::string &name, const std::string &doc_string)
void addClassDescription(const std::string &doc_string)
static InputParameters validParams()
static InputParameters validParams()

◆ velocity()

ADRealVectorValue INSFVAdvectionKernel::velocity ( ) const
inlineprotectedinherited
Returns
the advecting velocity

Definition at line 51 of file INSFVAdvectionKernel.h.

Referenced by INSFVMassAdvection::computeQpResidual(), INSFVScalarFieldAdvection::computeQpResidual(), INSFVEnergyAdvection::computeQpResidual(), and computeResidualsAndAData().

52 {
55 }
Moose::FV::InterpMethod _velocity_interp_method
The interpolation method to use for the velocity.
const FaceInfo * _face_info
Moose::StateArg determineState() const
THREAD_ID _tid
const RhieChowInterpolatorBase & _rc_vel_provider
The Rhie-Chow user object that provides us with the velocity.
virtual VectorValue< ADReal > getVelocity(const Moose::FV::InterpMethod m, const FaceInfo &fi, const Moose::StateArg &time, const THREAD_ID tid, bool subtract_mesh_velocity) const =0
Retrieve a face velocity.
virtual bool hasMaterialTimeDerivative() const =0

Member Data Documentation

◆ _advected_interp_method

Moose::FV::InterpMethod INSFVAdvectionKernel::_advected_interp_method
protectedinherited

◆ _ae

ADReal INSFVMomentumAdvection::_ae = 0
protected

The a coefficient for the element.

Definition at line 72 of file INSFVMomentumAdvection.h.

Referenced by WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), computeResidualsAndAData(), and gatherRCData().

◆ _all_boundaries

std::set<BoundaryID> INSFVBCInterface::_all_boundaries
protectedinherited

All the BoundaryIDs covered by our different types of INSFVBCs.

Definition at line 58 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::setupBoundaries(), and INSFVBCInterface::setupFlowBoundaries().

◆ _an

ADReal INSFVMomentumAdvection::_an = 0
protected

The a coefficient for the neighbor.

Definition at line 75 of file INSFVMomentumAdvection.h.

Referenced by WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), computeResidualsAndAData(), and gatherRCData().

◆ _approximate_as

const bool INSFVMomentumAdvection::_approximate_as
protected

Whether to approximately calculate the 'a' coefficients.

Definition at line 63 of file INSFVMomentumAdvection.h.

Referenced by WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), and computeResidualsAndAData().

◆ _cs

const Real INSFVMomentumAdvection::_cs
protected

Characteristic speed.

Definition at line 66 of file INSFVMomentumAdvection.h.

Referenced by WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), and computeResidualsAndAData().

◆ _elem_residual

ADReal INSFVMomentumAdvection::_elem_residual = 0
protected

◆ _flow_boundaries

std::set<BoundaryID> INSFVBCInterface::_flow_boundaries
protectedinherited

Flow Boundary IDs.

Definition at line 49 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::setupFlowBoundaries(), and INSFVAdvectionKernel::skipForBoundary().

◆ _free_surface_boundaries

std::set<BoundaryID> INSFVBCInterface::_free_surface_boundaries
protectedinherited

Boundary IDs for free surfaces.

Definition at line 46 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::initialSetup().

◆ _fully_developed_flow_boundaries

std::set<BoundaryID> INSFVBCInterface::_fully_developed_flow_boundaries
protectedinherited

Fully Developed Flow Boundary IDs. This is a subset of _flow_boundaries.

Definition at line 52 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::setupFlowBoundaries().

◆ _index

const unsigned int INSFVMomentumResidualObject::_index
protectedinherited

◆ _neighbor_residual

ADReal INSFVMomentumAdvection::_neighbor_residual = 0
protected

◆ _no_slip_wall_boundaries

std::set<BoundaryID> INSFVBCInterface::_no_slip_wall_boundaries
protectedinherited

Boundary IDs with no slip walls.

Definition at line 40 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::initialSetup().

◆ _rc_uo

RhieChowInterpolatorBase& INSFVMomentumResidualObject::_rc_uo
protectedinherited

◆ _rc_vel_provider

const RhieChowInterpolatorBase& INSFVAdvectionKernel::_rc_vel_provider
protectedinherited

The Rhie-Chow user object that provides us with the velocity.

Definition at line 47 of file INSFVAdvectionKernel.h.

Referenced by INSFVTurbulentAdvection::computeQpResidual(), and INSFVAdvectionKernel::velocity().

◆ _rho

const Moose::Functor<ADReal>& INSFVMomentumAdvection::_rho
protected

◆ _rho_u

std::unique_ptr<PiecewiseByBlockLambdaFunctor<ADReal> > INSFVMomentumAdvection::_rho_u
protected

Our local momentum functor.

Definition at line 69 of file INSFVMomentumAdvection.h.

Referenced by WCNSFV2PMomentumAdvectionSlip::computeResidualsAndAData(), computeResidualsAndAData(), and initialSetup().

◆ _slip_wall_boundaries

std::set<BoundaryID> INSFVBCInterface::_slip_wall_boundaries
protectedinherited

Boundary IDs with slip walls.

Definition at line 43 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::initialSetup().

◆ _symmetry_boundaries

std::set<BoundaryID> INSFVBCInterface::_symmetry_boundaries
protectedinherited

Symmetry Boundary IDs.

Definition at line 55 of file INSFVBCInterface.h.

Referenced by INSFVBCInterface::initialSetup().

◆ _unity_functor

const Moose::ConstantFunctor<ADReal> INSFVMomentumAdvection::_unity_functor {1}
protected

A unity functor used in the epsilon virtual method.

Definition at line 84 of file INSFVMomentumAdvection.h.

Referenced by epsilon().

◆ _velocity_interp_method

Moose::FV::InterpMethod INSFVAdvectionKernel::_velocity_interp_method
protectedinherited

The interpolation method to use for the velocity.

Definition at line 44 of file INSFVAdvectionKernel.h.

Referenced by INSFVTurbulentAdvection::computeQpResidual(), gatherRCData(), INSFVAdvectionKernel::INSFVAdvectionKernel(), and INSFVAdvectionKernel::velocity().


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