This class computes strong and weak components of the INS governing equations. More...

#include <INSBase.h>

Inheritance diagram for INSBase:

Public Types
enum	ResidualTagType { ResidualTagType::NonReference, ResidualTagType::Reference }

typedef DataFileName	DataFileParameterType

Public Member Functions
	INSBase (const InputParameters &parameters)

virtual	~INSBase ()

virtual void	computeResidual () override

virtual void	computeJacobian () override

virtual void	computeOffDiagJacobian (unsigned int jvar) override

virtual void	computeResidualAndJacobian () override

virtual void	computeOffDiagJacobianScalar (unsigned int jvar) override

virtual const MooseVariable &	variable () const override

virtual void	computeNonlocalJacobian ()

virtual void	computeNonlocalOffDiagJacobian (unsigned int)

const SubProblem &	subProblem () const

virtual void	prepareShapes (unsigned int var_num)

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 &	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 &name, const std::string *param=nullptr) const

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	subdomainSetup ()

virtual void	customSetup (const ExecFlagType &)

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

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

void	setRandomResetFrequency (ExecFlagType exec_flag)

unsigned long	getRandomLong () const

Real	getRandomReal () const

unsigned int	getSeed (std::size_t id)

unsigned int	getMasterSeed () const

bool	isNodal () const

ExecFlagType	getResetOnTime () const

void	setRandomDataPointer (RandomData *random_data)

virtual void	meshChanged ()

void	useVectorTag (const TagName &tag_name, VectorTagsKey)

void	useVectorTag (TagID tag_id, VectorTagsKey)

void	useMatrixTag (const TagName &tag_name, MatrixTagsKey)

void	useMatrixTag (TagID tag_id, MatrixTagsKey)

bool	isVectorTagged ()

bool	isMatrixTagged ()

bool	hasVectorTags () const

const std::set< TagID > &	getVectorTags (VectorTagsKey) const

const std::set< TagID > &	getMatrixTags (MatrixTagsKey) const

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

unsigned int	numBlocks () const

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

unsigned int	blocksMaxDimension () const

bool	hasBlocks (const SubdomainName &name) const

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

bool	hasBlocks (SubdomainID id) const

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

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

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

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

bool	hasBlockMaterialProperty (const std::string &prop_name)

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

virtual bool	blockRestricted () const

virtual void	checkVariable (const MooseVariableFieldBase &variable) const

const std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> > &	getCoupledVars () const

const std::vector< MooseVariableFieldBase *> &	getCoupledMooseVars () const

const std::vector< MooseVariable *> &	getCoupledStandardMooseVars () const

const std::vector< VectorMooseVariable *> &	getCoupledVectorMooseVars () const

const std::vector< ArrayMooseVariable *> &	getCoupledArrayMooseVars () const

void	addFEVariableCoupleableVectorTag (TagID tag)

void	addFEVariableCoupleableMatrixTag (TagID tag)

std::set< TagID > &	getFEVariableCoupleableVectorTags ()

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

std::set< TagID > &	getFEVariableCoupleableMatrixTags ()

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

auto &	getWritableCoupledVariables () const

bool	hasWritableCoupledVariables () const

const 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 std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()

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

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

const std::set< MooseVariableFieldBase *> &	getMooseVariableDependencies () const

std::set< MooseVariableFieldBase *>	checkAllVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_omit={})

std::set< MooseVariableFieldBase *>	checkVariables (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)

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

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

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)

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

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

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)

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

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

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

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

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

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

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

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

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)

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

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

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

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

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

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

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

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

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

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

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

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

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

std::set< 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 std::unordered_set< unsigned int > &	getMatPropDependencies () const

virtual void	resolveOptionalProperties ()

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

virtual unsigned int	getElementIDIndex (const std::string &id_parameter_name, unsigned int comp=0) const

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

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

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

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

virtual const dof_id_type &	getElementIDByName (const std::string &id_name) const

virtual const dof_id_type &	getElementIDNeighborByName (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

MooseVariableBase *	mooseVariableBase () const

MooseVariableField< Real > &	mooseVariableField ()

MooseVariableFE< Real > *	mooseVariable () const

MooseVariableFV< Real > *	mooseVariableFV () const

MooseLinearVariableFV< Real > *	mooseLinearVariableFV () const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

const 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

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)

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

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

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

MaterialBase &	getMaterialByName (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)

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 ()

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 Real	computeQpResidual ()=0

virtual Real	computeQpJacobian ()=0

virtual Real	computeQpOffDiagJacobian (unsigned jvar)=0

virtual RealVectorValue	convectiveTerm ()

virtual RealVectorValue	dConvecDUComp (unsigned comp)

virtual RealVectorValue	strongViscousTermLaplace ()

virtual RealVectorValue	strongViscousTermTraction ()

virtual RealVectorValue	dStrongViscDUCompLaplace (unsigned comp)

virtual RealVectorValue	dStrongViscDUCompTraction (unsigned comp)

virtual RealVectorValue	weakViscousTermLaplace (unsigned comp)

virtual RealVectorValue	weakViscousTermTraction (unsigned comp)

virtual RealVectorValue	dWeakViscDUCompLaplace ()

virtual RealVectorValue	dWeakViscDUCompTraction ()

virtual RealVectorValue	strongPressureTerm ()

virtual Real	weakPressureTerm ()

virtual RealVectorValue	dStrongPressureDPressure ()

virtual Real	dWeakPressureDPressure ()

virtual RealVectorValue	gravityTerm ()

virtual RealVectorValue	timeDerivativeTerm ()

virtual RealVectorValue	dTimeDerivativeDUComp (unsigned comp)

virtual Real	tau ()

virtual Real	dTauDUComp (unsigned comp)

virtual Real	tauNodal ()
	Provides tau which yields superconvergence for 1D advection-diffusion. More...

RealVectorValue	relativeVelocity () const
	Compute the velocity. More...

RealVectorValue	strongViscousTermLaplaceRZ () const
	Computes the additional RZ terms for the Laplace form of the strong viscous term. More...

RealVectorValue	dStrongViscDUCompLaplaceRZ (const unsigned int comp) const
	Computes the Jacobian for the additional RZ terms for the Laplace form of the strong viscous term for the given velocity component `comp`. More...

RealVectorValue	strongViscousTermTractionRZ () const
	Computes the additional RZ terms for the Traction form of the strong viscous term. More...

RealVectorValue	dStrongViscDUCompTractionRZ (const unsigned int comp) const
	Computes the Jacobian for the additional RZ terms for the Traction form of the strong viscous term for the given velocity component `comp`. More...

virtual Real	computeQpOffDiagJacobian (unsigned int)

virtual Real	computeQpOffDiagJacobianScalar (unsigned int)

virtual RealEigenVector	computeQpOffDiagJacobianArray (const ArrayMooseVariable &jvar)

virtual void	precalculateResidual ()

virtual void	precalculateJacobian ()

virtual void	precalculateOffDiagJacobian (unsigned int)

