Adds the axisymmetric viscous source term \(\mu u_r / r^2\) that appears in the vector Laplacian of cylindrical coordinates. More...

#include <LinearFVRZViscousSource.h>

Inheritance diagram for LinearFVRZViscousSource:

[legend]

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

typedef DataFileName	DataFileParameterType

Public Member Functions
	LinearFVRZViscousSource (const InputParameters &params)

virtual void	addMatrixContribution () override

virtual void	addRightHandSideContribution () override

virtual void	setCurrentElemInfo (const ElemInfo *elem_info)

void	setCurrentElemVolume (const Real volume)

virtual const MooseLinearVariableFV< Real > &	variable () const override

void	linkTaggedVectorsAndMatrices (const std::set< TagID > &vector_tags, const std::set< TagID > &matrix_tags)

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

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

bool	isKokkosObject () const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () 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	haveParameter (const std::string &name) const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &name) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (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

std::string	messagePrefix (const bool hit_prefix=true) const

std::string	errorPrefix (const std::string &) const

void	mooseError (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseDeprecatedNoTrace (Args &&... args) const

void	mooseInfo (Args &&... args) const

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

virtual void	initialSetup ()

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	subdomainSetup ()

virtual void	customSetup (const ExecFlagType &)

const 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

Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const

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

Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const

const T &	getKokkosFunctionByName (const FunctionName &name) const

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

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

const UserObjectBase &	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

void	setRandomResetFrequency (ExecFlagType exec_flag)

unsigned long	getRandomLong () const

Real	getRandomReal () const

unsigned int	getSeed (std::size_t id)

unsigned int	getMasterSeed () const

bool	isNodal () const

ExecFlagType	getResetOnTime () const

void	setRandomDataPointer (RandomData *random_data)

virtual void	meshChanged ()

void	useVectorTag (const TagName &tag_name, VectorTagsKey)

void	useVectorTag (TagID tag_id, VectorTagsKey)

void	useMatrixTag (const TagName &tag_name, MatrixTagsKey)

void	useMatrixTag (TagID tag_id, MatrixTagsKey)

bool	isVectorTagged ()

bool	isMatrixTagged ()

bool	hasVectorTags () const

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

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

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

unsigned int	numBlocks () const

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

unsigned int	blocksMaxDimension () const

bool	hasBlocks (const SubdomainName &name) const

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

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

bool	hasBlocks (SubdomainID id) const

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

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

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

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

bool	hasBlockMaterialProperty (const std::string &prop_name)

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

virtual bool	blockRestricted () const

virtual void	checkVariable (const MooseVariableFieldBase &variable) const

MooseVariableBase *	mooseVariableBase () const

MooseVariableField< Real > &	mooseVariableField ()

MooseVariableFE< Real > *	mooseVariable () const

MooseVariableFV< Real > *	mooseVariableFV () const

MooseLinearVariableFV< Real > *	mooseLinearVariableFV () 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)

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 Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

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

static void	setRMParamsAdvection (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short conditional_extended_layers)

static void	setRMParamsDiffusion (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short conditional_extended_layers)

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

Protected Member Functions
Real	computeMatrixContribution () override

Real	computeRightHandSideContribution () override

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

void	requestVariableCellGradient (const std::string &variable_name)

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

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

virtual void	addUserObjectDependencyHelper (const UserObjectBase &) 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 Residuals &residuals, const Indices &dof_indices, const std::vector< Real > &scaling_factors)

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

void	addResiduals (Assembly &assembly, const ADResidualsPacket &packet)

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

void	addResidualsAndJacobian (Assembly &assembly, const ADResidualsPacket &packet)

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

void	addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, const std::vector< Real > &scaling_factors)

void	addJacobian (Assembly &assembly, const ADResidualsPacket &packet)

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

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

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

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

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

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

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

void	setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)

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

void	initializeBlockRestrictable (const MooseObject *moose_object)

void	initializeKokkosBlockRestrictable ()

Moose::CoordinateSystemType	getBlockCoordSystem ()

KOKKOS_FUNCTION dof_id_type	numKokkosBlockElements () const

KOKKOS_FUNCTION dof_id_type	numKokkosBlockNodes () const

KOKKOS_FUNCTION dof_id_type	numKokkosBlockSides () const

KOKKOS_FUNCTION ContiguousElementID	kokkosBlockElementID (Moose::Kokkos::ThreadID tid) const

KOKKOS_FUNCTION ContiguousElementID	kokkosBlockNodeID (Moose::Kokkos::ThreadID tid) const

KOKKOS_FUNCTION auto	kokkosBlockElementSideID (Moose::Kokkos::ThreadID tid) 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)

const Moose::Functor< T > &	getFunctorByName (const std::string &name)

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

Protected Attributes
const ElemInfo *	_current_elem_info

Real	_current_elem_volume

dof_id_type	_dof_id

MooseLinearVariableFV< Real > &	_var

const unsigned int	_var_num

const unsigned int	_sys_num

FEProblemBase &	_fe_problem

