https://mooseframework.inl.gov
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Static Protected Member Functions | Protected Attributes | Static Protected Attributes | List of all members
PCNSFVKT Class Reference

Implements the centered Kurganov-Tadmor discretization of advective fluxes. More...

#include <PCNSFVKT.h>

Inheritance diagram for PCNSFVKT:
[legend]

Public Types

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

Public Member Functions

 PCNSFVKT (const InputParameters &params)
 
void computeResidual () override
 
void computeResidual (const FaceInfo &fi) override
 
void computeJacobian () override
 
void computeJacobian (const FaceInfo &fi) 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 initialSetup ()
 
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 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
 
std::pair< ADReal, ADRealcomputeAlphaAndOmega (const ADReal &u_elem_normal, const ADReal &u_neighbor_normal, const ADReal &c_elem, const ADReal &c_neighbor) const
 
std::string deduceFunctorName (const std::string &name) const
 
virtual ADReal gradUDotNormal (const Moose::StateArg &time, const bool correct_skewness) const
 
virtual bool skipForBoundary (const FaceInfo &fi) 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)
 

Static Protected Member Functions

static ADReal computeFaceFlux (const ADReal &alpha, const ADReal &omega, const ADReal &sup_vel_elem_normal, const ADReal &sup_vel_neighbor_normal, const ADReal &adv_quant_elem, const ADReal &adv_quant_neighbor)
 

Protected Attributes

const SinglePhaseFluidProperties_fluid
 fluid properties More...
 
const MooseEnum _eqn
 The equation we are solving, e.g. mass, momentum, fluid energy, or passive scalar. More...
 
const unsigned int _index
 When solving the momentum equation, the momentum component we are solving for. More...
 
std::unique_ptr< Moose::FV::Limiter< ADReal > > _limiter
 The slope limiter we will apply when interpolating from cell centroids to faces. More...
 
const bool _knp_for_omega
 Whether to use the Kurganov, Noelle, and Petrova method to compute the omega parameter for stabilization. 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
 
const ADMaterialProperty< Real > & _sup_vel_x_elem
 superficial velocities left == elem, right == neighbor More...
 
const ADMaterialProperty< Real > & _sup_vel_x_neighbor
 
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_x_elem
 
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_x_neighbor
 
const ADMaterialProperty< Real > & _sup_vel_y_elem
 
const ADMaterialProperty< Real > & _sup_vel_y_neighbor
 
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_y_elem
 
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_y_neighbor
 
const ADMaterialProperty< Real > & _sup_vel_z_elem
 
const ADMaterialProperty< Real > & _sup_vel_z_neighbor
 
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_z_elem
 
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_z_neighbor
 
const ADMaterialProperty< Real > & _T_fluid_elem
 fluid temperature left == elem, right == neighbor More...
 
const ADMaterialProperty< Real > & _T_fluid_neighbor
 
const ADMaterialProperty< RealVectorValue > & _grad_T_fluid_elem
 
const ADMaterialProperty< RealVectorValue > & _grad_T_fluid_neighbor
 
const ADMaterialProperty< Real > & _pressure_elem
 pressure left == elem, right == neighbor More...
 
const ADMaterialProperty< Real > & _pressure_neighbor
 
const ADMaterialProperty< RealVectorValue > & _grad_pressure_elem
 
const ADMaterialProperty< RealVectorValue > & _grad_pressure_neighbor
 
const MaterialProperty< Real > & _eps_elem
 porosity left == elem, right == neighbor More...
 
const MaterialProperty< Real > & _eps_neighbor
 
const ADVariableValue_scalar_elem
 passive scalar values left == elem, right == neighbor More...
 
const ADVariableValue_scalar_neighbor
 
const ADVariableGradient *const _grad_scalar_elem
 
const ADVariableGradient *const _grad_scalar_neighbor
 

Static Protected Attributes

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

Detailed Description

Implements the centered Kurganov-Tadmor discretization of advective fluxes.

Definition at line 29 of file PCNSFVKT.h.

Constructor & Destructor Documentation

◆ PCNSFVKT()

PCNSFVKT::PCNSFVKT ( const InputParameters params)

Definition at line 42 of file PCNSFVKT.C.