const MooseVariableFieldBase &	getVariable (unsigned int jvar_num) const

virtual void	addUserObjectDependencyHelper (const UserObject &) const

Moose::StateArg	determineState () const

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &) const

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const

T &	declareRestartableData (const std::string &data_name, Args &&... args)

ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

const T &	getRestartableData (const std::string &data_name) const

T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

T &	declareRecoverableData (const std::string &data_name, Args &&... args)

T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)

T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)

std::string	restartableName (const std::string &data_name) const

void	prepareVectorTag (Assembly &assembly, unsigned int ivar)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)

void	prepareVectorTagNeighbor (Assembly &assembly, unsigned int ivar)

void	prepareVectorTagLower (Assembly &assembly, unsigned int ivar)

void	prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar)

void	prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar, DenseMatrix< Number > &k) const

void	prepareMatrixTagNonlocal (Assembly &assembly, unsigned int ivar, unsigned int jvar)

void	prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)

void	prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, DenseMatrix< Number > &k) const

void	prepareMatrixTagLower (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)

void	accumulateTaggedLocalResidual ()

void	assignTaggedLocalResidual ()

void	accumulateTaggedLocalMatrix ()

void	accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, const DenseMatrix< Number > &k)

void	accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, const DenseMatrix< Number > &k)

void	accumulateTaggedNonlocalMatrix ()

void	assignTaggedLocalMatrix ()

void	addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResiduals (Assembly &assembly, const DenseVector< T > &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobian (Assembly &assembly, DenseMatrix< Real > &local_k, const std::vector< dof_id_type > &row_indices, const std::vector< dof_id_type > &column_indices, Real scaling_factor)

void	addResidualsWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addResidualsAndJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)

void	addJacobianElement (Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)

void	setResidual (SystemBase &sys, const T &residual, MooseVariableFE< T > &var)

void	setResidual (SystemBase &sys, Real residual, dof_id_type dof_index)

void	setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)

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

void	initializeBlockRestrictable (const MooseObject *moose_object)

Moose::CoordinateSystemType	getBlockCoordSystem ()

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

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

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

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

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

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

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

std::vector< unsigned 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

std::vector< const GenericVariableValue< is_ad > *>	coupledGenericValues (const std::string &var_name) const

std::vector< const GenericVariableValue< false > *>	coupledGenericValues (const std::string &var_name) const

std::vector< const GenericVariableValue< true > *>	coupledGenericValues (const std::string &var_name) const

const GenericVariableValue< is_ad > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericDofValue (const std::string &var_name, unsigned int comp) const

virtual const 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

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

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

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

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

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

Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const

const 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

void	addMooseVariableDependency (MooseVariableFieldBase *var)

void	addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)

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

void	markMatPropRequested (const std::string &)

MaterialPropertyName	getMaterialPropertyName (const std::string &name) const

void	checkExecutionStage ()

virtual const OutputTools< Real >::VariableValue &	value ()

virtual const OutputTools< Real >::VariableValue &	valueOld ()

virtual const OutputTools< Real >::VariableValue &	valueOlder ()

virtual const OutputTools< Real >::VariableValue &	dot ()

virtual const OutputTools< Real >::VariableValue &	dotDot ()

virtual const OutputTools< Real >::VariableValue &	dotOld ()

virtual const OutputTools< Real >::VariableValue &	dotDotOld ()

virtual const VariableValue &	dotDu ()

virtual const VariableValue &	dotDotDu ()

virtual const OutputTools< Real >::VariableGradient &	gradient ()

virtual const OutputTools< Real >::VariableGradient &	gradientOld ()

virtual const OutputTools< Real >::VariableGradient &	gradientOlder ()

virtual const OutputTools< Real >::VariableSecond &	second ()

virtual const OutputTools< Real >::VariableSecond &	secondOld ()

virtual const OutputTools< Real >::VariableSecond &	secondOlder ()

virtual const OutputTools< Real >::VariableTestSecond &	secondTest ()

virtual const OutputTools< Real >::VariableTestSecond &	secondTestFace ()

virtual const OutputTools< Real >::VariablePhiSecond &	secondPhi ()

virtual const OutputTools< Real >::VariablePhiSecond &	secondPhiFace ()

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 VariablePhiSecond &	_second_phi
	second derivatives of the shape function More...

const VariableValue &	_u_vel

const VariableValue &	_v_vel

const VariableValue &	_w_vel

const VariableValue &	_p

const bool	_picard

const VariableValue *const	_u_vel_previous_nl

const VariableValue *const	_v_vel_previous_nl

const VariableValue *const	_w_vel_previous_nl

const VariableGradient &	_grad_u_vel

const VariableGradient &	_grad_v_vel

const VariableGradient &	_grad_w_vel

const VariableGradient &	_grad_p

const VariableSecond &	_second_u_vel

const VariableSecond &	_second_v_vel

const VariableSecond &	_second_w_vel

const VariableValue &	_u_vel_dot

const VariableValue &	_v_vel_dot

const VariableValue &	_w_vel_dot

const VariableValue &	_d_u_vel_dot_du

const VariableValue &	_d_v_vel_dot_dv

const VariableValue &	_d_w_vel_dot_dw

unsigned	_u_vel_var_number

unsigned	_v_vel_var_number

unsigned	_w_vel_var_number

unsigned	_p_var_number

RealVectorValue	_gravity

const MaterialProperty< Real > &	_mu

const MaterialProperty< Real > &	_rho

const Real &	_alpha

bool	_laplace

bool	_convective_term

bool	_transient_term

const bool	_disps_provided
	Whether displacements are provided. More...

const VariableValue &	_disp_x_dot
	Time derivative of the x-displacement, mesh velocity in the x-direction. More...

const VariableValue &	_disp_y_dot
	Time derivative of the y-displacement, mesh velocity in the y-direction. More...

const VariableValue &	_disp_z_dot
	Time derivative of the z-displacement, mesh velocity in the z-direction. More...

const unsigned int	_rz_radial_coord
	The radial coordinate index for RZ coordinate systems. More...

MooseVariable &	_var

const VariableTestValue &	_test

const VariableTestGradient &	_grad_test

const VariablePhiValue &	_phi

const VariablePhiGradient &	_grad_phi

const VariableValue &	_u

const VariableGradient &	_grad_u

const Elem *const &	_current_elem

const Real &	_current_elem_volume

unsigned int	_qp

const MooseArray< Point > &	_q_point

const QBase *const &	_qrule

const MooseArray< Real > &	_JxW

const MooseArray< Real > &	_coord

unsigned int	_i

unsigned int	_j

bool	_has_save_in

std::vector< MooseVariableFEBase *>	_save_in

std::vector< AuxVariableName >	_save_in_strings

bool	_has_diag_save_in

std::vector< MooseVariableFEBase *>	_diag_save_in

std::vector< AuxVariableName >	_diag_save_in_strings

std::vector< unsigned int >	_displacements

