This is the base class all materials should use if you are trying to use the Navier-Stokes Kernels. More...

#include <NavierStokesMaterial.h>

Inheritance diagram for NavierStokesMaterial:

Public Types
enum	ConstantTypeEnum { ConstantTypeEnum::NONE, ConstantTypeEnum::ELEMENT, ConstantTypeEnum::SUBDOMAIN }

enum	TEST_TYPE

typedef DataFileName	DataFileParameterType

Public Member Functions
	NavierStokesMaterial (const InputParameters &parameters)

virtual const dof_id_type &	getElementID (const std::string &id_parameter_name, unsigned int comp=0) const override

dof_id_type	getElementID (const Elem *elem, unsigned int elem_id_index) const

virtual const dof_id_type &	getElementIDNeighbor (const std::string &id_parameter_name, unsigned int comp=0) const override

virtual const dof_id_type &	getElementIDByName (const std::string &id_parameter_name) const override

virtual const dof_id_type &	getElementIDNeighborByName (const std::string &id_parameter_name) const override

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false, bool no_dep=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

virtual bool	isBoundaryMaterial () const override

virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const override

virtual void	subdomainSetup () override

bool	ghostable () const override final

virtual void	resolveOptionalProperties () override

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

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 GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

virtual void	initStatefulProperties (unsigned int n_points)

virtual bool	isInterfaceMaterial ()

virtual void	resetProperties ()

virtual void	computePropertiesAtQp (unsigned int qp)

const MaterialProperty< T > &	getZeroMaterialPropertyByName (Ts... args)

virtual const std::set< std::string > &	getRequestedItems () override

virtual const std::set< std::string > &	getSuppliedItems () override

const std::set< unsigned int > &	getSuppliedPropIDs ()

void	checkStatefulSanity () const

std::set< OutputName >	getOutputs ()

bool	hasStatefulProperties () const

void	setFaceInfo (const FaceInfo &fi)

void	setActiveProperties (const std::unordered_set< unsigned int > &needed_props)

bool	forceStatefulInit () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

MooseApp &	getMooseApp () 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 InputParameters &	parameters () 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

T	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

const std::vector< SubdomainName > &	blocks () const

unsigned int	numBlocks () const

virtual const std::set< SubdomainID > &	blockIDs () const

unsigned int	blocksMaxDimension () const

bool	hasBlocks (const SubdomainName &name) const

bool	hasBlocks (const std::vector< SubdomainName > &names) const

bool	hasBlocks (const std::set< SubdomainName > &names) const

bool	hasBlocks (SubdomainID id) const

bool	hasBlocks (const std::vector< SubdomainID > &ids) const

bool	hasBlocks (const std::set< SubdomainID > &ids) const

bool	isBlockSubset (const std::set< SubdomainID > &ids) const

bool	isBlockSubset (const std::vector< SubdomainID > &ids) const

bool	hasBlockMaterialProperty (const std::string &prop_name)

const std::set< SubdomainID > &	meshBlockIDs () const

virtual bool	blockRestricted () const

virtual void	checkVariable (const MooseVariableFieldBase &variable) const

virtual const std::set< BoundaryID > &	boundaryIDs () const

const std::vector< BoundaryName > &	boundaryNames () const

unsigned int	numBoundaryIDs () const

bool	hasBoundary (const BoundaryName &name) const

bool	hasBoundary (const std::vector< BoundaryName > &names) const

bool	hasBoundary (const BoundaryID &id) const

bool	hasBoundary (const std::vector< BoundaryID > &ids, TEST_TYPE type=ALL) const

bool	hasBoundary (const std::set< BoundaryID > &ids, TEST_TYPE type=ALL) const

bool	isBoundarySubset (const std::set< BoundaryID > &ids) const

bool	isBoundarySubset (const std::vector< BoundaryID > &ids) const

bool	hasBoundaryMaterialProperty (const std::string &prop_name) const

virtual bool	boundaryRestricted () const

const std::set< BoundaryID > &	meshBoundaryIDs () const

virtual bool	checkVariableBoundaryIntegrity () const

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () 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)

const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()

const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const

const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const

const Function &	getFunction (const std::string &name) const

const Function &	getFunctionByName (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 UserObject &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObject &	getUserObjectBaseByName (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 VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (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

virtual void	meshChanged ()

void	buildOutputHideVariableList (std::set< std::string > variable_names)

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 unsigned int	getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const

virtual unsigned int	getElementIDIndexByName (const std::string &id_name) const

bool	hasElementID (const std::string &id_name) const

dof_id_type	maxElementID (unsigned int elem_id_index) const

dof_id_type	minElementID (unsigned int elem_id_index) const

bool	areElemIDsIdentical (const std::string &id_name1, const std::string &id_name2) const

std::unordered_map< dof_id_type, std::set< dof_id_type > >	getElemIDMapping (const std::string &id_name1, const std::string &id_name2) const

std::set< dof_id_type >	getAllElemIDs (unsigned int elem_id_index) const

std::set< dof_id_type >	getElemIDsOnBlocks (unsigned int elem_id_index, const std::set< SubdomainID > &blks) 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 ADVariableValue *	getADDefaultValue (const std::string &var_name) const

const ADVectorVariableValue *	getADDefaultVectorValue (const std::string &var_name) const

const ADVariableGradient &	getADDefaultGradient () const

const ADVectorVariableGradient &	getADDefaultVectorGradient () const

const ADVariableSecond &	getADDefaultSecond () const

const ADVectorVariableCurl &	getADDefaultCurl () 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 std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const std::string &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 std::string &prop_name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const std::string &prop_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 std::string &prop_name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const std::string &prop_name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const std::string &prop_name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const std::string &prop_name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const std::string &prop_name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const std::string &prop_name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)

std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)

std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)

std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)

std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)

void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

void	statefulPropertiesAllowed (bool)

bool	getMaterialPropertyCalled () const

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

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)

MaterialProperty< T > &	declarePropertyByName (const std::string &prop_name)

MaterialProperty< T > &	declarePropertyByName (const std::string &prop_name)

MaterialProperty< T > &	declareProperty (const std::string &name)

MaterialProperty< T > &	declareProperty (const std::string &name)

ADMaterialProperty< T > &	declareADPropertyByName (const std::string &prop_name)

ADMaterialProperty< T > &	declareADPropertyByName (const std::string &prop_name)