43  : FVFluxKernel(params),
44  _fluid(getUserObject<SinglePhaseFluidProperties>(NS::fluid)),
45 
46  _sup_vel_x_elem(getADMaterialProperty<Real>(NS::superficial_velocity_x)),
47  _sup_vel_x_neighbor(getNeighborADMaterialProperty<Real>(NS::superficial_velocity_x)),
49  getADMaterialProperty<RealVectorValue>(NS::grad(NS::superficial_velocity_x))),
51  getNeighborADMaterialProperty<RealVectorValue>(NS::grad(NS::superficial_velocity_x))),
52  _sup_vel_y_elem(getADMaterialProperty<Real>(NS::superficial_velocity_y)),
53  _sup_vel_y_neighbor(getNeighborADMaterialProperty<Real>(NS::superficial_velocity_y)),
55  getADMaterialProperty<RealVectorValue>(NS::grad(NS::superficial_velocity_y))),
57  getNeighborADMaterialProperty<RealVectorValue>(NS::grad(NS::superficial_velocity_y))),
58  _sup_vel_z_elem(getADMaterialProperty<Real>(NS::superficial_velocity_z)),
59  _sup_vel_z_neighbor(getNeighborADMaterialProperty<Real>(NS::superficial_velocity_z)),
61  getADMaterialProperty<RealVectorValue>(NS::grad(NS::superficial_velocity_z))),
63  getNeighborADMaterialProperty<RealVectorValue>(NS::grad(NS::superficial_velocity_z))),
64 
65  _T_fluid_elem(getADMaterialProperty<Real>(NS::T_fluid)),
66  _T_fluid_neighbor(getNeighborADMaterialProperty<Real>(NS::T_fluid)),
67  _grad_T_fluid_elem(getADMaterialProperty<RealVectorValue>(NS::grad(NS::T_fluid))),
68  _grad_T_fluid_neighbor(getNeighborADMaterialProperty<RealVectorValue>(NS::grad(NS::T_fluid))),
69  _pressure_elem(getADMaterialProperty<Real>(NS::pressure)),
70  _pressure_neighbor(getNeighborADMaterialProperty<Real>(NS::pressure)),
71  _grad_pressure_elem(getADMaterialProperty<RealVectorValue>(NS::grad(NS::pressure))),
72  _grad_pressure_neighbor(getNeighborADMaterialProperty<RealVectorValue>(NS::grad(NS::pressure))),
73  _eps_elem(getMaterialProperty<Real>(NS::porosity)),
74  _eps_neighbor(getNeighborMaterialProperty<Real>(NS::porosity)),
75  _eqn(getParam<MooseEnum>("eqn")),
76  _index(getParam<MooseEnum>("momentum_component")),
79  _grad_scalar_elem((_eqn == "scalar") ? &_var.adGradSln() : nullptr),
80  _grad_scalar_neighbor((_eqn == "scalar") ? &_var.adGradSlnNeighbor() : nullptr),
81  _limiter(Limiter<ADReal>::build(LimiterType(int(getParam<MooseEnum>("limiter"))))),
82  _knp_for_omega(getParam<bool>("knp_for_omega"))
83 {
84  if ((_eqn == "momentum") && !isParamValid("momentum_component"))
85  paramError("eqn",
86  "If 'momentum' is specified for 'eqn', then you must provide a parameter "
87  "value for 'momentum_component'");
88 }
const SinglePhaseFluidProperties & _fluid
fluid properties
Definition: PCNSFVKT.h:49
const unsigned int _index
When solving the momentum equation, the momentum component we are solving for.
Definition: PCNSFVKT.h:89
const ADMaterialProperty< RealVectorValue > & _grad_T_fluid_elem
Definition: PCNSFVKT.h:69
const ADVariableGradient *const _grad_scalar_neighbor
Definition: PCNSFVKT.h:95
const MooseEnum _eqn
The equation we are solving, e.g. mass, momentum, fluid energy, or passive scalar.
Definition: PCNSFVKT.h:86
static const std::string fluid
Definition: NS.h:87
const ADMaterialProperty< Real > & _pressure_elem
pressure left == elem, right == neighbor
Definition: PCNSFVKT.h:74
const MaterialProperty< Real > & _eps_neighbor
Definition: PCNSFVKT.h:82
const ADMaterialProperty< Real > & _sup_vel_x_elem
superficial velocities left == elem, right == neighbor
Definition: PCNSFVKT.h:52
const ADTemplateVariableGradient< Real > & adGradSlnNeighbor() const override
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_z_elem
Definition: PCNSFVKT.h:62
const ADVariableValue & _u_elem
const ADMaterialProperty< Real > & _sup_vel_y_neighbor
Definition: PCNSFVKT.h:57
bool isParamValid(const std::string &name) const
FVFluxKernel(const InputParameters &params)
static const std::string porosity
Definition: NS.h:104
const bool _knp_for_omega
Whether to use the Kurganov, Noelle, and Petrova method to compute the omega parameter for stabilizat...
Definition: PCNSFVKT.h:103
const ADMaterialProperty< Real > & _sup_vel_y_elem
Definition: PCNSFVKT.h:56
const ADMaterialProperty< Real > & _sup_vel_z_neighbor
Definition: PCNSFVKT.h:61
static const std::string T_fluid
Definition: NS.h:106
const ADVariableGradient *const _grad_scalar_elem
Definition: PCNSFVKT.h:94
static const std::string superficial_velocity_y
Definition: NS.h:51
const ADMaterialProperty< Real > & _T_fluid_neighbor
Definition: PCNSFVKT.h:68
LimiterType
std::unique_ptr< Moose::FV::Limiter< ADReal > > _limiter
The slope limiter we will apply when interpolating from cell centroids to faces.
Definition: PCNSFVKT.h:99
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_y_neighbor
Definition: PCNSFVKT.h:59
const ADVariableValue & _u_neighbor
void paramError(const std::string &param, Args... args) const
const ADMaterialProperty< Real > & _pressure_neighbor
Definition: PCNSFVKT.h:75
std::string grad(const std::string &var)
Definition: NS.h:91
const MaterialProperty< Real > & _eps_elem
porosity left == elem, right == neighbor
Definition: PCNSFVKT.h:81
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_x_elem
Definition: PCNSFVKT.h:54
static const std::string pressure
Definition: NS.h:56
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_x_neighbor
Definition: PCNSFVKT.h:55
MooseVariableFV< Real > & _var
const ADMaterialProperty< RealVectorValue > & _grad_pressure_neighbor
Definition: PCNSFVKT.h:77
const ADMaterialProperty< RealVectorValue > & _grad_pressure_elem
Definition: PCNSFVKT.h:76
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_z_neighbor
Definition: PCNSFVKT.h:63
const ADVariableValue & _scalar_elem
passive scalar values left == elem, right == neighbor
Definition: PCNSFVKT.h:92
const ADMaterialProperty< RealVectorValue > & _grad_T_fluid_neighbor
Definition: PCNSFVKT.h:70
const ADMaterialProperty< Real > & _T_fluid_elem
fluid temperature left == elem, right == neighbor
Definition: PCNSFVKT.h:67
const ADMaterialProperty< Real > & _sup_vel_z_elem
Definition: PCNSFVKT.h:60
const ADMaterialProperty< Real > & _sup_vel_x_neighbor
Definition: PCNSFVKT.h:53
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_y_elem
Definition: PCNSFVKT.h:58
static const std::string superficial_velocity_z
Definition: NS.h:52
const ADVariableValue & _scalar_neighbor
Definition: PCNSFVKT.h:93
static const std::string superficial_velocity_x
Definition: NS.h:50
const ADTemplateVariableGradient< Real > & adGradSln() const override