const bool	_use_displaced_mesh

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

THREAD_ID	_tid

Assembly &	_assembly

MooseMesh &	_mesh

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

DenseVector< Number >	_local_re

DenseMatrix< Number >	_local_ke

DenseMatrix< Number >	_nonlocal_ke

const MaterialData *	_blk_material_data

const InputParameters &	_c_parameters

const std::string &	_c_name

const std::string &	_c_type

FEProblemBase &	_c_fe_problem

const SystemBase *const	_c_sys

std::unordered_map< std::string, std::vector< MooseVariableFieldBase *> >	_coupled_vars

std::vector< MooseVariableFieldBase *>	_coupled_moose_vars

std::vector< MooseVariable *>	_coupled_standard_moose_vars

std::vector< VectorMooseVariable *>	_coupled_vector_moose_vars

std::vector< ArrayMooseVariable *>	_coupled_array_moose_vars

std::vector< MooseVariableFV< Real > *>	_coupled_standard_fv_moose_vars

std::vector< MooseLinearVariableFV< Real > *>	_coupled_standard_linear_fv_moose_vars

const std::unordered_map< std::string, std::string > &	_new_to_deprecated_coupled_vars

bool	_c_nodal

bool	_c_is_implicit

const bool	_c_allow_element_to_nodal_coupling

THREAD_ID	_c_tid

std::unordered_map< std::string, std::vector< std::unique_ptr< VariableValue > > >	_default_value

std::unordered_map< std::string, std::unique_ptr< MooseArray< DualReal > > >	_ad_default_value

std::unordered_map< std::string, std::unique_ptr< VectorVariableValue > >	_default_vector_value

std::unordered_map< std::string, std::unique_ptr< ArrayVariableValue > >	_default_array_value

std::unordered_map< std::string, std::unique_ptr< MooseArray< ADRealVectorValue > > >	_ad_default_vector_value

VariableValue	_default_value_zero

VariableGradient	_default_gradient

MooseArray< ADRealVectorValue >	_ad_default_gradient

MooseArray< ADRealTensorValue >	_ad_default_vector_gradient

VariableSecond	_default_second

MooseArray< ADRealTensorValue >	_ad_default_second

const VariableValue &	_zero

const VariablePhiValue &	_phi_zero

const MooseArray< DualReal > &	_ad_zero

const VariableGradient &	_grad_zero

const MooseArray< ADRealVectorValue > &	_ad_grad_zero

const VariablePhiGradient &	_grad_phi_zero

const VariableSecond &	_second_zero

const MooseArray< ADRealTensorValue > &	_ad_second_zero

const VariablePhiSecond &	_second_phi_zero

const VectorVariableValue &	_vector_zero

const VectorVariableCurl &	_vector_curl_zero

VectorVariableValue	_default_vector_value_zero

VectorVariableGradient	_default_vector_gradient

VectorVariableCurl	_default_vector_curl

ArrayVariableValue	_default_array_value_zero

ArrayVariableGradient	_default_array_gradient

ArrayVariableCurl	_default_array_curl

bool	_coupleable_neighbor

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

const InputParameters &	_mi_params

const std::string	_mi_name

const MooseObjectName	_mi_moose_object_name

FEProblemBase &	_mi_feproblem

SubProblem &	_mi_subproblem

const THREAD_ID	_mi_tid

const Moose::MaterialDataType	_material_data_type

MaterialData &	_material_data

bool	_stateful_allowed

bool	_get_material_property_called

std::vector< std::unique_ptr< PropertyValue > >	_default_properties

std::unordered_set< unsigned int >	_material_property_dependencies

const MaterialPropertyName	_get_suffix

const bool	_use_interpolated_state

GeometricSearchData &	_geometric_search_data

bool	_requires_geometric_search

bool	_nodal

MooseVariableFE< Real > *	_variable

MooseVariableFV< Real > *	_fv_variable

MooseLinearVariableFV< Real > *	_linear_fv_variable

MooseVariableField< Real > *	_field_variable

Assembly *	_mvi_assembly

const Parallel::Communicator &	_communicator

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Detailed Description

This class computes strong and weak components of the INS governing equations.

These terms can then be assembled in child classes

Definition at line 18 of file INSBase.h.

Constructor & Destructor Documentation

◆ INSBase()

INSBase::INSBase ( const InputParameters & parameters )

Definition at line 52 of file INSBase.C.

   : Kernel(parameters),
     _second_phi(_assembly.secondPhi()),
 
     // Coupled variables
     _u_vel(coupledValue("u")),
     _v_vel(coupledValue("v")),
     _w_vel(coupledValue("w")),
     _p(coupledValue(NS::pressure)),
 
     _picard(getParam<bool>("picard")),
     _u_vel_previous_nl(_picard ? &coupledValuePreviousNL("u") : nullptr),
     _v_vel_previous_nl(_picard ? &coupledValuePreviousNL("v") : nullptr),
     _w_vel_previous_nl(_picard ? &coupledValuePreviousNL("w") : nullptr),
 
     // Gradients
     _grad_u_vel(coupledGradient("u")),
     _grad_v_vel(coupledGradient("v")),
     _grad_w_vel(coupledGradient("w")),
     _grad_p(coupledGradient(NS::pressure)),
 
     // second derivative tensors
     _second_u_vel(coupledSecond("u")),
     _second_v_vel(coupledSecond("v")),
     _second_w_vel(coupledSecond("w")),
 
     // time derivatives
     _u_vel_dot(_is_transient ? coupledDot("u") : _zero),
     _v_vel_dot(_is_transient ? coupledDot("v") : _zero),
     _w_vel_dot(_is_transient ? coupledDot("w") : _zero),
 
     // derivatives of time derivatives
     _d_u_vel_dot_du(_is_transient ? coupledDotDu("u") : _zero),
     _d_v_vel_dot_dv(_is_transient ? coupledDotDu("v") : _zero),
     _d_w_vel_dot_dw(_is_transient ? coupledDotDu("w") : _zero),
 
     // Variable numberings
     _u_vel_var_number(coupled("u")),
     _v_vel_var_number(coupled("v")),
     _w_vel_var_number(coupled("w")),
     _p_var_number(coupled(NS::pressure)),
 
     _gravity(getParam<RealVectorValue>("gravity")),
 
     // Material properties
     _mu(getMaterialProperty<Real>("mu_name")),
     _rho(getMaterialProperty<Real>("rho_name")),
 
     _alpha(getParam<Real>("alpha")),
     _laplace(getParam<bool>("laplace")),
     _convective_term(getParam<bool>("convective_term")),
     _transient_term(getParam<bool>("transient_term")),
 
     // Displacements for mesh velocity for ALE simulations
     _disps_provided(isParamValid("disp_x")),
     _disp_x_dot(isParamValid("disp_x") ? coupledDot("disp_x") : _zero),
     _disp_y_dot(isParamValid("disp_y") ? coupledDot("disp_y") : _zero),
     _disp_z_dot(isParamValid("disp_z") ? coupledDot("disp_z") : _zero),
 
     _rz_radial_coord(_mesh.getAxisymmetricRadialCoord())
 {
   if (_picard && _disps_provided)
     paramError("picard",
                "Picard is not currently supported for ALE-type simulations in which we subtract "
                "the mesh velocity from the velocity variables");
 }

