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

Kernel that adds contributions to the source and the sink of the turbulent kinetic energy discretized using the finite volume method to a linear system. More...

#include <LinearFVTKESourceSink.h>

Inheritance diagram for LinearFVTKESourceSink:
[legend]

Public Types

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

Public Member Functions

 LinearFVTKESourceSink (const InputParameters &params)
 Class constructor. More...
 
virtual void initialSetup () override
 
virtual Real computeMatrixContribution () override
 
virtual Real computeRightHandSideContribution () override
 
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< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
MooseAppgetMooseApp () const
 
const std::string & type () const
 
virtual const std::string & name () const
 
std::string typeAndName () const
 
std::string errorPrefix (const std::string &error_type) const
 
void callMooseError (std::string msg, const bool with_prefix) const
 
MooseObjectParameterName uniqueParameterName (const std::string &parameter_name) const
 
const InputParametersparameters () const
 
MooseObjectName uniqueName () const
 
const T & getParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const T * queryParam (const std::string &name) const
 
const T & getRenamedParam (const std::string &old_name, const std::string &new_name) const
 
getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool isParamValid (const std::string &name) const
 
bool isParamSetByUser (const std::string &nm) const
 
void paramError (const std::string &param, Args... args) const
 
void paramWarning (const std::string &param, Args... args) const
 
void paramInfo (const std::string &param, Args... args) const
 
void connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const
 
void mooseError (Args &&... args) const
 
void mooseErrorNonPrefixed (Args &&... args) const
 
void mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const
 
void mooseWarning (Args &&... args) const
 
void mooseWarningNonPrefixed (Args &&... args) const
 
void mooseDeprecated (Args &&... args) const
 
void mooseInfo (Args &&... args) const
 
std::string getDataFileName (const std::string &param) const
 
std::string getDataFileNameByName (const std::string &relative_path) const
 
std::string getDataFilePath (const std::string &relative_path) const
 
virtual void timestepSetup ()
 
virtual void jacobianSetup ()
 
virtual void residualSetup ()
 
virtual void subdomainSetup ()
 
virtual void customSetup (const ExecFlagType &)
 
const ExecFlagEnumgetExecuteOnEnum () const
 
const FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (const FunctionName &name) const
 
bool hasFunction (const std::string &param_name) const
 
bool hasFunctionByName (const FunctionName &name) const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const T & getUserObject (const std::string &param_name, bool is_dependency=true) const
 
const T & getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectgetUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const
 
bool isImplicit ()
 
Moose::StateArg determineState () const
 
bool isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessor (const std::string &param_name, const unsigned int index=0) const
 
bool hasPostprocessorByName (const PostprocessorName &name) const
 
std::size_t coupledPostprocessors (const std::string &param_name) const
 