ADMaterialProperty< T > &	declareADProperty (const std::string &name)

ADMaterialProperty< T > &	declareADProperty (const std::string &name)

auto &	declareGenericProperty (const std::string &prop_name)

auto &	declareGenericProperty (const std::string &prop_name)

GenericMaterialProperty< T, is_ad > &	declareGenericPropertyByName (const std::string &prop_name)

GenericMaterialProperty< T, is_ad > &	declareGenericPropertyByName (const std::string &prop_name)

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &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	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 PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (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

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)

PenetrationLocator &	getPenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)

PenetrationLocator &	getQuadraturePenetrationLocator (const BoundaryName &primary, const BoundaryName &secondary, Order order)

NearestNodeLocator &	getNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)

NearestNodeLocator &	getQuadratureNearestNodeLocator (const BoundaryName &primary, const BoundaryName &secondary)

bool	requiresGeometricSearch () const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static std::deque< MaterialBase *>	buildRequiredMaterials (const Consumers &mat_consumers, const std::vector< std::shared_ptr< MaterialBase >> &mats, const bool allow_stateful)

static bool	restricted (const std::set< BoundaryID > &ids)

static void	sort (typename std::vector< T > &vector)

static void	sortDFS (typename std::vector< T > &vector)

static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)

static std::string	deduceFunctorName (const std::string &name, const InputParameters &params)

Public Attributes
	ALL

	ANY

const ConsoleStream	_console

Static Public Attributes
static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

Protected Types
enum	QP_Data_Type

Protected Member Functions
virtual void	computeProperties ()
	Must be called after the child class computes dynamic_viscocity. More...

virtual void	checkMaterialProperty (const std::string &name, const unsigned int state) override

virtual const MaterialData &	materialData () const override

virtual MaterialData &	materialData () override

virtual const QBase &	qRule () const override

virtual void	computeQpProperties ()

virtual void	resetQpProperties ()

virtual void	initQpStatefulProperties ()

virtual const FEProblemBase &	miProblem () const

virtual FEProblemBase &	miProblem ()

bool	isPropertyActive (const unsigned int prop_id) const

void	registerPropName (const std::string &prop_name, bool is_get, const unsigned int state)

void	checkExecutionStage ()

void	checkExecutionStage ()

virtual bool	hasBlockMaterialPropertyHelper (const std::string &prop_name)

void	initializeBlockRestrictable (const MooseObject *moose_object)

Moose::CoordinateSystemType	getBlockCoordSystem ()

bool	hasBoundaryMaterialPropertyHelper (const std::string &prop_name) 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 VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarValue (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 VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const

const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) 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

std::string	deduceFunctorName (const std::string &name) const

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)

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

virtual void	coupledCallback (const std::string &, bool) const

virtual bool	isCoupled (const std::string &var_name, unsigned int i=0) const

virtual bool	isCoupledConstant (const std::string &var_name) const

unsigned int	coupledComponents (const std::string &var_name) const

VariableName	coupledName (const std::string &var_name, unsigned int comp=0) const

std::vector< VariableName >	coupledNames (const std::string &var_name) const

virtual unsigned int	coupled (const std::string &var_name, unsigned int comp=0) const

std::vector< unsigned int >	coupledIndices (const std::string &var_name) const