Member Function Documentation

◆ computeAlphaAndOmega()

std::pair< ADReal, ADReal > PCNSFVKT::computeAlphaAndOmega ( const ADReal u_elem_normal,
const ADReal u_neighbor_normal,
const ADReal c_elem,
const ADReal c_neighbor 
) const
protected

Definition at line 91 of file PCNSFVKT.C.

Referenced by computeQpResidual().

95 {
96  // Equations from Greenshields sans multiplication by area (which will be done in
97  // computeResidual/Jacobian
98  const auto psi_elem =
99  std::max({c_elem + u_elem_normal, c_neighbor + u_neighbor_normal, ADReal(0)});
100  const auto psi_neighbor =
101  std::max({c_elem - u_elem_normal, c_neighbor - u_neighbor_normal, ADReal(0)});
102  auto alpha = _knp_for_omega ? psi_elem / (psi_elem + psi_neighbor) : ADReal(0.5);
103  auto omega = _knp_for_omega ? alpha * (1 - alpha) * (psi_elem + psi_neighbor)
104  : alpha * std::max(psi_elem, psi_neighbor);
105 
106  // Do this to avoid new nonzero mallocs
107  const auto dummy_quant = 0 * (c_elem + u_elem_normal + c_neighbor + u_neighbor_normal);
108 
109  alpha += dummy_quant;
110  omega += dummy_quant;
111  return std::make_pair(std::move(alpha), std::move(omega));
112 }
auto max(const L &left, const R &right)
const bool _knp_for_omega
Whether to use the Kurganov, Noelle, and Petrova method to compute the omega parameter for stabilizat...
Definition: PCNSFVKT.h:103
static const std::string alpha
Definition: NS.h:134