◆ ~INSBase()

virtual INSBase::~INSBase ( )

inlinevirtual

Definition at line 25 of file INSBase.h.

25 {}

Member Function Documentation

◆ computeQpJacobian()

virtual Real INSBase::computeQpJacobian ( )

protectedpure virtual

Reimplemented from Kernel.

Implemented in INSMomentumLaplaceFormRZ, INSMomentumTractionFormRZ, INSMass, INSMomentumBase, AdvectionSUPG, ConvectedMesh, and ConvectedMeshPSPG.

◆ computeQpOffDiagJacobian()

virtual Real INSBase::computeQpOffDiagJacobian ( unsigned jvar )

protectedpure virtual

Implemented in INSMassRZ, INSMomentumLaplaceFormRZ, INSMomentumTractionFormRZ, INSMass, INSMomentumBase, and AdvectionSUPG.

◆ computeQpResidual()

virtual Real INSBase::computeQpResidual ( )

protectedpure virtual

Implements Kernel.

Implemented in INSMassRZ, INSMomentumLaplaceFormRZ, INSMomentumTractionFormRZ, INSMass, INSMomentumBase, AdvectionSUPG, ConvectedMesh, and ConvectedMeshPSPG.

◆ convectiveTerm()

RealVectorValue INSBase::convectiveTerm ( )

protectedvirtual

Definition at line 132 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

 {
   const auto U = _picard ? RealVectorValue((*_u_vel_previous_nl)[_qp],
                                            (*_v_vel_previous_nl)[_qp],
                                            (*_w_vel_previous_nl)[_qp])
                          : RealVectorValue(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   return _rho[_qp] *
          RealVectorValue(U * _grad_u_vel[_qp], U * _grad_v_vel[_qp], U * _grad_w_vel[_qp]);
 }

◆ dConvecDUComp()

RealVectorValue INSBase::dConvecDUComp ( unsigned comp )

protectedvirtual

Definition at line 143 of file INSBase.C.

Referenced by INSMomentumBase::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), and INSMass::computeQpPGOffDiagJacobian().

 {
   if (_picard)
   {
     RealVectorValue U(
         (*_u_vel_previous_nl)[_qp], (*_v_vel_previous_nl)[_qp], (*_w_vel_previous_nl)[_qp]);
     RealVectorValue convective_term;
     convective_term(comp) = _rho[_qp] * U * _grad_phi[_j][_qp];
 
     return convective_term;
   }
   else
   {
     RealVectorValue U(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
     RealVectorValue d_U_d_comp(0, 0, 0);
     d_U_d_comp(comp) = _phi[_j][_qp];
 
     RealVectorValue convective_term = _rho[_qp] * RealVectorValue(d_U_d_comp * _grad_u_vel[_qp],
                                                                   d_U_d_comp * _grad_v_vel[_qp],
                                                                   d_U_d_comp * _grad_w_vel[_qp]);
     convective_term(comp) += _rho[_qp] * U * _grad_phi[_j][_qp];
 
     return convective_term;
   }
 }

◆ dStrongPressureDPressure()

RealVectorValue INSBase::dStrongPressureDPressure ( )

protectedvirtual

Definition at line 277 of file INSBase.C.

Referenced by INSMomentumBase::computeQpOffDiagJacobian(), and INSMass::computeQpPGJacobian().

 {
   return _grad_phi[_j][_qp];
 }

◆ dStrongViscDUCompLaplace()

RealVectorValue INSBase::dStrongViscDUCompLaplace ( unsigned comp )

protectedvirtual

Reimplemented in INSMassRZ, and INSMomentumLaplaceFormRZ.

Definition at line 185 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMomentumLaplaceFormRZ::dStrongViscDUCompLaplace(), and INSMassRZ::dStrongViscDUCompLaplace().

 {
   RealVectorValue viscous_term(0, 0, 0);
   viscous_term(comp) = -_mu[_qp] * _second_phi[_j][_qp].tr();
 
   return viscous_term;
 }

◆ dStrongViscDUCompLaplaceRZ()

RealVectorValue INSBase::dStrongViscDUCompLaplaceRZ ( const unsigned int comp ) const

protected

Computes the Jacobian for the additional RZ terms for the Laplace form of the strong viscous term for the given velocity component comp.

Definition at line 388 of file INSBase.C.

Referenced by INSMomentumLaplaceFormRZ::dStrongViscDUCompLaplace(), INSMassRZ::dStrongViscDUCompLaplace(), and dStrongViscDUCompTractionRZ().

 {
   const auto r = _q_point[_qp](_rz_radial_coord);
   RealVectorValue add_jac;
   add_jac(comp) = -_mu[_qp] * _grad_phi[_j][_qp](_rz_radial_coord) / r;
   if (comp == _rz_radial_coord)
     add_jac(comp) += _mu[_qp] * _phi[_j][_qp] / (r * r);
 
   return add_jac;
 }

◆ dStrongViscDUCompTraction()

RealVectorValue INSBase::dStrongViscDUCompTraction ( unsigned comp )

protectedvirtual

Reimplemented in INSMassRZ, and INSMomentumTractionFormRZ.

Definition at line 194 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMomentumTractionFormRZ::dStrongViscDUCompTraction(), and INSMassRZ::dStrongViscDUCompTraction().

 {
   RealVectorValue viscous_term(0, 0, 0);
   viscous_term(comp) = -_mu[_qp] * (_second_phi[_j][_qp](0, 0) + _second_phi[_j][_qp](1, 1) +
                                     _second_phi[_j][_qp](2, 2));
   for (unsigned i = 0; i < 3; i++)
     viscous_term(i) += -_mu[_qp] * _second_phi[_j][_qp](i, comp);
 
   return viscous_term;
 }

◆ dStrongViscDUCompTractionRZ()

RealVectorValue INSBase::dStrongViscDUCompTractionRZ ( const unsigned int comp ) const

protected

Computes the Jacobian for the additional RZ terms for the Traction form of the strong viscous term for the given velocity component comp.

Definition at line 415 of file INSBase.C.

Referenced by INSMomentumTractionFormRZ::dStrongViscDUCompTraction(), and INSMassRZ::dStrongViscDUCompTraction().

 {
   auto ret = dStrongViscDUCompLaplaceRZ(comp);
   if (comp != _rz_radial_coord)
     return ret;
 
   const auto r = _q_point[_qp](_rz_radial_coord);
 
   ret += -_mu[_qp] * _grad_phi[_j][_qp] / r;
   ret(_rz_radial_coord) += _mu[_qp] * _phi[_j][_qp] / (r * r);
 
   return ret;
 }

◆ dTauDUComp()

Real INSBase::dTauDUComp ( unsigned comp )

protectedvirtual

Definition at line 349 of file INSBase.C.

Referenced by ConvectedMesh::computePGVelocityJacobian(), ConvectedMeshPSPG::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), and INSMass::computeQpPGOffDiagJacobian().

 {
   Real nu = _mu[_qp] / _rho[_qp];
   const auto U = relativeVelocity();
   Real h = _current_elem->hmax();
   Real transient_part = _transient_term ? 4. / (_dt * _dt) : 0.;
   return -_alpha / 2. *
          std::pow(transient_part + (2. * U.norm() / h) * (2. * U.norm() / h) +
                       9. * (4. * nu / (h * h)) * (4. * nu / (h * h)),
                   -1.5) *
          2. * (2. * U.norm() / h) * 2. / h * U(comp) * _phi[_j][_qp] /
          (U.norm() + std::numeric_limits<double>::epsilon());
 }