virtual const VariableValue &	coupledValue (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 VariableValue &	coupledValueLower (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVariableValue *>	adCoupledValues (const std::string &var_name) const

const ADVariableValue &	adCoupledLowerValue (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableValue &	adCoupledVectorValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVectorVariableValue *>	adCoupledVectorValues (const std::string &var_name) const

virtual const VariableValue &	coupledVectorTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const VariableValue &	coupledVectorTagValue (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 ArrayVariableValue &	coupledVectorTagArrayValue (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const ArrayVariableValue &	coupledVectorTagArrayValue (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 VariableGradient &	coupledVectorTagGradient (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const VariableGradient &	coupledVectorTagGradient (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 ArrayVariableGradient &	coupledVectorTagArrayGradient (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const ArrayVariableGradient &	coupledVectorTagArrayGradient (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 VariableValue &	coupledVectorTagDofValue (const std::string &var_name, TagID tag, unsigned int index=0) const

virtual const VariableValue &	coupledVectorTagDofValue (const std::string &var_names, const std::string &tag_name, unsigned int index=0) const

const ArrayVariableValue &	coupledVectorTagArrayDofValue (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 VariableValue &	coupledMatrixTagValue (const std::string &var_names, TagID tag, unsigned int index=0) const

virtual const VariableValue &	coupledMatrixTagValue (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 VectorVariableValue &	coupledVectorValue (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayValue (const std::string &var_name, unsigned int comp=0) const

std::vector< const ArrayVariableValue *>	coupledArrayValues (const std::string &var_name) const

MooseWritableVariable &	writableVariable (const std::string &var_name, unsigned int comp=0)

virtual VariableValue &	writableCoupledValue (const std::string &var_name, unsigned int comp=0)

void	checkWritableVar (MooseWritableVariable *var)

virtual const VariableValue &	coupledValueOld (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 VariableValue &	coupledValueOlder (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledValuesOlder (const std::string &var_name) const

virtual const VariableValue &	coupledValuePreviousNL (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorValueOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorValueOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayValueOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayValueOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradient (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableGradient *>	coupledGradients (const std::string &var_name) const

const ADVariableGradient &	adCoupledGradient (const std::string &var_name, unsigned int comp=0) const

const ADVariableGradient &	adCoupledGradientDot (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 ADVectorVariableGradient &	adCoupledVectorGradient (const std::string &var_name, unsigned int comp=0) const

const ADVariableSecond &	adCoupledSecond (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableSecond &	adCoupledVectorSecond (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientOld (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableGradient *>	coupledGradientsOld (const std::string &var_name) const

virtual const VariableGradient &	coupledGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientPreviousNL (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableGradient &	coupledGradientDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorGradient (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableGradient &	coupledVectorGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradient (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradientOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradientOlder (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableGradient &	coupledArrayGradientDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableCurl &	coupledCurl (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableCurl &	coupledCurlOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableCurl &	coupledCurlOlder (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableCurl &	adCoupledCurl (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableDivergence &	coupledDiv (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableDivergence &	coupledDivOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableDivergence &	coupledDivOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecond (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecondOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecondOlder (const std::string &var_name, unsigned int comp=0) const

virtual const VariableSecond &	coupledSecondPreviousNL (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDot (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledDots (const std::string &var_name) const

virtual const VariableValue &	coupledDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotDotOld (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledDot (const std::string &var_name, unsigned int comp=0) const

std::vector< const ADVariableValue *>	adCoupledDots (const std::string &var_name) const

const ADVariableValue &	adCoupledDotDot (const std::string &var_name, unsigned int comp=0) const

const ADVectorVariableValue &	adCoupledVectorDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VectorVariableValue &	coupledVectorDotDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledVectorDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledVectorDotDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDot (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const ArrayVariableValue &	coupledArrayDotDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotDu (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDotDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledArrayDotDu (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 >::type &	adCoupledNodalValue (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 VariableValue &	coupledNodalDotDot (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNodalDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledNodalDotDotOld (const std::string &var_name, unsigned int comp=0) const

virtual const VariableValue &	coupledDofValues (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledAllDofValues (const std::string &var_name) const

virtual const VariableValue &	coupledDofValuesOld (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledAllDofValuesOld (const std::string &var_name) const

virtual const VariableValue &	coupledDofValuesOlder (const std::string &var_name, unsigned int comp=0) const

std::vector< const VariableValue *>	coupledAllDofValuesOlder (const std::string &var_name) const

virtual const ArrayVariableValue &	coupledArrayDofValues (const std::string &var_name, unsigned int comp=0) const

virtual const ADVariableValue &	adCoupledDofValues (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adZeroValue () const

const ADVariableGradient &	adZeroGradient () const

const ADVariableSecond &	adZeroSecond () 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 MooseVariableFieldBase *	getFEVar (const std::string &var_name, unsigned int comp) const

const MooseVariableFieldBase *	getFieldVar (const std::string &var_name, unsigned int comp) const

MooseVariableFieldBase *	getFieldVar (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)

MooseVariable *	getVar (const std::string &var_name, unsigned int comp)

const MooseVariable *	getVar (const std::string &var_name, unsigned int comp) const

VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp)

const VectorMooseVariable *	getVectorVar (const std::string &var_name, unsigned int comp) const

ArrayMooseVariable *	getArrayVar (const std::string &var_name, unsigned int comp)

const ArrayMooseVariable *	getArrayVar (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

void	markMatPropRequested (const std::string &)

MaterialPropertyName	getMaterialPropertyName (const std::string &name) 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 unsigned int	_mesh_dimension

const VariableGradient &	_grad_u

const VariableGradient &	_grad_v

const VariableGradient &	_grad_w

MaterialProperty< RealTensorValue > &	_viscous_stress_tensor

MaterialProperty< Real > &	_thermal_conductivity

MaterialProperty< Real > &	_dynamic_viscosity

MaterialProperty< std::vector< RealTensorValue > > &	_calA

MaterialProperty< std::vector< RealTensorValue > > &	_calC

MaterialProperty< std::vector< std::vector< RealTensorValue > > > &	_calE

std::vector< const VariableGradient * >	_vel_grads

const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

const VariableValue &	_temperature

const VariableValue &	_specific_total_enthalpy

const VariableValue &	_rho

const VariableValue &	_rho_u

const VariableValue &	_rho_v

const VariableValue &	_rho_w

const VariableValue &	_rho_et

const VariableValue &	_drho_dt

const VariableValue &	_drhou_dt

const VariableValue &	_drhov_dt

const VariableValue &	_drhow_dt

const VariableValue &	_drhoE_dt

const VariableGradient &	_grad_rho

const VariableGradient &	_grad_rho_u

const VariableGradient &	_grad_rho_v

const VariableGradient &	_grad_rho_w

const VariableGradient &	_grad_rho_et

MaterialProperty< Real > &	_hsupg

MaterialProperty< Real > &	_tauc

MaterialProperty< Real > &	_taum

MaterialProperty< Real > &	_taue

MaterialProperty< std::vector< Real > > &	_strong_residuals

const IdealGasFluidProperties &	_fp

	CURR

	PREV

bool	_bnd

bool	_neighbor

const MooseArray< Point > &	_q_point

const QBase *const &	_qrule

const MooseArray< Real > &	_JxW

const Elem *const &	_current_elem

const SubdomainID &	_current_subdomain_id

const unsigned int &	_current_side

const ConstantTypeEnum	_constant_option

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

THREAD_ID	_tid

Assembly &	_assembly

unsigned int	_qp

const MooseArray< Real > &	_coord

const MooseArray< Point > &	_normals

MooseMesh &	_mesh

const Moose::CoordinateSystemType &	_coord_sys

std::set< std::string >	_requested_props

std::set< std::string >	_supplied_props

std::set< unsigned int >	_supplied_prop_ids

std::unordered_set< unsigned int >	_active_prop_ids

const bool	_compute

std::unordered_map< unsigned int, unsigned int >	_props_to_min_states

std::vector< unsigned int >	_displacements

bool	_has_stateful_property

bool	_overrides_init_stateful_props

const FaceInfo *	_face_info

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const MaterialData *	_blk_material_data

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

GeometricSearchData &	_geometric_search_data

bool	_requires_geometric_search

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

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

const Parallel::Communicator &	_communicator

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Private Member Functions
void	computeHSUPG (unsigned int qp)

void	computeTau (unsigned int qp)

void	computeStrongResiduals (unsigned int qp)

Detailed Description

This is the base class all materials should use if you are trying to use the Navier-Stokes Kernels.

Note that the derived class just needs to compute dynamic_viscocity then call this class's computeProperties() function.

Also make sure that the derived class's validParams function just adds to this class's validParams.

Finally, note that this Material isn't registered with the MaterialFactory. The reason is that by itself this material doesn't work! You must derive from this material and compute dynamic_viscocity!

Definition at line 32 of file NavierStokesMaterial.h.

Constructor & Destructor Documentation

◆ NavierStokesMaterial()

NavierStokesMaterial::NavierStokesMaterial ( const InputParameters & parameters )

Definition at line 50 of file NavierStokesMaterial.C.

   : Material(parameters),
     _mesh_dimension(_mesh.dimension()),
     _grad_u(coupledGradient(nms::velocity_x)),
     _grad_v(_mesh_dimension >= 2 ? coupledGradient(nms::velocity_y) : _grad_zero),
     _grad_w(_mesh_dimension == 3 ? coupledGradient(nms::velocity_z) : _grad_zero),
 
     _viscous_stress_tensor(declareProperty<RealTensorValue>("viscous_stress_tensor")),
     _thermal_conductivity(declareProperty<Real>("thermal_conductivity")),
 
     // Declared here but _not_ calculated here
     // (See e.g. derived class, bighorn/include/materials/FluidTC1.h)
     _dynamic_viscosity(declareProperty<Real>("dynamic_viscosity")),
 
     // The momentum components of the inviscid flux Jacobians.
     _calA(declareProperty<std::vector<RealTensorValue>>("calA")),
 
     // "Velocity column" matrices
     _calC(declareProperty<std::vector<RealTensorValue>>("calC")),
 
     // Energy equation inviscid flux matrices, "cal E_{kl}" in the notes.
     _calE(declareProperty<std::vector<std::vector<RealTensorValue>>>("calE")),
     _vel_grads({&_grad_u, &_grad_v, &_grad_w}),
 
     // Coupled solution values needed for computing SUPG stabilization terms
     _u_vel(coupledValue(nms::velocity_x)),
     _v_vel(_mesh.dimension() >= 2 ? coupledValue(nms::velocity_y) : _zero),
     _w_vel(_mesh.dimension() == 3 ? coupledValue(nms::velocity_z) : _zero),
 
     _temperature(coupledValue(nms::temperature)),
     _specific_total_enthalpy(coupledValue(nms::specific_total_enthalpy)),
 
     // Coupled solution values
     _rho(coupledValue(nms::density)),
     _rho_u(coupledValue(nms::momentum_x)),
     _rho_v(_mesh.dimension() >= 2 ? coupledValue(nms::momentum_y) : _zero),
     _rho_w(_mesh.dimension() == 3 ? coupledValue(nms::momentum_z) : _zero),
     _rho_et(coupledValue(nms::total_energy_density)),
 
     // Time derivative values
     _drho_dt(coupledDot(nms::density)),
     _drhou_dt(coupledDot(nms::momentum_x)),
     _drhov_dt(_mesh.dimension() >= 2 ? coupledDot(nms::momentum_y) : _zero),
     _drhow_dt(_mesh.dimension() == 3 ? coupledDot(nms::momentum_z) : _zero),
     _drhoE_dt(coupledDot(nms::total_energy_density)),
 
     // Gradients
     _grad_rho(coupledGradient(nms::density)),
     _grad_rho_u(coupledGradient(nms::momentum_x)),
     _grad_rho_v(_mesh.dimension() >= 2 ? coupledGradient(nms::momentum_y) : _grad_zero),
     _grad_rho_w(_mesh.dimension() == 3 ? coupledGradient(nms::momentum_z) : _grad_zero),
     _grad_rho_et(coupledGradient(nms::total_energy_density)),
 
     // Material properties for stabilization
     _hsupg(declareProperty<Real>("hsupg")),
     _tauc(declareProperty<Real>("tauc")),
     _taum(declareProperty<Real>("taum")),
     _taue(declareProperty<Real>("taue")),
     _strong_residuals(declareProperty<std::vector<Real>>("strong_residuals")),
     _fp(getUserObject<IdealGasFluidProperties>("fluid_properties"))
 {
 }

Member Function Documentation

◆ computeHSUPG()

void NavierStokesMaterial::computeHSUPG ( unsigned int qp )

private

Definition at line 179 of file NavierStokesMaterial.C.

Referenced by computeProperties().

 {
   // // Grab reference to linear Lagrange finite element object pointer,
   // // currently this is always a linear Lagrange element, so this might need to
   // // be generalized if we start working with higher-order elements...
   // FEBase*& fe(_assembly.getFE(FEType(), _current_elem->dim()));
   //
   // // Grab references to FE object's mapping data from the _subproblem's FE object
   // const std::vector<Real> & dxidx(fe->get_dxidx());
   // const std::vector<Real> & dxidy(fe->get_dxidy());
   // const std::vector<Real> & dxidz(fe->get_dxidz());
   // const std::vector<Real> & detadx(fe->get_detadx());
   // const std::vector<Real> & detady(fe->get_detady());
   // const std::vector<Real> & detadz(fe->get_detadz());
   // const std::vector<Real> & dzetadx(fe->get_dzetadx()); // Empty in 2D
   // const std::vector<Real> & dzetady(fe->get_dzetady()); // Empty in 2D
   // const std::vector<Real> & dzetadz(fe->get_dzetadz()); // Empty in 2D
   //
   // // Bounds checking on element data
   // mooseAssert(qp < dxidx.size(), "Insufficient data in dxidx array!");
   // mooseAssert(qp < dxidy.size(), "Insufficient data in dxidy array!");
   // mooseAssert(qp < dxidz.size(), "Insufficient data in dxidz array!");
   //
   // mooseAssert(qp < detadx.size(), "Insufficient data in detadx array!");
   // mooseAssert(qp < detady.size(), "Insufficient data in detady array!");
   // mooseAssert(qp < detadz.size(), "Insufficient data in detadz array!");
   //
   // if (_mesh_dimension == 3)
   // {
   //   mooseAssert(qp < dzetadx.size(), "Insufficient data in dzetadx array!");
   //   mooseAssert(qp < dzetady.size(), "Insufficient data in dzetady array!");
   //   mooseAssert(qp < dzetadz.size(), "Insufficient data in dzetadz array!");
   // }
   //
   // // The velocity vector at this quadrature point.
   // RealVectorValue U(_u_vel[qp],_v_vel[qp],_w_vel[qp]);
   //
   // // Pull out element inverse map values at the current qp into a little dense matrix
   // Real dxi_dx[3][3] = {{0.,0.,0.}, {0.,0.,0.}, {0.,0.,0.}};
   //
   // dxi_dx[0][0] = dxidx[qp];  dxi_dx[0][1] = dxidy[qp];
   // dxi_dx[1][0] = detadx[qp]; dxi_dx[1][1] = detady[qp];
   //
   // // OK to access third entries on 2D elements if LIBMESH_DIM==3, though they
   // // may be zero...
   // if (LIBMESH_DIM == 3)
   // {
   //   /**/             /**/               dxi_dx[0][2] = dxidz[qp];
   //   /**/             /**/               dxi_dx[1][2] = detadz[qp];
   // }
   //
   // // The last row of entries available only for 3D elements.
   // if (_mesh_dimension == 3)
   // {
   //   dxi_dx[2][0] = dzetadx[qp];   dxi_dx[2][1] = dzetady[qp];   dxi_dx[2][2] = dzetadz[qp];
   // }
   //
   // // Construct the g_ij = d(xi_k)/d(x_j) * d(xi_k)/d(x_i) matrix
   // // from Ben and Bova's paper by summing over k...
   // Real g[3][3] = {{0.,0.,0.}, {0.,0.,0.}, {0.,0.,0.}};
   // for (unsigned int i = 0; i < 3; ++i)
   //   for (unsigned int j = 0; j < 3; ++j)
   //     for (unsigned int k = 0; k < 3; ++k)
   //       g[i][j] += dxi_dx[k][j] * dxi_dx[k][i];
   //
   // // Compute the denominator of the h_supg term: U * (g) * U
   // Real denom = 0.;
   // for (unsigned int i = 0; i < 3; ++i)
   //   for (unsigned int j = 0; j < 3; ++j)
   //     denom += U(j) * g[i][j] * U(i);
   //
   // // Compute h_supg.  Some notes:
   // // .) The 2 coefficient in this term should be a 1 if we are using tets/triangles.
   // // .) The denominator will be identically zero only if the velocity
   // //    is identically zero, in which case we can't divide by it.
   // if (denom != 0.0)
   //   _hsupg[qp] = 2.* sqrt( U.norm_sq() / denom );
   // else
   //   _hsupg[qp] = 0.;
 
   // Simple (and fast) implementation: Just use hmin for the element!
   _hsupg[qp] = _current_elem->hmin();
 }

◆ computeProperties()

void NavierStokesMaterial::computeProperties ( )

protectedvirtual

Must be called after the child class computes dynamic_viscocity.

Must be called after the child class computes dynamic_viscosity.

Reimplemented from Material.

Reimplemented in Air.

Definition at line 117 of file NavierStokesMaterial.C.

Referenced by Air::computeProperties().

 {
   for (unsigned int qp = 0; qp < _qrule->n_points(); ++qp)
   {
     /******* Viscous Stress Tensor *******/
     // Technically... this _is_ the transpose (since we are loading these by rows)
     // But it doesn't matter....
     RealTensorValue grad_outer_u(_grad_u[qp], _grad_v[qp], _grad_w[qp]);
 
     grad_outer_u += grad_outer_u.transpose();
 
     Real div_vel = 0.0;
     for (unsigned int i = 0; i < 3; ++i)
       div_vel += (*_vel_grads[i])[qp](i);
 
     // Add diagonal terms
     for (unsigned int i = 0; i < 3; ++i)
       grad_outer_u(i, i) -= 2.0 / 3.0 * div_vel;
 
     grad_outer_u *= _dynamic_viscosity[qp];
 
     _viscous_stress_tensor[qp] = grad_outer_u;
 
     // Tabulated values of thermal conductivity vs. Temperature for air (k increases slightly with
     // T):
     // T (K)    k (W/m-K)
     // 273      0.0243
     // 373      0.0314
     // 473      0.0386
     // 573      0.0454
     // 673      0.0515
 
     // Pr = (mu * cp) / k  ==>  k = (mu * cp) / Pr = (mu * gamma * cv) / Pr.
     // TODO: We are using a fixed value of the Prandtl number which is
     // valid for air, it may or may not depend on temperature?  Since
     // this is a property of the fluid, it could possibly be moved to
     // the FluidProperties module...
     const Real Pr = 0.71;
     _thermal_conductivity[qp] = (_dynamic_viscosity[qp] * _fp.cp()) / Pr;
 
     // Compute stabilization parameters:
 
     // .) Compute SUPG element length scale.
     computeHSUPG(qp);
     // Moose::out << "_hsupg[" << qp << "]=" << _hsupg[qp] << std::endl;
 
     // .) Compute SUPG parameter values.  (Must call this after computeHSUPG())
     computeTau(qp);
     // Moose::out << "_tauc[" << qp << "]=" << _tauc[qp] << ", ";
     // Moose::out << "_taum[" << qp << "]=" << _taum[qp] << ", ";
     // Moose::out << "_taue[" << qp << "]=" << _taue[qp] << std::endl;
 
     // .) Compute strong residual values.
     computeStrongResiduals(qp);
     // Moose::out << "_strong_residuals[" << qp << "]=";
     // for (unsigned i=0; i<_strong_residuals[qp].size(); ++i)
     //   Moose::out << _strong_residuals[qp][i] << " ";
     // Moose::out << std::endl;
   }
 }

◆ computeStrongResiduals()

void NavierStokesMaterial::computeStrongResiduals ( unsigned int qp )

private

Definition at line 354 of file NavierStokesMaterial.C.

Referenced by computeProperties().

 {
   // Create storage at this qp for the strong residuals of all the equations.
   // In 2D, the value for the z-velocity equation will just be zero.
   _strong_residuals[qp].resize(5);
 
   // The timestep is stored in the Problem object, which can be accessed through
   // the parent pointer of the SubProblem.  Don't need this if we are not
   // approximating time derivatives ourselves.
   // Real dt = _subproblem.parent()->dt();
   // Moose::out << "dt=" << dt << std::endl;
 
   // Vector object for the velocity
   RealVectorValue vel(_u_vel[qp], _v_vel[qp], _w_vel[qp]);
 
   // A VectorValue object containing all zeros.  Makes it easier to
   // construct type tensor objects
   RealVectorValue zero(0., 0., 0.);
 
   // Velocity vector magnitude squared
   Real velmag2 = vel.norm_sq();
 
   // Debugging: How large are the time derivative parts of the strong residuals?
   //  Moose::out << "drho_dt=" << _drho_dt
   //            << ", drhou_dt=" << _drhou_dt
   //            << ", drhov_dt=" << _drhov_dt
   //            << ", drhow_dt=" << _drhow_dt
   //            << ", drhoE_dt=" << _drhoE_dt
   //            << std::endl;
 
   // Momentum divergence
   Real divU = _grad_rho_u[qp](0) + _grad_rho_v[qp](1) + _grad_rho_w[qp](2);
 
   // Enough space to hold three space dimensions of velocity components at each qp,
   // regardless of what dimension we are actually running in.
   _calC[qp].resize(3);
 
   // Explicitly zero the calC
   for (unsigned int i = 0; i < 3; ++i)
     _calC[qp][i].zero();
 
   // x-column matrix
   _calC[qp][0](0, 0) = _u_vel[qp];
   _calC[qp][0](1, 0) = _v_vel[qp];
   _calC[qp][0](2, 0) = _w_vel[qp];
 
   // y-column matrix
   _calC[qp][1](0, 1) = _u_vel[qp];
   _calC[qp][1](1, 1) = _v_vel[qp];
   _calC[qp][1](2, 1) = _w_vel[qp];
 
   // z-column matrix (this assumes LIBMESH_DIM==3!)
   _calC[qp][2](0, 2) = _u_vel[qp];
   _calC[qp][2](1, 2) = _v_vel[qp];
   _calC[qp][2](2, 2) = _w_vel[qp];
 
   // The matrix S can be computed from any of the calC via calC_1*calC_1^T
   RealTensorValue calS = _calC[qp][0] * _calC[qp][0].transpose();
 
   // Enough space to hold five (=n_sd + 2) 3*3 calA matrices at this qp, regarless of dimension
   _calA[qp].resize(5);
 
   // 0.) _calA_0 = diag( (gam - 1)/2*|u|^2 ) - S
   _calA[qp][0].zero(); // zero this calA entry
   _calA[qp][0](0, 0) = _calA[qp][0](1, 1) = _calA[qp][0](2, 2) =
       0.5 * (_fp.gamma() - 1.0) * velmag2; // set diag. entries
   _calA[qp][0] -= calS;
 
   for (unsigned int m = 1; m <= 3; ++m)
   {
     // Use m_local when indexing into matrices and vectors
     unsigned int m_local = m - 1;
 
     // For m=1,2,3, calA_m = C_m + C_m^T + diag( (1.-gam)*u_m )
     _calA[qp][m].zero(); // zero this calA entry
     _calA[qp][m](0, 0) = _calA[qp][m](1, 1) = _calA[qp][m](2, 2) =
         (1. - _fp.gamma()) * vel(m_local);          // set diag. entries
     _calA[qp][m] += _calC[qp][m_local];             // Note: use m_local for indexing into _calC!
     _calA[qp][m] += _calC[qp][m_local].transpose(); // Note: use m_local for indexing into _calC!
   }
 
   // 4.) calA_4 = diag(gam - 1)
   _calA[qp][4].zero(); // zero this calA entry
   _calA[qp][4](0, 0) = _calA[qp][4](1, 1) = _calA[qp][4](2, 2) = (_fp.gamma() - 1.0);
 
   // Enough space to hold the 3*5 "cal E" matrices which comprise the inviscid flux term
   // of the energy equation.  See notes for additional details
   _calE[qp].resize(3); // Three rows, 5 entries in each row
 
   for (unsigned int k = 0; k < 3; ++k)
   {
     // Make enough room to store all 5 E matrices for this k
     _calE[qp][k].resize(5);
 
     // Store and reuse the velocity column transpose matrix for the
     // current value of k.
     RealTensorValue Ck_T = _calC[qp][k].transpose();
 
     // E_{k0} (density gradient term)
     _calE[qp][k][0].zero();
     _calE[qp][k][0] = (0.5 * (_fp.gamma() - 1.0) * velmag2 - _specific_total_enthalpy[qp]) * Ck_T;
 
     for (unsigned int m = 1; m <= 3; ++m)
     {
       // Use m_local when indexing into matrices and vectors
       unsigned int m_local = m - 1;
 
       // E_{km} (momentum gradient terms)
       _calE[qp][k][m].zero();
       _calE[qp][k][m](k, m_local) = _specific_total_enthalpy[qp];  // H * D_{km}
       _calE[qp][k][m] += (1. - _fp.gamma()) * vel(m_local) * Ck_T; // (1-gam) * u_m * C_k^T
     }
 
     // E_{k4} (energy gradient term)
     _calE[qp][k][4].zero();
     _calE[qp][k][4] = _fp.gamma() * Ck_T;
   }
 
   // Compute the sum over ell of: A_ell grad(U_ell), store in DenseVector or Gradient object?
   // The gradient object might be more useful, since we are multiplying by VariableGradient
   // (which is a MooseArray of RealGradients) objects?
   RealVectorValue mom_resid = _calA[qp][0] * _grad_rho[qp] + _calA[qp][1] * _grad_rho_u[qp] +
                               _calA[qp][2] * _grad_rho_v[qp] + _calA[qp][3] * _grad_rho_w[qp] +
                               _calA[qp][4] * _grad_rho_et[qp];
 
   // No matrices/vectors for the energy residual strong form... just write it out like
   // the mass equation residual.  See "Momentum SUPG terms prop. to energy residual"
   // section of the notes.
   Real energy_resid =
       (0.5 * (_fp.gamma() - 1.0) * velmag2 - _specific_total_enthalpy[qp]) * (vel * _grad_rho[qp]) +
       _specific_total_enthalpy[qp] * divU +
       (1. - _fp.gamma()) * (vel(0) * (vel * _grad_rho_u[qp]) + vel(1) * (vel * _grad_rho_v[qp]) +
                             vel(2) * (vel * _grad_rho_w[qp])) +
       _fp.gamma() * (vel * _grad_rho_et[qp]);
 
   // Now for the actual residual values...
 
   // The density strong-residual
   _strong_residuals[qp][0] = _drho_dt[qp] + divU;
 
   // The x-momentum strong-residual, viscous terms neglected.
   // TODO: If we want to add viscous contributions back in, should this kernel
   // not inherit from NSViscousFluxBase so it can get tau values?  This would
   // also involve shape function second derivative values.
   _strong_residuals[qp][1] = _drhou_dt[qp] + mom_resid(0);
 
   // The y-momentum strong residual, viscous terms neglected.
   _strong_residuals[qp][2] = _drhov_dt[qp] + mom_resid(1);
 
   // The z-momentum strong residual, viscous terms neglected.
   if (_mesh_dimension == 3)
     _strong_residuals[qp][3] = _drhow_dt[qp] + mom_resid(2);
   else
     _strong_residuals[qp][3] = 0.;
 
   // The energy equation strong residual
   _strong_residuals[qp][4] = _drhoE_dt[qp] + energy_resid;
 }

◆ computeTau()

void NavierStokesMaterial::computeTau ( unsigned int qp )

private

Definition at line 264 of file NavierStokesMaterial.C.

Referenced by computeProperties().

 {
   Real velmag =
       std::sqrt(_u_vel[qp] * _u_vel[qp] + _v_vel[qp] * _v_vel[qp] + _w_vel[qp] * _w_vel[qp]);
 
   // Moose::out << "velmag=" << velmag << std::endl;
 
   // Make sure temperature >= 0 before trying to take sqrt
   // if (_temperature[qp] < 0.)
   // {
   //   Moose::err << "Negative temperature "
   //             << _temperature[qp]
   //             << " found at quadrature point "
   //             << qp
   //             << ", element "
   //             << _current_elem->id()
   //             << std::endl;
   //   mooseError("Can't continue, would be nice to throw an exception here?");
   // }
 
   // The speed of sound for an ideal gas, sqrt(gamma * R * T).  Not needed unless
   // we want to use a form of Tau that requires it.
   // Real soundspeed = _fp.c_from_v_e(_specific_volume[_qp], _internal_energy[_qp]);
 
   // If velmag == 0, then _hsupg should be zero as well.  Then tau
   // will have only the time-derivative contribution (or zero, if we
   // are not including dt terms in our taus!)  Note that using the
   // time derivative contribution in this way assumes we are solving
   // unsteady, and guarantees *some* stabilization is added even when
   // u -> 0 in certain regions of the flow.
   if (velmag == 0.)
   {
     // 1.) Tau without dt terms
     // _tauc[qp] = 0.;
     // _taum[qp] = 0.;
     // _taue[qp] = 0.;
 
     // 2.) Tau *with* dt terms
     _tauc[qp] = _taum[qp] = _taue[qp] = 0.5 * _dt;
   }
   else
   {
     // The element length parameter, squared
     Real h2 = _hsupg[qp] * _hsupg[qp];
 
     // The viscosity-based term
     Real visc_term = _dynamic_viscosity[qp] / _rho[qp] / h2;
 
     // The thermal conductivity-based term, cp = gamma * cv
     Real k_term = _thermal_conductivity[qp] / _rho[qp] / _fp.cp() / h2;
 
     // 1a.) Standard compressible flow tau.  Does not account for low Mach number
     // limit.
     //    _tauc[qp] = _hsupg[qp] / (velmag + soundspeed);
 
     // 1b.) Inspired by Hauke, the sum of the compressible and incompressible tauc.
     //    _tauc[qp] =
     //      _hsupg[qp] / (velmag + soundspeed) +
     //      _hsupg[qp] / (velmag);
 
     // 1c.) From Wong 2001.  This tau is O(M^2) for small M.  At small M,
     // tauc dominates the inverse square sums and basically makes
     // taum=taue=tauc.  However, all my flows occur at low Mach numbers,
     // so there would basically never be any stabilization...
     // _tauc[qp] = (_hsupg[qp] * velmag) / (velmag*velmag + soundspeed*soundspeed);
 
     // For use with option "1",
     // (tau_c)^{-2}
     //    Real taucm2 = 1./_tauc[qp]/_tauc[qp];
     //    _taum[qp] = 1. / std::sqrt(taucm2 + visc_term*visc_term);
     //    _taue[qp] = 1. / std::sqrt(taucm2 + k_term*k_term);
 
     // 2.) Tau with timestep dependence (guarantees stabilization even
     // in zero-velocity limit) incorporated via the "r-switch" method,
     // with r=2.
     Real sqrt_term = 4. / _dt / _dt + velmag * velmag / h2;
 
     // For use with option "2", i.e. the option that uses dt in the definition of tau
     _tauc[qp] = 1. / std::sqrt(sqrt_term);
     _taum[qp] = 1. / std::sqrt(sqrt_term + visc_term * visc_term);
     _taue[qp] = 1. / std::sqrt(sqrt_term + k_term * k_term);
   }
 
   // Debugging
   // Moose::out << "_tauc[" << qp << "]=" << _tauc[qp] << std::endl;
   // Moose::out << "_hsupg[" << qp << "]=" << _hsupg[qp] << std::endl;
   // Moose::out << "velmag[" << qp << "]=" << velmag << std::endl;
 }

◆ validParams()

InputParameters NavierStokesMaterial::validParams ( )

static

Definition at line 26 of file NavierStokesMaterial.C.

Referenced by Air::validParams().

 {
   InputParameters params = Material::validParams();
 
   params.addClassDescription(
       "This is the base class of all materials should use if you are trying to "
       "use the Navier-Stokes Kernels.");
   params.addRequiredCoupledVar(nms::velocity_x, "x-velocity");
   params.addCoupledVar(nms::velocity_y, "y-velocity"); // only required in >= 2D
   params.addCoupledVar(nms::velocity_z, "z-velocity"); // only required in 3D
 
   params.addRequiredCoupledVar(nms::temperature, "temperature");
   params.addRequiredCoupledVar(nms::specific_total_enthalpy, "specific total enthalpy");
 
   params.addRequiredCoupledVar(nms::density, "density");
   params.addRequiredCoupledVar(nms::momentum_x, "x-momentum");
   params.addCoupledVar(nms::momentum_y, "y-momentum"); // only required in >= 2D
   params.addCoupledVar(nms::momentum_z, "z-momentum"); // only required in 3D
   params.addRequiredCoupledVar(nms::total_energy_density, "energy");
   params.addRequiredParam<UserObjectName>("fluid_properties",
                                           "The name of the user object for fluid properties");
   return params;
 }

Member Data Documentation

◆ _calA

MaterialProperty<std::vector<RealTensorValue> >& NavierStokesMaterial::_calA

protected

Definition at line 59 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _calC

MaterialProperty<std::vector<RealTensorValue> >& NavierStokesMaterial::_calC

protected

Definition at line 63 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _calE

MaterialProperty<std::vector<std::vector<RealTensorValue> > >& NavierStokesMaterial::_calE

protected

Definition at line 68 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _drho_dt

const VariableValue& NavierStokesMaterial::_drho_dt

protected

Definition at line 94 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _drhoE_dt

const VariableValue& NavierStokesMaterial::_drhoE_dt

protected

Definition at line 98 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _drhou_dt

const VariableValue& NavierStokesMaterial::_drhou_dt

protected

Definition at line 95 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _drhov_dt

const VariableValue& NavierStokesMaterial::_drhov_dt

protected

Definition at line 96 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _drhow_dt

const VariableValue& NavierStokesMaterial::_drhow_dt

protected

Definition at line 97 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _dynamic_viscosity

MaterialProperty<Real>& NavierStokesMaterial::_dynamic_viscosity

protected

Definition at line 53 of file NavierStokesMaterial.h.

Referenced by Air::computeProperties(), computeProperties(), and computeTau().

◆ _fp

const IdealGasFluidProperties& NavierStokesMaterial::_fp

protected

Definition at line 119 of file NavierStokesMaterial.h.

Referenced by computeProperties(), computeStrongResiduals(), and computeTau().

◆ _grad_rho

const VariableGradient& NavierStokesMaterial::_grad_rho

protected

Definition at line 101 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _grad_rho_et

const VariableGradient& NavierStokesMaterial::_grad_rho_et

protected

Definition at line 105 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _grad_rho_u

const VariableGradient& NavierStokesMaterial::_grad_rho_u

protected

Definition at line 102 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _grad_rho_v

const VariableGradient& NavierStokesMaterial::_grad_rho_v

protected

Definition at line 103 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _grad_rho_w

const VariableGradient& NavierStokesMaterial::_grad_rho_w

protected

Definition at line 104 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _grad_u

const VariableGradient& NavierStokesMaterial::_grad_u

protected

Definition at line 47 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and NavierStokesMaterial().

◆ _grad_v

const VariableGradient& NavierStokesMaterial::_grad_v

protected

Definition at line 48 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and NavierStokesMaterial().

◆ _grad_w

const VariableGradient& NavierStokesMaterial::_grad_w

protected

Definition at line 49 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and NavierStokesMaterial().

◆ _hsupg

MaterialProperty<Real>& NavierStokesMaterial::_hsupg

protected

Definition at line 109 of file NavierStokesMaterial.h.

Referenced by computeHSUPG(), and computeTau().

◆ _mesh_dimension

const unsigned int NavierStokesMaterial::_mesh_dimension

protected

Definition at line 45 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _rho

const VariableValue& NavierStokesMaterial::_rho

protected

Definition at line 87 of file NavierStokesMaterial.h.

Referenced by computeTau().

◆ _rho_et

const VariableValue& NavierStokesMaterial::_rho_et

protected

Definition at line 91 of file NavierStokesMaterial.h.

◆ _rho_u

const VariableValue& NavierStokesMaterial::_rho_u

protected

Definition at line 88 of file NavierStokesMaterial.h.

◆ _rho_v

const VariableValue& NavierStokesMaterial::_rho_v

protected

Definition at line 89 of file NavierStokesMaterial.h.

◆ _rho_w

const VariableValue& NavierStokesMaterial::_rho_w

protected

Definition at line 90 of file NavierStokesMaterial.h.

◆ _specific_total_enthalpy

const VariableValue& NavierStokesMaterial::_specific_total_enthalpy

protected

Definition at line 84 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _strong_residuals

MaterialProperty<std::vector<Real> >& NavierStokesMaterial::_strong_residuals

protected

Definition at line 116 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals().

◆ _tauc

MaterialProperty<Real>& NavierStokesMaterial::_tauc

protected

Definition at line 110 of file NavierStokesMaterial.h.

Referenced by computeTau().

◆ _taue

MaterialProperty<Real>& NavierStokesMaterial::_taue

protected

Definition at line 112 of file NavierStokesMaterial.h.

Referenced by computeTau().

◆ _taum

MaterialProperty<Real>& NavierStokesMaterial::_taum

protected

Definition at line 111 of file NavierStokesMaterial.h.

Referenced by computeTau().

◆ _temperature

const VariableValue& NavierStokesMaterial::_temperature

protected

Definition at line 81 of file NavierStokesMaterial.h.

◆ _thermal_conductivity

MaterialProperty<Real>& NavierStokesMaterial::_thermal_conductivity

protected

Definition at line 52 of file NavierStokesMaterial.h.

Referenced by computeProperties(), and computeTau().

◆ _u_vel

const VariableValue& NavierStokesMaterial::_u_vel

protected

Definition at line 76 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals(), and computeTau().

◆ _v_vel

const VariableValue& NavierStokesMaterial::_v_vel

protected

Definition at line 77 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals(), and computeTau().

◆ _vel_grads

std::vector<const VariableGradient *> NavierStokesMaterial::_vel_grads

protected

Definition at line 72 of file NavierStokesMaterial.h.

Referenced by computeProperties().

◆ _viscous_stress_tensor

MaterialProperty<RealTensorValue>& NavierStokesMaterial::_viscous_stress_tensor

protected

Definition at line 51 of file NavierStokesMaterial.h.

Referenced by computeProperties().

◆ _w_vel

const VariableValue& NavierStokesMaterial::_w_vel

protected

Definition at line 78 of file NavierStokesMaterial.h.

Referenced by computeStrongResiduals(), and computeTau().

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

navier_stokes/include/materials/NavierStokesMaterial.h
navier_stokes/src/materials/NavierStokesMaterial.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ NavierStokesMaterial()

Member Function Documentation

◆ computeHSUPG()

◆ computeProperties()

◆ computeStrongResiduals()

◆ computeTau()

◆ validParams()

Member Data Documentation

◆ _calA

◆ _calC

◆ _calE

◆ _drho_dt

◆ _drhoE_dt

◆ _drhou_dt

◆ _drhov_dt

◆ _drhow_dt

◆ _dynamic_viscosity

◆ _fp

◆ _grad_rho

◆ _grad_rho_et

◆ _grad_rho_u

◆ _grad_rho_v

◆ _grad_rho_w

◆ _grad_u

◆ _grad_v

◆ _grad_w

◆ _hsupg

◆ _mesh_dimension

◆ _rho

◆ _rho_et

◆ _rho_u

◆ _rho_v

◆ _rho_w

◆ _specific_total_enthalpy

◆ _strong_residuals

◆ _tauc

◆ _taue

◆ _taum

◆ _temperature

◆ _thermal_conductivity

◆ _u_vel

◆ _v_vel

◆ _vel_grads

◆ _viscous_stress_tensor

◆ _w_vel