SystemBase &	_sys

libMesh::LinearImplicitSystem &	_linear_system

const THREAD_ID	_tid

MooseMesh &	_mesh

std::vector< NumericVector< Number > *>	_vectors

std::vector< SparseMatrix< Number > *>	_matrices

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

SubProblem &	_subproblem

DenseVector< Number >	_local_re

DenseMatrix< Number >	_local_ke

DenseMatrix< Number >	_nonlocal_ke

const MaterialData *	_blk_material_data

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

Private Member Functions
const MooseLinearVariableFVReal &	velocityVar (unsigned int dir) const
	Helper to access the velocity variable for a given direction. More...

Private Attributes
const Moose::Functor< Real > &	_mu
	Dynamic viscosity functor evaluated at each element. More...

const unsigned int	_component
	Momentum component this source acts on (should equal the radial direction) More...

const unsigned int	_rz_radial_coord
	Index of the radial coordinate for the current mesh (0 -> x, 1 -> y, ...) More...

const unsigned int	_dim
	Spatial dimension of the mesh. More...

const bool	_use_deviatoric_terms
	Whether the deviatoric correction (-2/3 div u) is requested. More...

const Moose::CoordinateSystemType	_coord_type
	Coordinate system of the active blocks (must be COORD_RZ) More...

const Real	_stress_multiplier
	Precomputed factor (1 or 2) multiplying the implicit hoop term. More...

std::array< MooseLinearVariableFVReal *, 2 >	_velocity_vars
	Cached pointers to the velocity components required to build divergence. More...

Detailed Description

Adds the axisymmetric viscous source term \(\mu u_r / r^2\) that appears in the vector Laplacian of cylindrical coordinates.

The contribution is only active when the block uses an RZ coordinate system and the momentum component matches the radial direction.

Definition at line 21 of file LinearFVRZViscousSource.h.

Constructor & Destructor Documentation

◆ LinearFVRZViscousSource()

LinearFVRZViscousSource::LinearFVRZViscousSource ( const InputParameters & params )

Definition at line 38 of file LinearFVRZViscousSource.C.

   : LinearFVElementalKernel(params),
     _mu(getFunctor<Real>(NS::mu)),
     _component(getParam<MooseEnum>("momentum_component")),
     _rz_radial_coord(_subproblem.getAxisymmetricRadialCoord()),
     _dim(_subproblem.mesh().dimension()),
     _use_deviatoric_terms(getParam<bool>("use_deviatoric_terms")),
     _coord_type(getBlockCoordSystem()),
     _stress_multiplier(_use_deviatoric_terms ? 2.0 : 1.0),
     _velocity_vars{nullptr, nullptr}
 {
   if (_coord_type != Moose::CoordinateSystemType::COORD_RZ)
     paramError("block", "LinearFVRZViscousSource is only valid on RZ coordinate systems.");
 
   if (_component != _rz_radial_coord)
     paramError("momentum_component", "LinearFVRZViscousSource must act on the radial component.");
 
   if (_use_deviatoric_terms)
     _var.computeCellGradients();
 
   const auto get_velocity_var =
       [this](const std::string & param_name) -> MooseLinearVariableFVReal *
   {
     auto & var = _fe_problem.getVariable(_tid, getParam<SolverVariableName>(param_name));
     auto ptr = dynamic_cast<MooseLinearVariableFVReal *>(&var);
     if (!ptr)
       paramError(param_name, "The supplied variable must be a MooseLinearVariableFVReal.");
     return ptr;
   };
 
   if (_use_deviatoric_terms)
   {
     if (isParamValid("u"))
       _velocity_vars[0] = get_velocity_var("u");
     if (isParamValid("v"))
       _velocity_vars[1] = get_velocity_var("v");
 
     if (!_velocity_vars[0])
       paramError("u", "The x-velocity must be provided when using deviatoric terms.");
     if (!_velocity_vars[1])
       paramError("v", "The y-velocity must be provided when using deviatoric terms.");
 
     for (const auto dir : make_range(_dim))
       _velocity_vars[dir]->computeCellGradients();
   }
 }

Member Function Documentation

◆ computeMatrixContribution()

Real LinearFVRZViscousSource::computeMatrixContribution ( )

overrideprotectedvirtual

Implements LinearFVElementalKernel.

Definition at line 86 of file LinearFVRZViscousSource.C.

 {
   const Real r = _current_elem_info->centroid()(_rz_radial_coord);
   mooseAssert(r > 0, "Axisymmetric control volumes should not sit on the axis (r = 0).");
 
   const Real mu = _mu(makeElemArg(_current_elem_info->elem()), determineState());
   return mu * _stress_multiplier * _current_elem_volume / (r * r);
 }

◆ computeRightHandSideContribution()

Real LinearFVRZViscousSource::computeRightHandSideContribution ( )

overrideprotectedvirtual

Implements LinearFVElementalKernel.