const PostprocessorName & getPostprocessorName (const std::string &param_name, const unsigned int index=0) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
bool hasVectorPostprocessor (const std::string &param_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool hasVectorPostprocessorByName (const VectorPostprocessorName &name) const
 
const VectorPostprocessorName & getVectorPostprocessorName (const std::string &param_name) const
 
void setRandomResetFrequency (ExecFlagType exec_flag)
 
unsigned long getRandomLong () const
 
Real getRandomReal () const
 
unsigned int getSeed (std::size_t id)
 
unsigned int getMasterSeed () const
 
bool isNodal () const
 
ExecFlagType getResetOnTime () const
 
void setRandomDataPointer (RandomData *random_data)
 
virtual void meshChanged ()
 
void useVectorTag (const TagName &tag_name, VectorTagsKey)
 
void useVectorTag (TagID tag_id, VectorTagsKey)
 
void useMatrixTag (const TagName &tag_name, MatrixTagsKey)
 
void useMatrixTag (TagID tag_id, MatrixTagsKey)
 
bool isVectorTagged ()
 
bool isMatrixTagged ()
 
bool hasVectorTags () const
 
const std::set< TagID > & getVectorTags (VectorTagsKey) const
 
const std::set< TagID > & getMatrixTags (MatrixTagsKey) const
 
const std::vector< SubdomainName > & blocks () const
 
unsigned int numBlocks () const
 
virtual const std::set< SubdomainID > & blockIDs () const
 
unsigned int blocksMaxDimension () const
 
bool hasBlocks (const SubdomainName &name) const
 
bool hasBlocks (const std::vector< SubdomainName > &names) const
 
bool hasBlocks (const std::set< SubdomainName > &names) const
 
bool hasBlocks (SubdomainID id) const
 
bool hasBlocks (const std::vector< SubdomainID > &ids) const
 
bool hasBlocks (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::set< SubdomainID > &ids) const
 
bool isBlockSubset (const std::vector< SubdomainID > &ids) const
 
bool hasBlockMaterialProperty (const std::string &prop_name)
 
const std::set< SubdomainID > & meshBlockIDs () const
 
virtual bool blockRestricted () const
 
virtual void checkVariable (const MooseVariableFieldBase &variable) const
 
MooseVariableBasemooseVariableBase () 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 PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributed (const std::string &param_name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
bool isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

static InputParameters validParams ()
 
static std::string deduceFunctorName (const std::string &name, const InputParameters &params)
 
static void setRMParamsAdvection (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short conditional_extended_layers)
 
static void setRMParamsDiffusion (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short conditional_extended_layers)
 

Public Attributes

const ConsoleStream _console
 

Protected Member Functions

std::string deduceFunctorName (const std::string &name) const
 
void requestVariableCellGradient (const std::string &variable_name)
 
virtual void addUserObjectDependencyHelper (const UserObject &) const
 
virtual void addPostprocessorDependencyHelper (const PostprocessorName &) const
 
virtual void addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &) const
 
T & declareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< T > declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const T & getRestartableData (const std::string &data_name) const
 
T & declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
T & declareRecoverableData (const std::string &data_name, Args &&... args)
 
T & declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
T & declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTag (Assembly &assembly, unsigned int ivar, ResidualTagType tag_type)
 
void prepareVectorTagNeighbor (Assembly &assembly, unsigned int ivar)
 
void prepareVectorTagLower (Assembly &assembly, unsigned int ivar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTag (Assembly &assembly, unsigned int ivar, unsigned int jvar, DenseMatrix< Number > &k) const
 
void prepareMatrixTagNonlocal (Assembly &assembly, unsigned int ivar, unsigned int jvar)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type)
 
void prepareMatrixTagNeighbor (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, DenseMatrix< Number > &k) const
 
void prepareMatrixTagLower (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::ConstraintJacobianType type)
 
void accumulateTaggedLocalResidual ()
 
void assignTaggedLocalResidual ()
 
void accumulateTaggedLocalMatrix ()
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, const DenseMatrix< Number > &k)
 
void accumulateTaggedLocalMatrix (Assembly &assembly, unsigned int ivar, unsigned int jvar, Moose::DGJacobianType type, const DenseMatrix< Number > &k)
 
void accumulateTaggedNonlocalMatrix ()
 
void assignTaggedLocalMatrix ()
 