◆ dTimeDerivativeDUComp()

RealVectorValue INSBase::dTimeDerivativeDUComp ( unsigned comp )

protectedvirtual

Definition at line 301 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), and INSMass::computeQpPGOffDiagJacobian().

 {
   Real base = _rho[_qp] * _phi[_j][_qp];
   switch (comp)
   {
     case 0:
       return RealVectorValue(base * _d_u_vel_dot_du[_qp], 0, 0);
 
     case 1:
       return RealVectorValue(0, base * _d_v_vel_dot_dv[_qp], 0);
 
     case 2:
       return RealVectorValue(0, 0, base * _d_w_vel_dot_dw[_qp]);
 
     default:
       mooseError("comp must be 0, 1, or 2");
   }
 }

◆ dWeakPressureDPressure()

Real INSBase::dWeakPressureDPressure ( )

protectedvirtual

Definition at line 283 of file INSBase.C.

Referenced by INSMomentumBase::computeQpOffDiagJacobian().

 {
   return -_phi[_j][_qp];
 }

◆ dWeakViscDUCompLaplace()

RealVectorValue INSBase::dWeakViscDUCompLaplace ( )

protectedvirtual

Definition at line 253 of file INSBase.C.

 {
   return _mu[_qp] * _grad_phi[_j][_qp];
 }

◆ dWeakViscDUCompTraction()

RealVectorValue INSBase::dWeakViscDUCompTraction ( )

protectedvirtual

Definition at line 259 of file INSBase.C.

 {
   return _mu[_qp] * _grad_phi[_j][_qp];
 }

◆ gravityTerm()

RealVectorValue INSBase::gravityTerm ( )

protectedvirtual

Definition at line 289 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

 {
   return -_rho[_qp] * _gravity;
 }

◆ relativeVelocity()

RealVectorValue INSBase::relativeVelocity ( ) const

protected

Compute the velocity.

If displacements are provided, then this routine will subtract the time derivative of the displacements, e.g. the mesh velocity

Definition at line 120 of file INSBase.C.

Referenced by ConvectedMesh::computePGVelocityJacobian(), ConvectedMesh::computeQpJacobian(), ConvectedMesh::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMomentumBase::computeQpPGResidual(), ConvectedMesh::computeQpResidual(), dTauDUComp(), tau(), and tauNodal().

 {
   auto U = _picard ? RealVectorValue((*_u_vel_previous_nl)[_qp],
                                      (*_v_vel_previous_nl)[_qp],
                                      (*_w_vel_previous_nl)[_qp])
                    : RealVectorValue(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
   if (_disps_provided)
     U -= RealVectorValue{_disp_x_dot[_qp], _disp_y_dot[_qp], _disp_z_dot[_qp]};
   return U;
 }

◆ strongPressureTerm()

RealVectorValue INSBase::strongPressureTerm ( )

protectedvirtual

Definition at line 265 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), and INSMomentumBase::computeQpResidual().

 {
   return _grad_p[_qp];
 }

◆ strongViscousTermLaplace()

RealVectorValue INSBase::strongViscousTermLaplace ( )

protectedvirtual

Reimplemented in INSMassRZ, and INSMomentumLaplaceFormRZ.

Definition at line 170 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), and strongViscousTermTraction().

 {
   return -_mu[_qp] *
          RealVectorValue(_second_u_vel[_qp].tr(), _second_v_vel[_qp].tr(), _second_w_vel[_qp].tr());
 }

◆ strongViscousTermLaplaceRZ()

RealVectorValue INSBase::strongViscousTermLaplaceRZ ( ) const

protected

Computes the additional RZ terms for the Laplace form of the strong viscous term.

Definition at line 364 of file INSBase.C.

Referenced by INSMomentumLaplaceFormRZ::strongViscousTermLaplace(), INSMassRZ::strongViscousTermLaplace(), and strongViscousTermTractionRZ().

 {
   // To understand the code below, visit
   // https://en.wikipedia.org/wiki/Del_in_cylindrical_and_spherical_coordinates.
   // The u_r / r^2 term comes from the vector Laplacian. The -du_r/dr * 1/r term comes from
   // the scalar Laplacian. The scalar Laplacian in axisymmetric cylindrical coordinates is
   // equivalent to the Cartesian Laplacian plus a 1/r * df/dr term. And of course we are
   // applying a minus sign here because the strong form is -\nabala^2 * \vec{u}
 
   const auto r = _q_point[_qp](_rz_radial_coord);
   RealVectorValue rz_term;
   rz_term(0) = -_mu[_qp] * _grad_u_vel[_qp](_rz_radial_coord) / r;
   rz_term(1) = -_mu[_qp] * _grad_v_vel[_qp](_rz_radial_coord) / r;
   mooseAssert((_rz_radial_coord == 0) || (_rz_radial_coord == 1),
               "We expect X or Y as the possible radial coordinate");
   if (_rz_radial_coord == 0)
     rz_term(0) += _mu[_qp] * _u_vel[_qp] / (r * r);
   else
     rz_term(1) += _mu[_qp] * _v_vel[_qp] / (r * r);
 
   return rz_term;
 }

◆ strongViscousTermTraction()

RealVectorValue INSBase::strongViscousTermTraction ( )

protectedvirtual

Reimplemented in INSMassRZ, and INSMomentumTractionFormRZ.

Definition at line 177 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSMomentumTractionFormRZ::strongViscousTermTraction(), and INSMassRZ::strongViscousTermTraction().

 {
   return INSBase::strongViscousTermLaplace() -
          _mu[_qp] *
              (_second_u_vel[_qp].row(0) + _second_v_vel[_qp].row(1) + _second_w_vel[_qp].row(2));
 }

◆ strongViscousTermTractionRZ()

RealVectorValue INSBase::strongViscousTermTractionRZ ( ) const

protected

Computes the additional RZ terms for the Traction form of the strong viscous term.

Definition at line 400 of file INSBase.C.