Definition at line 96 of file LinearFVRZViscousSource.C.

 {
   if (!_use_deviatoric_terms)
     return 0.0;
 
   const auto state = determineState();
   Real divergence = 0.0;
   for (const auto dir : make_range(_dim))
     divergence += velocityVar(dir).gradSln(*_current_elem_info, state)(dir);
 
   const Real r = _current_elem_info->centroid()(_rz_radial_coord);
   mooseAssert(r > 0, "Axisymmetric control volumes should not sit on the axis (r = 0).");
 
   const Real radial_value = velocityVar(_rz_radial_coord).getElemValue(*_current_elem_info, state);
   divergence += radial_value / r;
 
   const auto elem_arg = makeElemArg(_current_elem_info->elem());
   const Real mu = _mu(elem_arg, state);
 
   return (2.0 / 3.0) * mu * divergence * _current_elem_volume / r;
 }

◆ validParams()

InputParameters LinearFVRZViscousSource::validParams ( )

static

Definition at line 20 of file LinearFVRZViscousSource.C.

 {
   InputParameters params = LinearFVElementalKernel::validParams();
   params.addClassDescription(
       "Adds the axisymmetric viscous source term mu * u_r / r^2 to the FV momentum equations.");
   params.addRequiredParam<MooseFunctorName>(NS::mu, "Dynamic viscosity functor.");
   MooseEnum component("x=0 y=1 z=2");
   params.addRequiredParam<MooseEnum>(
       "momentum_component", component, "Momentum component this kernel contributes to.");
   params.addParam<SolverVariableName>("u", "The velocity in the x direction.");
   params.addParam<SolverVariableName>("v", "The velocity in the y direction.");
   params.addParam<bool>(
       "use_deviatoric_terms",
       false,
       "Include the deviatoric correction (-2/3 div(u)) in the axisymmetric term.");
   return params;
 }

◆ velocityVar()

const MooseLinearVariableFVReal & LinearFVRZViscousSource::velocityVar ( unsigned int dir ) const

private

Helper to access the velocity variable for a given direction.

Definition at line 119 of file LinearFVRZViscousSource.C.

Referenced by computeRightHandSideContribution().

 {
   mooseAssert(dir < _velocity_vars.size() && _velocity_vars[dir],
               "Velocity variable for requested direction is not available.");
   return *_velocity_vars[dir];
 }

Member Data Documentation

◆ _component

const unsigned int LinearFVRZViscousSource::_component

private

Momentum component this source acts on (should equal the radial direction)

Definition at line 36 of file LinearFVRZViscousSource.h.

◆ _coord_type

const Moose::CoordinateSystemType LinearFVRZViscousSource::_coord_type

private

Coordinate system of the active blocks (must be COORD_RZ)

Definition at line 44 of file LinearFVRZViscousSource.h.

◆ _dim

const unsigned int LinearFVRZViscousSource::_dim

private

Spatial dimension of the mesh.

Definition at line 40 of file LinearFVRZViscousSource.h.

Referenced by computeRightHandSideContribution().

◆ _mu

const Moose::Functor<Real>& LinearFVRZViscousSource::_mu

private

Dynamic viscosity functor evaluated at each element.

Definition at line 34 of file LinearFVRZViscousSource.h.

Referenced by computeMatrixContribution(), and computeRightHandSideContribution().

◆ _rz_radial_coord

const unsigned int LinearFVRZViscousSource::_rz_radial_coord

private

Index of the radial coordinate for the current mesh (0 -> x, 1 -> y, ...)

Definition at line 38 of file LinearFVRZViscousSource.h.

Referenced by computeMatrixContribution(), and computeRightHandSideContribution().

◆ _stress_multiplier

const Real LinearFVRZViscousSource::_stress_multiplier

private

Precomputed factor (1 or 2) multiplying the implicit hoop term.

Definition at line 46 of file LinearFVRZViscousSource.h.

Referenced by computeMatrixContribution().

◆ _use_deviatoric_terms

const bool LinearFVRZViscousSource::_use_deviatoric_terms

private

Whether the deviatoric correction (-2/3 div u) is requested.

Definition at line 42 of file LinearFVRZViscousSource.h.

Referenced by computeRightHandSideContribution().

◆ _velocity_vars

std::array<MooseLinearVariableFVReal *, 2> LinearFVRZViscousSource::_velocity_vars

private

Cached pointers to the velocity components required to build divergence.

Definition at line 49 of file LinearFVRZViscousSource.h.

Referenced by velocityVar().

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

navier_stokes/include/linearfvkernels/LinearFVRZViscousSource.h
navier_stokes/src/linearfvkernels/LinearFVRZViscousSource.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ LinearFVRZViscousSource()

Member Function Documentation

◆ computeMatrixContribution()

◆ computeRightHandSideContribution()

◆ validParams()

◆ velocityVar()

Member Data Documentation

◆ _component

◆ _coord_type

◆ _dim

◆ _mu

◆ _rz_radial_coord

◆ _stress_multiplier

◆ _use_deviatoric_terms

◆ _velocity_vars