◆ computeFaceFlux()

ADReal PCNSFVKT::computeFaceFlux ( const ADReal alpha,
const ADReal omega,
const ADReal sup_vel_elem_normal,
const ADReal sup_vel_neighbor_normal,
const ADReal adv_quant_elem,
const ADReal adv_quant_neighbor 
)
staticprotected

Definition at line 115 of file PCNSFVKT.C.

Referenced by computeQpResidual().

121 {
122  return alpha * (sup_vel_elem_normal * adv_quant_elem) +
123  (1 - alpha) * sup_vel_neighbor_normal * adv_quant_neighbor -
124  omega * (adv_quant_neighbor - adv_quant_elem);
125 }
static const std::string alpha
Definition: NS.h:134

◆ computeQpResidual()

ADReal PCNSFVKT::computeQpResidual ( )
overrideprotectedvirtual

Implements FVFluxKernel.

Reimplemented in PCNSFVKTDC.

Definition at line 128 of file PCNSFVKT.C.

Referenced by PCNSFVKTDC::computeQpResidual().

129 {
130  // Perform primitive interpolations
131  const auto pressure_elem = interpolate(*_limiter,
135  *_face_info,
136  /*elem_is_up=*/true);
137  const auto pressure_neighbor = interpolate(*_limiter,
141  *_face_info,
142  /*elem_is_up=*/false);
143  const auto T_fluid_elem = interpolate(*_limiter,
147  *_face_info,
148  /*elem_is_up=*/true);
149  const auto T_fluid_neighbor = interpolate(*_limiter,
153  *_face_info,
154  /*elem_is_up=*/false);
155  const auto sup_vel_x_elem = interpolate(*_limiter,
159  *_face_info,
160  /*elem_is_up=*/true);
161  const auto sup_vel_x_neighbor = interpolate(*_limiter,
165  *_face_info,
166  /*elem_is_up=*/false);
167  const auto sup_vel_y_elem = interpolate(*_limiter,
171  *_face_info,
172  /*elem_is_up=*/true);
173  const auto sup_vel_y_neighbor = interpolate(*_limiter,
177  *_face_info,
178  /*elem_is_up=*/false);
179  const auto sup_vel_z_elem = interpolate(*_limiter,
183  *_face_info,
184  /*elem_is_up=*/true);
185  const auto sup_vel_z_neighbor = interpolate(*_limiter,
189  *_face_info,
190  /*elem_is_up=*/false);
191 
192  const auto sup_vel_elem = VectorValue<ADReal>(sup_vel_x_elem, sup_vel_y_elem, sup_vel_z_elem);
193  const auto u_elem = sup_vel_elem / _eps_elem[_qp];
194  const auto rho_elem = _fluid.rho_from_p_T(pressure_elem, T_fluid_elem);
195  const auto specific_volume_elem = 1. / rho_elem;
196  const auto e_elem = _fluid.e_from_T_v(T_fluid_elem, specific_volume_elem);
197  const auto sup_vel_neighbor =
198  VectorValue<ADReal>(sup_vel_x_neighbor, sup_vel_y_neighbor, sup_vel_z_neighbor);
199  const auto u_neighbor = sup_vel_neighbor / _eps_neighbor[_qp];
200  const auto rho_neighbor = _fluid.rho_from_p_T(pressure_neighbor, T_fluid_neighbor);
201  const auto specific_volume_neighbor = 1. / rho_neighbor;
202  const auto e_neighbor = _fluid.e_from_T_v(T_fluid_neighbor, specific_volume_neighbor);
203 
204  const auto c_elem = _fluid.c_from_v_e(specific_volume_elem, e_elem);
205  const auto c_neighbor = _fluid.c_from_v_e(specific_volume_neighbor, e_neighbor);
206 
207  const auto sup_vel_elem_normal = sup_vel_elem * _face_info->normal();
208  const auto sup_vel_neighbor_normal = sup_vel_neighbor * _face_info->normal();
209  const auto u_elem_normal = u_elem * _face_info->normal();
210  const auto u_neighbor_normal = u_neighbor * _face_info->normal();
211 
212  const auto pr = computeAlphaAndOmega(u_elem_normal, u_neighbor_normal, c_elem, c_neighbor);
213  const auto & alpha = pr.first;
214  const auto & omega = pr.second;
215 
216  if (_eqn == "mass")
217  return computeFaceFlux(
218  alpha, omega, sup_vel_elem_normal, sup_vel_neighbor_normal, rho_elem, rho_neighbor);
219  else if (_eqn == "momentum")
220  {
221  const auto rhou_elem = u_elem(_index) * rho_elem;
222  const auto rhou_neighbor = u_neighbor(_index) * rho_neighbor;
223  return computeFaceFlux(alpha,
224  omega,
225  sup_vel_elem_normal,
226  sup_vel_neighbor_normal,
227  rhou_elem,
228  rhou_neighbor) +
229  _face_info->normal()(_index) * (alpha * _eps_elem[_qp] * pressure_elem +
230  (1 - alpha) * _eps_neighbor[_qp] * pressure_neighbor);
231  }
232  else if (_eqn == "energy")
233  {
234  const auto ht_elem = e_elem + 0.5 * u_elem * u_elem + pressure_elem / rho_elem;
235  const auto ht_neighbor =
236  e_neighbor + 0.5 * u_neighbor * u_neighbor + pressure_neighbor / rho_neighbor;
237  const auto rho_ht_elem = rho_elem * ht_elem;
238  const auto rho_ht_neighbor = rho_neighbor * ht_neighbor;
239  return computeFaceFlux(
240  alpha, omega, sup_vel_elem_normal, sup_vel_neighbor_normal, rho_ht_elem, rho_ht_neighbor);
241  }
242  else if (_eqn == "scalar")
243  {
244  const auto scalar_elem = interpolate(*_limiter,
245  _scalar_elem[_qp],
247  &(*_grad_scalar_elem)[_qp],
248  *_face_info,
249  true);
250  const auto scalar_neighbor = interpolate(*_limiter,
252  _scalar_elem[_qp],
254  *_face_info,
255  false);
256  const auto rhos_elem = rho_elem * scalar_elem;
257  const auto rhos_neighbor = rho_neighbor * scalar_neighbor;
258  return computeFaceFlux(
259  alpha, omega, sup_vel_elem_normal, sup_vel_neighbor_normal, rhos_elem, rhos_neighbor);
260  }
261  else
262  mooseError("Unrecognized enum type ", _eqn);
263 }
const SinglePhaseFluidProperties & _fluid
fluid properties
Definition: PCNSFVKT.h:49
const unsigned int _index
When solving the momentum equation, the momentum component we are solving for.
Definition: PCNSFVKT.h:89
const ADMaterialProperty< RealVectorValue > & _grad_T_fluid_elem
Definition: PCNSFVKT.h:69
const ADVariableGradient *const _grad_scalar_neighbor
Definition: PCNSFVKT.h:95
const FaceInfo * _face_info
const MooseEnum _eqn
The equation we are solving, e.g. mass, momentum, fluid energy, or passive scalar.
Definition: PCNSFVKT.h:86
const ADMaterialProperty< Real > & _pressure_elem
pressure left == elem, right == neighbor
Definition: PCNSFVKT.h:74
const MaterialProperty< Real > & _eps_neighbor
Definition: PCNSFVKT.h:82
const ADMaterialProperty< Real > & _sup_vel_x_elem
superficial velocities left == elem, right == neighbor
Definition: PCNSFVKT.h:52
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_z_elem
Definition: PCNSFVKT.h:62
const ADMaterialProperty< Real > & _sup_vel_y_neighbor
Definition: PCNSFVKT.h:57
static ADReal computeFaceFlux(const ADReal &alpha, const ADReal &omega, const ADReal &sup_vel_elem_normal, const ADReal &sup_vel_neighbor_normal, const ADReal &adv_quant_elem, const ADReal &adv_quant_neighbor)
Definition: PCNSFVKT.C:115
const ADMaterialProperty< Real > & _sup_vel_y_elem
Definition: PCNSFVKT.h:56
const ADMaterialProperty< Real > & _sup_vel_z_neighbor
Definition: PCNSFVKT.h:61
const ADVariableGradient *const _grad_scalar_elem
Definition: PCNSFVKT.h:94
const unsigned int _qp
const ADMaterialProperty< Real > & _T_fluid_neighbor
Definition: PCNSFVKT.h:68
std::unique_ptr< Moose::FV::Limiter< ADReal > > _limiter
The slope limiter we will apply when interpolating from cell centroids to faces.
Definition: PCNSFVKT.h:99
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_y_neighbor
Definition: PCNSFVKT.h:59
const Point & normal() const
const ADMaterialProperty< Real > & _pressure_neighbor
Definition: PCNSFVKT.h:75
const MaterialProperty< Real > & _eps_elem
porosity left == elem, right == neighbor
Definition: PCNSFVKT.h:81
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_x_elem
Definition: PCNSFVKT.h:54
static const std::string alpha
Definition: NS.h:134
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_x_neighbor
Definition: PCNSFVKT.h:55
void mooseError(Args &&... args) const
const ADMaterialProperty< RealVectorValue > & _grad_pressure_neighbor
Definition: PCNSFVKT.h:77
const ADMaterialProperty< RealVectorValue > & _grad_pressure_elem
Definition: PCNSFVKT.h:76
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_z_neighbor
Definition: PCNSFVKT.h:63
const ADVariableValue & _scalar_elem
passive scalar values left == elem, right == neighbor
Definition: PCNSFVKT.h:92
const ADMaterialProperty< RealVectorValue > & _grad_T_fluid_neighbor
Definition: PCNSFVKT.h:70
const ADMaterialProperty< Real > & _T_fluid_elem
fluid temperature left == elem, right == neighbor
Definition: PCNSFVKT.h:67
void interpolate(InterpMethod m, T &result, const T2 &value1, const T3 &value2, const FaceInfo &fi, const bool one_is_elem)
const ADMaterialProperty< Real > & _sup_vel_z_elem
Definition: PCNSFVKT.h:60
const ADMaterialProperty< Real > & _sup_vel_x_neighbor
Definition: PCNSFVKT.h:53
const ADMaterialProperty< RealVectorValue > & _grad_sup_vel_y_elem
Definition: PCNSFVKT.h:58
const ADVariableValue & _scalar_neighbor
Definition: PCNSFVKT.h:93
std::pair< ADReal, ADReal > computeAlphaAndOmega(const ADReal &u_elem_normal, const ADReal &u_neighbor_normal, const ADReal &c_elem, const ADReal &c_neighbor) const
Definition: PCNSFVKT.C:91