Referenced by INSMomentumTractionFormRZ::strongViscousTermTraction(), and INSMassRZ::strongViscousTermTraction().

 {
   auto ret = strongViscousTermLaplaceRZ();
 
   const auto r = _q_point[_qp](_rz_radial_coord);
 
   const auto & grad_r_vel = (_rz_radial_coord == 0) ? _grad_u_vel[_qp] : _grad_v_vel[_qp];
   const auto & r_vel = (_rz_radial_coord == 0) ? _u_vel[_qp] : _v_vel[_qp];
   ret += -_mu[_qp] * grad_r_vel / r;
   ret(_rz_radial_coord) += _mu[_qp] * r_vel / (r * r);
 
   return ret;
 }

◆ tau()

Real INSBase::tau ( )

protectedvirtual

Definition at line 321 of file INSBase.C.

Referenced by ConvectedMesh::computePGVelocityJacobian(), ConvectedMesh::computeQpJacobian(), AdvectionSUPG::computeQpJacobian(), ConvectedMesh::computeQpOffDiagJacobian(), ConvectedMeshPSPG::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMass::computeQpPGJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), ConvectedMeshPSPG::computeQpResidual(), ConvectedMesh::computeQpResidual(), and AdvectionSUPG::computeQpResidual().

 {
   Real nu = _mu[_qp] / _rho[_qp];
   const auto U = relativeVelocity();
   Real h = _current_elem->hmax();
   Real transient_part = _transient_term ? 4. / (_dt * _dt) : 0.;
   return _alpha / std::sqrt(transient_part + (2. * U.norm() / h) * (2. * U.norm() / h) +
                             9. * (4. * nu / (h * h)) * (4. * nu / (h * h)));
 }

◆ tauNodal()

Real INSBase::tauNodal ( )

protectedvirtual

Provides tau which yields superconvergence for 1D advection-diffusion.

Definition at line 332 of file INSBase.C.

Referenced by AdvectionSUPG::computeQpJacobian(), and AdvectionSUPG::computeQpResidual().

 {
   Real nu = _mu[_qp] / _rho[_qp];
   const auto U = relativeVelocity();
   Real h = _current_elem->hmax();
   Real xi;
   if (nu < std::numeric_limits<Real>::epsilon())
     xi = 1;
   else
   {
     Real alpha = U.norm() * h / (2 * nu);
     xi = 1. / std::tanh(alpha) - 1. / alpha;
   }
   return h / (2 * U.norm()) * xi;
 }

◆ timeDerivativeTerm()

RealVectorValue INSBase::timeDerivativeTerm ( )

protectedvirtual

Definition at line 295 of file INSBase.C.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), and INSMomentumBase::computeQpPGResidual().

 {
   return _rho[_qp] * RealVectorValue(_u_vel_dot[_qp], _v_vel_dot[_qp], _w_vel_dot[_qp]);
 }

◆ validParams()

InputParameters INSBase::validParams ( )

static

Definition at line 15 of file INSBase.C.

Referenced by ConvectedMesh::validParams(), ConvectedMeshPSPG::validParams(), AdvectionSUPG::validParams(), INSMass::validParams(), and INSMomentumBase::validParams().

 {
   InputParameters params = Kernel::validParams();
 
   params.addClassDescription("This class computes various strong and weak components of the "
                              "incompressible navier stokes equations which can then be assembled "
                              "together in child classes.");
   // Coupled variables
   params.addRequiredCoupledVar("u", "x-velocity");
   params.addCoupledVar("v", 0, "y-velocity"); // only required in 2D and 3D
   params.addCoupledVar("w", 0, "z-velocity"); // only required in 3D
   params.addRequiredCoupledVar(NS::pressure, "pressure");
 
   params.addParam<RealVectorValue>(
       "gravity", RealVectorValue(0, 0, 0), "Direction of the gravity vector");
 
   params.addParam<MaterialPropertyName>("mu_name", "mu", "The name of the dynamic viscosity");
   params.addParam<MaterialPropertyName>("rho_name", "rho", "The name of the density");
 
   params.addParam<Real>("alpha", 1., "Multiplicative factor on the stabilization parameter tau.");
   params.addParam<bool>(
       "laplace", true, "Whether the viscous term of the momentum equations is in laplace form.");
   params.addParam<bool>("convective_term", true, "Whether to include the convective term.");
   params.addParam<bool>("transient_term",
                         false,
                         "Whether there should be a transient term in the momentum residuals.");
   params.addCoupledVar("disp_x", "The x displacement");
   params.addCoupledVar("disp_y", "The y displacement");
   params.addCoupledVar("disp_z", "The z displacement");
   params.addParam<bool>("picard",
                         false,
                         "Whether we are applying a Picard strategy in which case we will linearize "
                         "the nonlinear convective term.");
 
   return params;
 }

◆ weakPressureTerm()

Real INSBase::weakPressureTerm ( )

protectedvirtual

Definition at line 271 of file INSBase.C.

Referenced by INSMomentumBase::computeQpResidual().

 {
   return -_p[_qp];
 }

◆ weakViscousTermLaplace()

RealVectorValue INSBase::weakViscousTermLaplace ( unsigned comp )

protectedvirtual

Definition at line 206 of file INSBase.C.

 {
   switch (comp)
   {
     case 0:
       return _mu[_qp] * _grad_u_vel[_qp];
 
     case 1:
       return _mu[_qp] * _grad_v_vel[_qp];
 
     case 2:
       return _mu[_qp] * _grad_w_vel[_qp];
 
     default:
       return _zero[_qp];
   }
 }

◆ weakViscousTermTraction()

RealVectorValue INSBase::weakViscousTermTraction ( unsigned comp )

protectedvirtual

Definition at line 225 of file INSBase.C.

 {
   switch (comp)
   {
     case 0:
     {
       RealVectorValue transpose(_grad_u_vel[_qp](0), _grad_v_vel[_qp](0), _grad_w_vel[_qp](0));
       return _mu[_qp] * _grad_u_vel[_qp] + _mu[_qp] * transpose;
     }
 
     case 1:
     {
       RealVectorValue transpose(_grad_u_vel[_qp](1), _grad_v_vel[_qp](1), _grad_w_vel[_qp](1));
       return _mu[_qp] * _grad_v_vel[_qp] + _mu[_qp] * transpose;
     }
 
     case 2:
     {
       RealVectorValue transpose(_grad_u_vel[_qp](2), _grad_v_vel[_qp](2), _grad_w_vel[_qp](2));
       return _mu[_qp] * _grad_w_vel[_qp] + _mu[_qp] * transpose;
     }
 
     default:
       return _zero[_qp];
   }
 }

Member Data Documentation

◆ _alpha

const Real& INSBase::_alpha

protected

Definition at line 136 of file INSBase.h.

Referenced by dTauDUComp(), and tau().

◆ _convective_term

bool INSBase::_convective_term

protected

Definition at line 138 of file INSBase.h.

Referenced by INSMomentumBase::computeQpJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), INSMomentumBase::computeQpResidual(), and INSMomentumBase::INSMomentumBase().