void addResiduals (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const DenseVector< T > &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResiduals (Assembly &assembly, const ADResidualsPacket &packet)
 
void addResidualsAndJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobian (Assembly &assembly, const ADResidualsPacket &packet)
 
void addJacobian (Assembly &assembly, DenseMatrix< Real > &local_k, const std::vector< dof_id_type > &row_indices, const std::vector< dof_id_type > &column_indices, Real scaling_factor)
 
void addResidualsWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addResidualsAndJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianWithoutConstraints (Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
 
void addJacobianElement (Assembly &assembly, Real value, dof_id_type row_index, dof_id_type column_index, Real scaling_factor)
 
void setResidual (SystemBase &sys, const T &residual, MooseVariableFE< T > &var)
 
void setResidual (SystemBase &sys, Real residual, dof_id_type dof_index)
 
void setResidual (SystemBase &sys, SetResidualFunctor set_residual_functor)
 
virtual bool hasBlockMaterialPropertyHelper (const std::string &prop_name)
 
void initializeBlockRestrictable (const MooseObject *moose_object)
 
Moose::CoordinateSystemType getBlockCoordSystem ()
 
const Moose::Functor< T > & getFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctor (const std::string &name, THREAD_ID tid)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem)
 
const Moose::Functor< T > & getFunctor (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
bool isFunctor (const std::string &name) const
 
bool isFunctor (const std::string &name, const SubProblem &subproblem) const
 
Moose::ElemArg makeElemArg (const Elem *elem, bool correct_skewnewss=false) const
 
void checkFunctorSupportsSideIntegration (const std::string &name, bool qp_integration)
 
virtual const OutputTools< Real >::VariableValuevalue ()
 
virtual const OutputTools< Real >::VariableValuevalueOld ()
 
virtual const OutputTools< Real >::VariableValuevalueOlder ()
 
virtual const OutputTools< Real >::VariableValuedot ()
 
virtual const OutputTools< Real >::VariableValuedotDot ()
 
virtual const OutputTools< Real >::VariableValuedotOld ()
 
virtual const OutputTools< Real >::VariableValuedotDotOld ()
 
virtual const VariableValuedotDu ()
 
virtual const VariableValuedotDotDu ()
 
virtual const OutputTools< Real >::VariableGradientgradient ()
 
virtual const OutputTools< Real >::VariableGradientgradientOld ()
 
virtual const OutputTools< Real >::VariableGradientgradientOlder ()
 
virtual const OutputTools< Real >::VariableSecondsecond ()
 
virtual const OutputTools< Real >::VariableSecondsecondOld ()
 
virtual const OutputTools< Real >::VariableSecondsecondOlder ()
 
virtual const OutputTools< Real >::VariableTestSecondsecondTest ()
 
virtual const OutputTools< Real >::VariableTestSecondsecondTestFace ()
 
virtual const OutputTools< Real >::VariablePhiSecondsecondPhi ()
 
virtual const OutputTools< Real >::VariablePhiSecondsecondPhiFace ()
 

Protected Attributes

const unsigned int _dim
 The dimension of the domain. More...
 
const Moose::Functor< Real > & _u_var
 x-velocity More...
 
const Moose::Functor< Real > * _v_var
 y-velocity More...
 
const Moose::Functor< Real > * _w_var
 z-velocity More...
 
const Moose::Functor< Real > & _epsilon
 epsilon - dissipation rate of TKE More...
 
const Moose::Functor< Real > & _rho
 Density. More...
 
const Moose::Functor< Real > & _mu
 Dynamic viscosity. More...
 
const Moose::Functor< Real > & _mu_t
 Turbulent dynamic viscosity. More...
 
const std::vector< BoundaryName > & _wall_boundary_names
 Wall boundaries. More...
 
const bool _linearized_model
 Linearized model? More...
 
const NS::WallTreatmentEnum _wall_treatment
 Method used for wall treatment. More...
 
const Real _C_mu
 C_mu constant. More...
 
const Real _C_pl
 Production Limiter Constant. More...
 
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
 
const std::string _type
 
const std::string _name
 
const InputParameters_pars
 
Factory_factory
 
ActionFactory_action_factory
 
const ExecFlagEnum_execute_enum
 
const ExecFlagType_current_execute_flag
 
const InputParameters_ti_params
 
FEProblemBase_ti_feproblem
 
bool _is_implicit
 
Real_t
 
const Real_t_old
 
int_t_step
 
Real_dt
 
Real_dt_old
 
bool _is_transient
 
MooseApp_restartable_app
 
const std::string _restartable_system_name
 
const THREAD_ID _restartable_tid
 
const bool _restartable_read_only
 
FEProblemBase_mci_feproblem
 
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
 
std::unordered_set< const Elem * > _wall_bounded
 
std::map< const Elem *, std::vector< Real > > _dist
 
std::map< const Elem *, std::vector< const FaceInfo * > > _face_infos
 

Detailed Description

Kernel that adds contributions to the source and the sink of the turbulent kinetic energy discretized using the finite volume method to a linear system.

Definition at line 19 of file LinearFVTKESourceSink.h.

Constructor & Destructor Documentation

◆ LinearFVTKESourceSink()

LinearFVTKESourceSink::LinearFVTKESourceSink ( const InputParameters params)

Class constructor.

Parameters
paramsThe InputParameters for the kernel.

Definition at line 50 of file LinearFVTKESourceSink.C.

51  : LinearFVElementalKernel(params),
53  _u_var(getFunctor<Real>("u")),
54  _v_var(params.isParamValid("v") ? &(getFunctor<Real>("v")) : nullptr),
55  _w_var(params.isParamValid("w") ? &(getFunctor<Real>("w")) : nullptr),
56  _epsilon(getFunctor<Real>(NS::TKED)),
57  _rho(getFunctor<Real>(NS::density)),
58  _mu(getFunctor<Real>(NS::mu)),
59  _mu_t(getFunctor<Real>(NS::mu_t)),
60  _wall_boundary_names(getParam<std::vector<BoundaryName>>("walls")),
61  _linearized_model(getParam<bool>("linearized_model")),
62  _wall_treatment(getParam<MooseEnum>("wall_treatment").getEnum<NS::WallTreatmentEnum>()),
63  _C_mu(getParam<Real>("C_mu")),
64  _C_pl(getParam<Real>("C_pl"))
65 {
66  if (_dim >= 2 && !_v_var)
67  paramError("v", "In two or more dimensions, the v velocity must be supplied!");
68 
69  if (_dim >= 3 && !_w_var)
70  paramError("w", "In three or more dimensions, the w velocity must be supplied!");
71 
72  // // Strain tensor term requires velocity gradients;
73  if (dynamic_cast<const MooseLinearVariableFV<Real> *>(&_u_var))
74  requestVariableCellGradient(getParam<MooseFunctorName>("u"));
75  if (dynamic_cast<const MooseLinearVariableFV<Real> *>(_v_var))
76  requestVariableCellGradient(getParam<MooseFunctorName>("v"));
77  if (dynamic_cast<const MooseLinearVariableFV<Real> *>(_w_var))
78  requestVariableCellGradient(getParam<MooseFunctorName>("w"));
79 }
virtual MooseMesh & mesh()=0
const Moose::Functor< Real > * _w_var
z-velocity
static const std::string mu_t
Definition: NS.h:125
SubProblem & _subproblem
void requestVariableCellGradient(const std::string &variable_name)
const Real _C_pl
Production Limiter Constant.
const Moose::Functor< Real > & _u_var
x-velocity
static const std::string density
Definition: NS.h:33
WallTreatmentEnum
Wall treatment options.
Definition: NS.h:182
const NS::WallTreatmentEnum _wall_treatment
Method used for wall treatment.
const Moose::Functor< Real > & _rho
Density.
const std::vector< BoundaryName > & _wall_boundary_names
Wall boundaries.
const Real _C_mu
C_mu constant.
const Moose::Functor< Real > & _epsilon
epsilon - dissipation rate of TKE
const unsigned int _dim
The dimension of the domain.
virtual unsigned int dimension() const
static const std::string mu
Definition: NS.h:123
const Moose::Functor< Real > & _mu
Dynamic viscosity.
const T & getParam(const std::string &name) const
void paramError(const std::string &param, Args... args) const
LinearFVElementalKernel(const InputParameters &params)
const bool _linearized_model
Linearized model?
static const std::string TKED
Definition: NS.h:177
const Moose::Functor< Real > & _mu_t
Turbulent dynamic viscosity.
bool isParamValid(const std::string &name) const
const Moose::Functor< Real > * _v_var
y-velocity

Member Function Documentation

◆ computeMatrixContribution()

Real LinearFVTKESourceSink::computeMatrixContribution ( )
overridevirtual

Implements LinearFVElementalKernel.

Definition at line 92 of file LinearFVTKESourceSink.C.

93 {
94  /*
95  Matrix contribution:
96  - Computes near-wall TKE destruction
97  - Computes bulk TKE destruction
98  */
99 
100  // Useful variables
101  const auto state = determineState();
102  const auto elem_arg = makeElemArg(_current_elem_info->elem());
103  const Real rho = _rho(elem_arg, state);
104 
105  if (_wall_bounded.find(_current_elem_info->elem()) != _wall_bounded.end())
106  {
107 
108  // Useful variables
109  const Real mu = _mu(elem_arg, state);
110  const Real TKE = _var.getElemValue(*_current_elem_info, state);
111 
112  // Assembly variables
113  Real production = 0.0;
114  Real destruction = 0.0;
115  Real tot_weight = 0.0;
116  std::vector<Real> y_plus_vec, velocity_grad_norm_vec;
117 
118  // Get near-wall element velocity vector
119  RealVectorValue velocity(_u_var(elem_arg, state));
120  if (_v_var)
121  velocity(1) = (*_v_var)(elem_arg, state);
122  if (_w_var)
123  velocity(2) = (*_w_var)(elem_arg, state);
124 
125  // Get faceInfo and distance vector for between all walls and elements
126  const auto & face_info_vec = libmesh_map_find(_face_infos, _current_elem_info->elem());
127  const auto & distance_vec = libmesh_map_find(_dist, _current_elem_info->elem());
128  mooseAssert(distance_vec.size(), "Should have found a distance vector");
129  mooseAssert(distance_vec.size() == face_info_vec.size(),
130  "Should be as many distance vectors as face info vectors");
131 
132  // Populate y+, near wall velocity gradient, and update weights
133  for (unsigned int i = 0; i < distance_vec.size(); i++)
134  {
135  const auto parallel_speed = NS::computeSpeed<Real>(
136  velocity - velocity * face_info_vec[i]->normal() * face_info_vec[i]->normal());
137  const auto distance = distance_vec[i];
138 
139  Real y_plus;
140  if (_wall_treatment == NS::WallTreatmentEnum::NEQ) // Non-equilibrium --> Non-iterative
141  y_plus = distance * std::sqrt(std::sqrt(_C_mu) * TKE) * rho / mu;
142  else // Equilibrium --> Iterative
143  y_plus = NS::findyPlus<Real>(mu, rho, std::max(parallel_speed, 1e-10), distance);
144 
145  y_plus_vec.push_back(y_plus);
146  const Real velocity_grad_norm = parallel_speed / distance;
147  velocity_grad_norm_vec.push_back(velocity_grad_norm);
148  tot_weight += 1.0;
149  }
150 
151  // Compute near-wall production and destruction
152  for (unsigned int i = 0; i < y_plus_vec.size(); i++)
153  {
154  const auto y_plus = y_plus_vec[i];
155 
156  const auto fi = face_info_vec[i];
157  const bool defined_on_elem_side = _var.hasFaceSide(*fi, true);
158  const Elem * const loc_elem = defined_on_elem_side ? &fi->elem() : fi->neighborPtr();
159  const Moose::FaceArg facearg = {
160  fi, Moose::FV::LimiterType::CentralDifference, false, false, loc_elem, nullptr};
161 
162  const Real wall_mut = _mu_t(facearg, state);
163  const Real wall_mu = _mu(facearg, state);
164  const auto tau_w = (wall_mut + wall_mu) * velocity_grad_norm_vec[i];
165  const auto destruction_visc = 2.0 * wall_mu / Utility::pow<2>(distance_vec[i]) / tot_weight;
166  const auto destruction_log = std::pow(_C_mu, 0.75) * rho * std::pow(TKE, 0.5) /
167  (NS::von_karman_constant * distance_vec[i]) / tot_weight;
168 
169  if (y_plus < 11.25)
170  destruction += destruction_visc;
171  else
172  {
173  destruction += destruction_log;
174  production += tau_w * std::pow(_C_mu, 0.25) / std::sqrt(TKE) /
175  (NS::von_karman_constant * distance_vec[i]) / tot_weight;
176  }
177  }
178 
179  // Assign terms to matrix to solve implicitly (they get multiplied by TKE)
180  return (destruction - production) * _current_elem_volume;
181  }
182  else
183  {
184  // Compute destruction
185  const auto destruction =
186  rho * _epsilon(elem_arg, state) / _var.getElemValue(*_current_elem_info, state);
187 
188  // Assign to matrix to solve implicitly
189  return destruction * _current_elem_volume;
190  }
191 }
const ElemInfo * _current_elem_info
const Moose::Functor< Real > * _w_var
z-velocity
std::map< const Elem *, std::vector< Real > > _dist
static constexpr Real von_karman_constant
Definition: NS.h:191
Moose::StateArg determineState() const
const Elem * elem() const
const Moose::Functor< Real > & _u_var
x-velocity
MooseLinearVariableFV< Real > & _var
static const std::string TKE
Definition: NS.h:176
const NS::WallTreatmentEnum _wall_treatment
Method used for wall treatment.
const Moose::Functor< Real > & _rho
Density.
Real distance(const Point &p)
Moose::ElemArg makeElemArg(const Elem *elem, bool correct_skewnewss=false) const
template Real findyPlus< Real >(const Real &mu, const Real &rho, const Real &u, Real dist)
const Real _C_mu
C_mu constant.
const Moose::Functor< Real > & _epsilon
epsilon - dissipation rate of TKE
std::unordered_set< const Elem * > _wall_bounded
static const std::string mu
Definition: NS.h:123
const Moose::Functor< Real > & _mu
Dynamic viscosity.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Real getElemValue(const ElemInfo &elem_info, const StateArg &state) const
static const std::string velocity
Definition: NS.h:45
template Real computeSpeed< Real >(const libMesh::VectorValue< Real > &velocity)
std::map< const Elem *, std::vector< const FaceInfo * > > _face_infos
MooseUnits pow(const MooseUnits &, int)
const Moose::Functor< Real > & _mu_t
Turbulent dynamic viscosity.
virtual bool hasFaceSide(const FaceInfo &fi, const bool fi_elem_side) const override
const Moose::Functor< Real > * _v_var
y-velocity

◆ computeRightHandSideContribution()

Real LinearFVTKESourceSink::computeRightHandSideContribution ( )
overridevirtual

Implements LinearFVElementalKernel.

Definition at line 194 of file LinearFVTKESourceSink.C.

195 {
196 
197  // Useful variables
198  const auto state = determineState();
199  const auto elem_arg = makeElemArg(_current_elem_info->elem());
200 
201  if (_wall_bounded.find(_current_elem_info->elem()) != _wall_bounded.end()) // Wall bounded
202  return 0.0; // Do nothing
203  else // Not wall bounded
204  {
205  // Compute TKE production
206  const auto symmetric_strain_tensor_sq_norm =
208 
209  auto production = _mu_t(elem_arg, state) * symmetric_strain_tensor_sq_norm;
210 
211  // k-Production limiter (needed for flows with stagnation zones)
212  const Real production_limit = _C_pl * _rho(elem_arg, state) * _epsilon(elem_arg, state);
213 
214  // Apply production limiter
215  production = std::min(production, production_limit);
216 
217  // Assign production to RHS
218  return production * _current_elem_volume;
219  }
220 }
const ElemInfo * _current_elem_info
const Moose::Functor< Real > * _w_var
z-velocity
const Real _C_pl
Production Limiter Constant.
Moose::StateArg determineState() const
const Elem * elem() const
const Moose::Functor< Real > & _u_var
x-velocity
const Moose::Functor< Real > & _rho
Density.
Moose::ElemArg makeElemArg(const Elem *elem, bool correct_skewnewss=false) const
const Moose::Functor< Real > & _epsilon
epsilon - dissipation rate of TKE
std::unordered_set< const Elem * > _wall_bounded
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
const Moose::Functor< Real > & _mu_t
Turbulent dynamic viscosity.
template Real computeShearStrainRateNormSquared< Real >(const Moose::Functor< Real > &u, const Moose::Functor< Real > *v, const Moose::Functor< Real > *w, const Moose::ElemArg &elem_arg, const Moose::StateArg &state)
const Moose::Functor< Real > * _v_var
y-velocity

◆ initialSetup()

void LinearFVTKESourceSink::initialSetup ( )
overridevirtual

Reimplemented from LinearFVElementalKernel.

Definition at line 82 of file LinearFVTKESourceSink.C.

83 {
89 }
std::map< const Elem *, std::vector< Real > > _dist
SubProblem & _subproblem
void getWallBoundedElements(const std::vector< BoundaryName > &wall_boundary_name, const FEProblemBase &fe_problem, const SubProblem &subproblem, const std::set< SubdomainID > &block_ids, std::unordered_set< const Elem *> &wall_bounded)
Map marking wall bounded elements The map passed in wall_bounded_map gets cleared and re-populated...
virtual const std::set< SubdomainID > & blockIDs() const
const std::vector< BoundaryName > & _wall_boundary_names
Wall boundaries.
std::unordered_set< const Elem * > _wall_bounded
void getWallDistance(const std::vector< BoundaryName > &wall_boundary_name, const FEProblemBase &fe_problem, const SubProblem &subproblem, const std::set< SubdomainID > &block_ids, std::map< const Elem *, std::vector< Real >> &dist_map)
Map storing wall ditance for near-wall marked elements The map passed in dist_map gets cleared and re...
void getElementFaceArgs(const std::vector< BoundaryName > &wall_boundary_name, const FEProblemBase &fe_problem, const SubProblem &subproblem, const std::set< SubdomainID > &block_ids, std::map< const Elem *, std::vector< const FaceInfo *>> &face_info_map)
Map storing face arguments to wall bounded faces The map passed in face_info_map gets cleared and re-...
FEProblemBase & _fe_problem
virtual void initialSetup()
std::map< const Elem *, std::vector< const FaceInfo * > > _face_infos

◆ validParams()

InputParameters LinearFVTKESourceSink::validParams ( )
static

Definition at line 18 of file LinearFVTKESourceSink.C.

19 {
21  params.addClassDescription("Elemental kernel to compute the production and destruction "
22  " terms of turbulent kinetic energy (TKE).");
23  params.addRequiredParam<MooseFunctorName>("u", "The velocity in the x direction.");
24  params.addParam<MooseFunctorName>("v", "The velocity in the y direction.");
25  params.addParam<MooseFunctorName>("w", "The velocity in the z direction.");
26  params.addRequiredParam<MooseFunctorName>(NS::TKED,
27  "Coupled turbulent kinetic energy dissipation rate.");
28  params.addRequiredParam<MooseFunctorName>(NS::density, "Fluid density");
29  params.addRequiredParam<MooseFunctorName>(NS::mu, "Dynamic viscosity.");
30  params.addRequiredParam<MooseFunctorName>(NS::mu_t, "Turbulent viscosity.");
31 
32  params.addParam<std::vector<BoundaryName>>(
33  "walls", {}, "Boundaries that correspond to solid walls.");
34  params.addParam<bool>(
35  "linearized_model",
36  true,
37  "Boolean to determine if the problem should be use in a linear or nonlinear solve.");
38  MooseEnum wall_treatment("eq_newton eq_incremental eq_linearized neq", "neq");
39 
40  params.addParam<MooseEnum>("wall_treatment",
41  wall_treatment,
42  "The method used for computing the wall functions "
43  "'eq_newton', 'eq_incremental', 'eq_linearized', 'neq'");
44  params.addParam<Real>("C_mu", 0.09, "Coupled turbulent kinetic energy closure.");
45  params.addParam<Real>("C_pl", 10.0, "Production Limiter Constant Multiplier.");
46 
47  return params;
48 }
static const std::string mu_t
Definition: NS.h:125
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
static InputParameters validParams()
static const std::string density
Definition: NS.h:33
void addRequiredParam(const std::string &name, const std::string &doc_string)
static const std::string mu
Definition: NS.h:123
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
static const std::string TKED
Definition: NS.h:177
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _C_mu

const Real LinearFVTKESourceSink::_C_mu
protected

C_mu constant.

Definition at line 69 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution().

◆ _C_pl

const Real LinearFVTKESourceSink::_C_pl
protected

Production Limiter Constant.

Definition at line 72 of file LinearFVTKESourceSink.h.

Referenced by computeRightHandSideContribution().

◆ _dim

const unsigned int LinearFVTKESourceSink::_dim
protected

The dimension of the domain.

Definition at line 38 of file LinearFVTKESourceSink.h.

Referenced by LinearFVTKESourceSink().

◆ _dist

std::map<const Elem *, std::vector<Real> > LinearFVTKESourceSink::_dist
protected

Definition at line 77 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution(), and initialSetup().

◆ _epsilon

const Moose::Functor<Real>& LinearFVTKESourceSink::_epsilon
protected

epsilon - dissipation rate of TKE

Definition at line 48 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution(), and computeRightHandSideContribution().

◆ _face_infos

std::map<const Elem *, std::vector<const FaceInfo *> > LinearFVTKESourceSink::_face_infos
protected

Definition at line 78 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution(), and initialSetup().

◆ _linearized_model

const bool LinearFVTKESourceSink::_linearized_model
protected

Linearized model?

Definition at line 63 of file LinearFVTKESourceSink.h.

◆ _mu

const Moose::Functor<Real>& LinearFVTKESourceSink::_mu
protected

Dynamic viscosity.

Definition at line 54 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution().

◆ _mu_t

const Moose::Functor<Real>& LinearFVTKESourceSink::_mu_t
protected

Turbulent dynamic viscosity.

Definition at line 57 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution(), and computeRightHandSideContribution().

◆ _rho

const Moose::Functor<Real>& LinearFVTKESourceSink::_rho
protected

Density.

Definition at line 51 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution(), and computeRightHandSideContribution().

◆ _u_var

const Moose::Functor<Real>& LinearFVTKESourceSink::_u_var
protected

◆ _v_var

const Moose::Functor<Real>* LinearFVTKESourceSink::_v_var
protected

◆ _w_var

const Moose::Functor<Real>* LinearFVTKESourceSink::_w_var
protected

◆ _wall_boundary_names

const std::vector<BoundaryName>& LinearFVTKESourceSink::_wall_boundary_names
protected

Wall boundaries.

Definition at line 60 of file LinearFVTKESourceSink.h.

Referenced by initialSetup().

◆ _wall_bounded

std::unordered_set<const Elem *> LinearFVTKESourceSink::_wall_bounded
protected

Maps for wall treatement

Definition at line 76 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution(), computeRightHandSideContribution(), and initialSetup().

◆ _wall_treatment

const NS::WallTreatmentEnum LinearFVTKESourceSink::_wall_treatment
protected

Method used for wall treatment.

Definition at line 66 of file LinearFVTKESourceSink.h.

Referenced by computeMatrixContribution().


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