◆ validParams()

InputParameters PCNSFVKT::validParams ( )
static

Definition at line 20 of file PCNSFVKT.C.

Referenced by PCNSFVKTDC::validParams().

21 {
23  params.addClassDescription("Computes the residual of advective term using finite volume method.");
24  params.addRequiredParam<UserObjectName>(NS::fluid, "Fluid userobject");
25  MooseEnum eqn("mass momentum energy scalar");
26  params.addRequiredParam<MooseEnum>("eqn", eqn, "The equation you're solving.");
27  MooseEnum momentum_component("x=0 y=1 z=2");
28  params.addParam<MooseEnum>("momentum_component",
29  momentum_component,
30  "The component of the momentum equation that this kernel applies to.");
31  params.addParam<MooseEnum>(
32  "limiter", moose_limiter_type, "The limiter to apply during interpolation.");
33  params.set<unsigned short>("ghost_layers") = 2;
34  params.addParam<bool>(
35  "knp_for_omega",
36  true,
37  "Whether to use the Kurganov, Noelle, and Petrova method to compute the omega parameter for "
38  "stabilization. If false, then the Kurganov-Tadmor method will be used.");
39  return params;
40 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
T & set(const std::string &name, bool quiet_mode=false)
static const std::string fluid
Definition: NS.h:87
void addRequiredParam(const std::string &name, const std::string &doc_string)
static InputParameters validParams()
const MooseEnum moose_limiter_type
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _eps_elem

const MaterialProperty<Real>& PCNSFVKT::_eps_elem
protected

porosity left == elem, right == neighbor

Definition at line 81 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _eps_neighbor

const MaterialProperty<Real>& PCNSFVKT::_eps_neighbor
protected

Definition at line 82 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _eqn

const MooseEnum PCNSFVKT::_eqn
protected

The equation we are solving, e.g. mass, momentum, fluid energy, or passive scalar.

Definition at line 86 of file PCNSFVKT.h.

Referenced by computeQpResidual(), and PCNSFVKT().

◆ _fluid

const SinglePhaseFluidProperties& PCNSFVKT::_fluid
protected

fluid properties

Definition at line 49 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_pressure_elem

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_pressure_elem
protected

Definition at line 76 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_pressure_neighbor

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_pressure_neighbor
protected

Definition at line 77 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_scalar_elem

const ADVariableGradient* const PCNSFVKT::_grad_scalar_elem
protected

Definition at line 94 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_scalar_neighbor

const ADVariableGradient* const PCNSFVKT::_grad_scalar_neighbor
protected

Definition at line 95 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_sup_vel_x_elem

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_sup_vel_x_elem
protected

Definition at line 54 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_sup_vel_x_neighbor

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_sup_vel_x_neighbor
protected

Definition at line 55 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_sup_vel_y_elem

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_sup_vel_y_elem
protected

Definition at line 58 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_sup_vel_y_neighbor

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_sup_vel_y_neighbor
protected

Definition at line 59 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_sup_vel_z_elem

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_sup_vel_z_elem
protected

Definition at line 62 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_sup_vel_z_neighbor

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_sup_vel_z_neighbor
protected

Definition at line 63 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_T_fluid_elem

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_T_fluid_elem
protected

Definition at line 69 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _grad_T_fluid_neighbor

const ADMaterialProperty<RealVectorValue>& PCNSFVKT::_grad_T_fluid_neighbor
protected

Definition at line 70 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _index

const unsigned int PCNSFVKT::_index
protected

When solving the momentum equation, the momentum component we are solving for.

Definition at line 89 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _knp_for_omega

const bool PCNSFVKT::_knp_for_omega
protected

Whether to use the Kurganov, Noelle, and Petrova method to compute the omega parameter for stabilization.

If false, then the Kurganov-Tadmor method will be used

Definition at line 103 of file PCNSFVKT.h.

Referenced by computeAlphaAndOmega().

◆ _limiter

std::unique_ptr<Moose::FV::Limiter<ADReal> > PCNSFVKT::_limiter
protected

The slope limiter we will apply when interpolating from cell centroids to faces.

Definition at line 99 of file PCNSFVKT.h.

Referenced by computeQpResidual(), and PCNSFVKTDC::computeQpResidual().

◆ _pressure_elem

const ADMaterialProperty<Real>& PCNSFVKT::_pressure_elem
protected

pressure left == elem, right == neighbor

Definition at line 74 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _pressure_neighbor

const ADMaterialProperty<Real>& PCNSFVKT::_pressure_neighbor
protected

Definition at line 75 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _scalar_elem

const ADVariableValue& PCNSFVKT::_scalar_elem
protected

passive scalar values left == elem, right == neighbor

Definition at line 92 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _scalar_neighbor

const ADVariableValue& PCNSFVKT::_scalar_neighbor
protected

Definition at line 93 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _sup_vel_x_elem

const ADMaterialProperty<Real>& PCNSFVKT::_sup_vel_x_elem
protected

superficial velocities left == elem, right == neighbor

Definition at line 52 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _sup_vel_x_neighbor

const ADMaterialProperty<Real>& PCNSFVKT::_sup_vel_x_neighbor
protected

Definition at line 53 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _sup_vel_y_elem

const ADMaterialProperty<Real>& PCNSFVKT::_sup_vel_y_elem
protected

Definition at line 56 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _sup_vel_y_neighbor

const ADMaterialProperty<Real>& PCNSFVKT::_sup_vel_y_neighbor
protected

Definition at line 57 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _sup_vel_z_elem

const ADMaterialProperty<Real>& PCNSFVKT::_sup_vel_z_elem
protected

Definition at line 60 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _sup_vel_z_neighbor

const ADMaterialProperty<Real>& PCNSFVKT::_sup_vel_z_neighbor
protected

Definition at line 61 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _T_fluid_elem

const ADMaterialProperty<Real>& PCNSFVKT::_T_fluid_elem
protected

fluid temperature left == elem, right == neighbor

Definition at line 67 of file PCNSFVKT.h.

Referenced by computeQpResidual().

◆ _T_fluid_neighbor

const ADMaterialProperty<Real>& PCNSFVKT::_T_fluid_neighbor
protected

Definition at line 68 of file PCNSFVKT.h.

Referenced by computeQpResidual().


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