◆ _d_u_vel_dot_du

const VariableValue& INSBase::_d_u_vel_dot_du

protected

Definition at line 120 of file INSBase.h.

Referenced by dTimeDerivativeDUComp().

◆ _d_v_vel_dot_dv

const VariableValue& INSBase::_d_v_vel_dot_dv

protected

Definition at line 121 of file INSBase.h.

Referenced by dTimeDerivativeDUComp().

◆ _d_w_vel_dot_dw

const VariableValue& INSBase::_d_w_vel_dot_dw

protected

Definition at line 122 of file INSBase.h.

Referenced by dTimeDerivativeDUComp().

◆ _disp_x_dot

const VariableValue& INSBase::_disp_x_dot

protected

Time derivative of the x-displacement, mesh velocity in the x-direction.

Definition at line 144 of file INSBase.h.

Referenced by relativeVelocity().

◆ _disp_y_dot

const VariableValue& INSBase::_disp_y_dot

protected

Time derivative of the y-displacement, mesh velocity in the y-direction.

Definition at line 146 of file INSBase.h.

Referenced by relativeVelocity().

◆ _disp_z_dot

const VariableValue& INSBase::_disp_z_dot

protected

Time derivative of the z-displacement, mesh velocity in the z-direction.

Definition at line 148 of file INSBase.h.

Referenced by relativeVelocity().

◆ _disps_provided

const bool INSBase::_disps_provided

protected

Whether displacements are provided.

Definition at line 142 of file INSBase.h.

Referenced by INSBase(), and relativeVelocity().

◆ _grad_p

const VariableGradient& INSBase::_grad_p

protected

Definition at line 107 of file INSBase.h.

Referenced by strongPressureTerm().

◆ _grad_u_vel

const VariableGradient& INSBase::_grad_u_vel

protected

Definition at line 104 of file INSBase.h.

Referenced by INSMass::computeQpResidual(), INSMomentumTractionForm::computeQpResidualViscousPart(), convectiveTerm(), dConvecDUComp(), ConvectedMeshPSPG::dStrongResidualDDisp(), ConvectedMeshPSPG::strongResidual(), strongViscousTermLaplaceRZ(), strongViscousTermTractionRZ(), weakViscousTermLaplace(), and weakViscousTermTraction().

◆ _grad_v_vel

const VariableGradient& INSBase::_grad_v_vel

protected

Definition at line 105 of file INSBase.h.

Referenced by INSMass::computeQpResidual(), INSMomentumTractionForm::computeQpResidualViscousPart(), convectiveTerm(), dConvecDUComp(), ConvectedMeshPSPG::dStrongResidualDDisp(), ConvectedMeshPSPG::strongResidual(), strongViscousTermLaplaceRZ(), strongViscousTermTractionRZ(), weakViscousTermLaplace(), and weakViscousTermTraction().

◆ _grad_w_vel

const VariableGradient& INSBase::_grad_w_vel

protected

Definition at line 106 of file INSBase.h.

Referenced by INSMass::computeQpResidual(), INSMomentumTractionForm::computeQpResidualViscousPart(), convectiveTerm(), dConvecDUComp(), ConvectedMeshPSPG::dStrongResidualDDisp(), ConvectedMeshPSPG::strongResidual(), weakViscousTermLaplace(), and weakViscousTermTraction().

◆ _gravity

RealVectorValue INSBase::_gravity

protected

Definition at line 130 of file INSBase.h.

Referenced by gravityTerm().

◆ _laplace

bool INSBase::_laplace

protected

Definition at line 137 of file INSBase.h.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), and INSMomentumBase::computeQpPGResidual().

◆ _mu

const MaterialProperty<Real>& INSBase::_mu

protected

Definition at line 133 of file INSBase.h.

Referenced by INSMomentumTractionFormRZ::computeQpJacobian(), INSMomentumLaplaceFormRZ::computeQpJacobian(), INSMomentumTractionForm::computeQpJacobianViscousPart(), INSMomentumLaplaceForm::computeQpJacobianViscousPart(), INSMomentumTractionForm::computeQpOffDiagJacobianViscousPart(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), INSMomentumLaplaceForm::computeQpResidualViscousPart(), INSMomentumTractionForm::computeQpResidualViscousPart(), dStrongViscDUCompLaplace(), dStrongViscDUCompLaplaceRZ(), dStrongViscDUCompTraction(), dStrongViscDUCompTractionRZ(), dTauDUComp(), dWeakViscDUCompLaplace(), dWeakViscDUCompTraction(), strongViscousTermLaplace(), strongViscousTermLaplaceRZ(), strongViscousTermTraction(), strongViscousTermTractionRZ(), tau(), tauNodal(), weakViscousTermLaplace(), and weakViscousTermTraction().

◆ _p

const VariableValue& INSBase::_p

protected

Definition at line 96 of file INSBase.h.

Referenced by INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), and weakPressureTerm().

◆ _p_var_number

unsigned INSBase::_p_var_number

protected

Definition at line 128 of file INSBase.h.

Referenced by INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumTractionFormRZ::computeQpOffDiagJacobian(), and INSMomentumLaplaceFormRZ::computeQpOffDiagJacobian().

◆ _picard

const bool INSBase::_picard

protected

Definition at line 98 of file INSBase.h.

Referenced by convectiveTerm(), dConvecDUComp(), INSBase(), and relativeVelocity().

◆ _rho

const MaterialProperty<Real>& INSBase::_rho

protected

Definition at line 134 of file INSBase.h.

Referenced by INSMass::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), convectiveTerm(), dConvecDUComp(), dTauDUComp(), dTimeDerivativeDUComp(), gravityTerm(), tau(), tauNodal(), and timeDerivativeTerm().

◆ _rz_radial_coord

const unsigned int INSBase::_rz_radial_coord

protected

The radial coordinate index for RZ coordinate systems.

Definition at line 151 of file INSBase.h.

Referenced by INSMomentumTractionFormRZ::computeQpJacobian(), INSMomentumLaplaceFormRZ::computeQpJacobian(), INSMomentumTractionFormRZ::computeQpOffDiagJacobian(), INSMomentumLaplaceFormRZ::computeQpOffDiagJacobian(), INSMassRZ::computeQpOffDiagJacobian(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), INSMassRZ::computeQpResidual(), dStrongViscDUCompLaplaceRZ(), dStrongViscDUCompTractionRZ(), strongViscousTermLaplaceRZ(), and strongViscousTermTractionRZ().

◆ _second_phi

const VariablePhiSecond& INSBase::_second_phi

protected

second derivatives of the shape function

Definition at line 90 of file INSBase.h.

Referenced by dStrongViscDUCompLaplace(), and dStrongViscDUCompTraction().

◆ _second_u_vel

const VariableSecond& INSBase::_second_u_vel

protected

Definition at line 110 of file INSBase.h.

Referenced by strongViscousTermLaplace(), and strongViscousTermTraction().

◆ _second_v_vel

const VariableSecond& INSBase::_second_v_vel

protected

Definition at line 111 of file INSBase.h.

Referenced by strongViscousTermLaplace(), and strongViscousTermTraction().

◆ _second_w_vel

const VariableSecond& INSBase::_second_w_vel

protected

Definition at line 112 of file INSBase.h.

Referenced by strongViscousTermLaplace(), and strongViscousTermTraction().

◆ _transient_term

bool INSBase::_transient_term

protected

Definition at line 139 of file INSBase.h.

Referenced by INSMomentumBase::computeQpPGJacobian(), INSMass::computeQpPGOffDiagJacobian(), INSMass::computeQpPGResidual(), INSMomentumBase::computeQpPGResidual(), dTauDUComp(), and tau().

◆ _u_vel

const VariableValue& INSBase::_u_vel

protected

Definition at line 93 of file INSBase.h.

Referenced by AdvectionSUPG::computeQpJacobian(), AdvectionSUPG::computeQpResidual(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), INSMassRZ::computeQpResidual(), convectiveTerm(), dConvecDUComp(), relativeVelocity(), strongViscousTermLaplaceRZ(), and strongViscousTermTractionRZ().

◆ _u_vel_dot

const VariableValue& INSBase::_u_vel_dot

protected

Definition at line 115 of file INSBase.h.

Referenced by timeDerivativeTerm().

◆ _u_vel_previous_nl

const VariableValue* const INSBase::_u_vel_previous_nl

protected

Definition at line 99 of file INSBase.h.

Referenced by convectiveTerm(), dConvecDUComp(), and relativeVelocity().

◆ _u_vel_var_number

unsigned INSBase::_u_vel_var_number

protected

Definition at line 125 of file INSBase.h.

Referenced by ConvectedMesh::computeQpOffDiagJacobian(), ConvectedMeshPSPG::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMass::computeQpOffDiagJacobian(), INSMassRZ::computeQpOffDiagJacobian(), INSMomentumTractionForm::computeQpOffDiagJacobianViscousPart(), and ConvectedMesh::ConvectedMesh().

◆ _v_vel

const VariableValue& INSBase::_v_vel

protected

Definition at line 94 of file INSBase.h.

Referenced by AdvectionSUPG::computeQpJacobian(), AdvectionSUPG::computeQpResidual(), INSMomentumTractionFormRZ::computeQpResidual(), INSMomentumLaplaceFormRZ::computeQpResidual(), INSMassRZ::computeQpResidual(), convectiveTerm(), dConvecDUComp(), relativeVelocity(), strongViscousTermLaplaceRZ(), and strongViscousTermTractionRZ().

◆ _v_vel_dot

const VariableValue& INSBase::_v_vel_dot

protected

Definition at line 116 of file INSBase.h.

Referenced by timeDerivativeTerm().

◆ _v_vel_previous_nl

const VariableValue* const INSBase::_v_vel_previous_nl

protected

Definition at line 100 of file INSBase.h.

Referenced by convectiveTerm(), dConvecDUComp(), and relativeVelocity().

◆ _v_vel_var_number

unsigned INSBase::_v_vel_var_number

protected

Definition at line 126 of file INSBase.h.

Referenced by ConvectedMesh::computeQpOffDiagJacobian(), ConvectedMeshPSPG::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMass::computeQpOffDiagJacobian(), INSMassRZ::computeQpOffDiagJacobian(), INSMomentumTractionForm::computeQpOffDiagJacobianViscousPart(), and ConvectedMesh::ConvectedMesh().

◆ _w_vel

const VariableValue& INSBase::_w_vel

protected

Definition at line 95 of file INSBase.h.

Referenced by AdvectionSUPG::computeQpJacobian(), AdvectionSUPG::computeQpResidual(), convectiveTerm(), dConvecDUComp(), and relativeVelocity().

◆ _w_vel_dot

const VariableValue& INSBase::_w_vel_dot

protected

Definition at line 117 of file INSBase.h.

Referenced by timeDerivativeTerm().

◆ _w_vel_previous_nl

const VariableValue* const INSBase::_w_vel_previous_nl

protected

Definition at line 101 of file INSBase.h.

Referenced by convectiveTerm(), dConvecDUComp(), and relativeVelocity().

◆ _w_vel_var_number

unsigned INSBase::_w_vel_var_number

protected

Definition at line 127 of file INSBase.h.

Referenced by ConvectedMesh::computeQpOffDiagJacobian(), ConvectedMeshPSPG::computeQpOffDiagJacobian(), INSMass::computeQpOffDiagJacobian(), INSMomentumBase::computeQpOffDiagJacobian(), INSMomentumTractionForm::computeQpOffDiagJacobianViscousPart(), and ConvectedMesh::ConvectedMesh().

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

navier_stokes/include/kernels/INSBase.h
navier_stokes/src/kernels/INSBase.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ INSBase()

◆ ~INSBase()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ convectiveTerm()

◆ dConvecDUComp()

◆ dStrongPressureDPressure()

◆ dStrongViscDUCompLaplace()

◆ dStrongViscDUCompLaplaceRZ()

◆ dStrongViscDUCompTraction()

◆ dStrongViscDUCompTractionRZ()

◆ dTauDUComp()

◆ dTimeDerivativeDUComp()

◆ dWeakPressureDPressure()

◆ dWeakViscDUCompLaplace()

◆ dWeakViscDUCompTraction()

◆ gravityTerm()

◆ relativeVelocity()

◆ strongPressureTerm()

◆ strongViscousTermLaplace()

◆ strongViscousTermLaplaceRZ()

◆ strongViscousTermTraction()

◆ strongViscousTermTractionRZ()

◆ tau()

◆ tauNodal()

◆ timeDerivativeTerm()

◆ validParams()

◆ weakPressureTerm()

◆ weakViscousTermLaplace()

◆ weakViscousTermTraction()

Member Data Documentation

◆ _alpha

◆ _convective_term

◆ _d_u_vel_dot_du

◆ _d_v_vel_dot_dv

◆ _d_w_vel_dot_dw

◆ _disp_x_dot

◆ _disp_y_dot

◆ _disp_z_dot

◆ _disps_provided

◆ _grad_p

◆ _grad_u_vel

◆ _grad_v_vel

◆ _grad_w_vel

◆ _gravity

◆ _laplace

◆ _mu

◆ _p

◆ _p_var_number

◆ _picard

◆ _rho

◆ _rz_radial_coord

◆ _second_phi

◆ _second_u_vel

◆ _second_v_vel

◆ _second_w_vel

◆ _transient_term

◆ _u_vel

◆ _u_vel_dot

◆ _u_vel_previous_nl

◆ _u_vel_var_number

◆ _v_vel

◆ _v_vel_dot

◆ _v_vel_previous_nl

◆ _v_vel_var_number

◆ _w_vel

◆ _w_vel_dot

◆ _w_vel_previous_nl