https://mooseframework.inl.gov
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Static Public Member Functions | Public Attributes | Static Public Attributes | Protected Member Functions | Protected Attributes | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
LinearFVTurbulentAdvection Class Reference

An advection kernel that implements the advection term for the turbulent variables limited for the first cells near the wall. More...

#include <LinearFVTurbulentAdvection.h>

Inheritance diagram for LinearFVTurbulentAdvection:
[legend]

Public Types

enum  ResidualTagType
 
typedef DataFileName DataFileParameterType
 

Public Member Functions

 LinearFVTurbulentAdvection (const InputParameters &params)
 
virtual void addMatrixContribution () override
 
virtual void addRightHandSideContribution () override
 
virtual void initialSetup () override
 
virtual Real computeElemMatrixContribution () override
 
virtual Real computeNeighborMatrixContribution () override
 
virtual Real computeElemRightHandSideContribution () override
 
virtual Real computeNeighborRightHandSideContribution () override
 
virtual Real computeBoundaryMatrixContribution (const LinearFVBoundaryCondition &bc) override
 
virtual Real computeBoundaryRHSContribution (const LinearFVBoundaryCondition &bc) override
 
virtual void setupFaceData (const FaceInfo *face_info) override
 
virtual bool hasFaceSide (const FaceInfo &fi, bool fi_elem_side) const override
 
void setCurrentFaceArea (const Real area)
 
virtual Real computeBoundaryFlux (const LinearFVBoundaryCondition &bc)
 
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 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 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
 
const Moose::Functor< Real > * defaultFunctor (const std::string &name)
 
const Moose::Functor< ADReal > * defaultFunctor (const std::string &name)
 
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={})
 
virtual std::set< MooseVariableFieldBase * > checkVariables (const libMesh::Node &node, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
std::set< MooseVariableFieldBase * > checkVariables (const libMesh::Elem &element, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
void addMooseVariableDependency (MooseVariableFieldBase *var)
 
void addMooseVariableDependency (const std::vector< MooseVariableFieldBase * > &vars)
 
Moose::FaceArg makeFace (const FaceInfo &fi, const Moose::FV::LimiterType limiter_type, const bool elem_is_upwind, const bool correct_skewness=false, const Moose::StateArg *state_limiter=nullptr) const
 
Moose::FaceArg makeCDFace (const FaceInfo &fi, const bool correct_skewness=false) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObject (const std::string &param_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
bool hasUserObjectByName (const UserObjectName &object_name) const
 
const 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
 
const FVFaceInterpolationMethod & getFVFaceInterpolationMethod (const InterpolationMethodName &name) const
 
const FVAdvectedInterpolationMethod & getFVAdvectedInterpolationMethod (const InterpolationMethodName &name) const
 
bool hasFVInterpolationMethod (const InterpolationMethodName &name) 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

std::string deduceFunctorName (const std::string &name) const
 
Moose::FaceArg singleSidedFaceArg (const FaceInfo *fi, Moose::FV::LimiterType limiter_type=Moose::FV::LimiterType::CentralDifference, bool correct_skewness=false) 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, const 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 RhieChowMassFlux & _mass_flux_provider
 The Rhie-Chow user object that provides us with the face velocity.
 
const FaceInfo * _current_face_info
 
Real _current_face_area
 
FaceInfo::VarFaceNeighbors _current_face_type
 
bool _cached_matrix_contribution
 
bool _cached_rhs_contribution
 
const bool _force_boundary_execution
 
DenseVector< dof_id_type > _dof_indices
 
DenseMatrix< Real > _matrix_contribution
 
DenseVector< Real > _rhs_contribution
 
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
 
const InputParameters & _fvim_params
 
FEProblemBase & _fvim_feproblem
 
const THREAD_ID _fvim_tid
 

Private Member Functions

const Moose::Functor< T > * defaultFunctor (const std::string &name)
 
const Moose::Functor< Real > * defaultFunctor (const std::string &name)
 
const Moose::Functor< ADReal > * defaultFunctor (const std::string &name)
 
const Moose::Functor< T > & getFunctorByName (const std::string &name, SubProblem &subproblem, THREAD_ID tid)
 
const Moose::FunctionBase * getKokkosFunctionByNameHelper (const FunctionName &name) const
 
const UserObjectBase & getUserObjectFromFEProblem (const UserObjectName &object_name, const THREAD_ID tid=0) const
 
const T & castUserObject (const UserObjectBase &uo_base, const std::string &param_name="") const
 
void mooseObjectError (const std::string &param_name, std::stringstream &oss) const
 
const std::string & userObjectType (const UserObjectBase &uo) const
 
const std::string & userObjectName (const UserObjectBase &uo) const
 
const PostprocessorName & getPostprocessorNameInternal (const std::string &param_name, const unsigned int index, const bool allow_default_value=true) const
 
bool isDefaultPostprocessorValueByName (const PostprocessorName &name) const
 
PostprocessorValue getDefaultPostprocessorValueByName (const PostprocessorName &name) const
 
void checkParam (const std::string &param_name, const unsigned int index=std::numeric_limits< unsigned int >::max()) const
 
bool postprocessorsAdded () const
 
const VectorPostprocessorValue & getVectorPostprocessorByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name, bool broadcast, std::size_t t_index) const
 
const VectorPostprocessorContext< VectorPostprocessorValue > & getVectorPostprocessorContextByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name) const
 
bool vectorPostprocessorsAdded () const
 
RestartableDataValue & registerRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const
 
void registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)
 
RestartableData< T > & declareRestartableDataHelper (const std::string &data_name, void *context, Args &&... args) const
 
void prepareVectorTagInternal (Assembly &assembly, unsigned int ivar, const std::set< TagID > &vector_tags, const std::set< TagID > &absolute_value_vector_tags)
 
void checkForNans () const
 
virtual bool isADObject () const override
 
std::set< MooseVariableFieldBase * > checkVariablesHelper (const DofObjectType &dof_object, const std::set< MooseVariableFieldBase * > &vars_to_check)
 
const PostprocessorValue & getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const
 
const PostprocessorValue & getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const
 
const PostprocessorValue & getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const
 
const PostprocessorValue & getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const
 
void possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
void possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 

Static Private Member Functions

static const hit::Node * getHitNode (const InputParameters &params)
 
static std::string messagePrefix (const InputParameters &params, const bool hit_prefix)
 
static void setRMParams (const InputParameters &obj_params, InputParameters &rm_params, const unsigned short ghost_layers)
 
static void parameterError (const InputParameters &obj_params, const std::string &parameter_name, const std::string &function_name, const std::string &description)
 

Private Attributes

const FVAdvectedInterpolationMethod & _adv_interp_method
 The interpolation method to use for the advected quantity.
 
FVAdvectedInterpolationMethod::AdvectedSystemContribution _adv_interp_result
 Cached weights/correction for the current face (refreshed in setupFaceData)
 
VectorValue< Real > _elem_grad_storage
 Reusable gradient storage used when advected interpolation requires gradients.
 
VectorValue< Real > _neighbor_grad_storage
 
Real _mass_face_flux
 Container for the mass flux on the face which will be reused in the advection term's matrix and right hand side contribution.
 
const std::vector< BoundaryName > & _wall_boundary_names
 Wall boundaries.
 
std::unordered_set< const Elem * > _wall_bounded
 List for wall bounded elements.
 
const ParallelParamObject & _parent
 
const MooseBase & _si_moose_base
 
const FEProblemBase * _si_problem
 
ExecFlagEnum _empty_execute_enum
 
const MooseObject & _fni_object
 
const InputParameters & _fni_params
 
FEProblemBase & _fni_feproblem
 
const THREAD_ID _fni_tid
 
const MooseObject & _uoi_moose_object
 
const FEProblemBase & _uoi_feproblem
 
const THREAD_ID _uoi_tid
 
const std::string _ti_name
 
const MooseObject & _ppi_moose_object
 
const InputParameters & _ppi_params
 
const FEProblemBase & _ppi_feproblem
 
std::map< PostprocessorName, std::unique_ptr< PostprocessorValue > > _default_values
 
const bool _broadcast_by_default
 
const MooseObject & _vpi_moose_object
 
const FEProblemBase & _vpi_feproblem
 
const THREAD_ID _vpi_tid
 
RandomData * _random_data
 
MooseRandom * _generator
 
FEProblemBase & _ri_problem
 
const std::string _ri_name
 
unsigned int _master_seed
 
bool _is_nodal
 
ExecFlagType _reset_on
 
const Node *const & _curr_node
 
const Elem *const & _curr_element
 
const RestartableDataMapName _metaname
 
std::string _restartable_name
 
std::set< TagID > _vector_tags
 
std::set< TagID > _abs_vector_tags
 
std::set< TagID > _matrix_tags
 
std::set< TagID > _non_ref_vector_tags
 
std::set< TagID > _non_ref_abs_vector_tags
 
std::set< TagID > _ref_vector_tags
 
std::set< TagID > _ref_abs_vector_tags
 
const MooseObject & _moose_object
 
const InputParameters & _tag_params
 
std::vector< DenseVector< Number > * > _re_blocks
 
std::vector< DenseVector< Number > * > _absre_blocks
 
std::vector< DenseMatrix< Number > * > _ke_blocks
 
std::vector< Real > _absolute_residuals
 
std::set< SubdomainID > _blk_ids
 
std::vector< SubdomainID > _vec_ids
 
std::vector< SubdomainName > _blocks
 
const bool _blk_dual_restrictable
 
FEProblemBase * _blk_feproblem
 
const MooseMesh * _blk_mesh
 
const std::set< BoundaryID > _empty_boundary_ids
 
const std::set< BoundaryID > & _boundary_ids
 
THREAD_ID _blk_tid
 
const std::string & _blk_name
 
unsigned int _blk_dim
 
Moose::Kokkos::Array< ContiguousElementID > _kokkos_element_ids
 
Moose::Kokkos::Array< ContiguousNodeID > _kokkos_node_ids
 
Moose::Kokkos::Array< Moose::Kokkos::Pair< ContiguousElementID, unsigned int > > _kokkos_element_side_ids
 
const InputParameters & _fi_params
 
const std::string _fi_name
 
SubProblem *const _fi_subproblem
 
const THREAD_ID _fi_tid
 
std::vector< std::unique_ptr< Moose::Functor< Real > > > _default_real_functors
 
std::vector< std::unique_ptr< Moose::Functor< ADReal > > > _default_ad_real_functors
 
std::set< MooseVariableFieldBase * > _moose_variable_dependencies
 

Detailed Description

An advection kernel that implements the advection term for the turbulent variables limited for the first cells near the wall.

Definition at line 22 of file LinearFVTurbulentAdvection.h.

Constructor & Destructor Documentation

◆ LinearFVTurbulentAdvection()

LinearFVTurbulentAdvection::LinearFVTurbulentAdvection ( const InputParameters &  params)

Definition at line 38 of file LinearFVTurbulentAdvection.C.

39 : LinearFVFluxKernel(params),
40 FVInterpolationMethodInterface(this),
41 _mass_flux_provider(getUserObject<RhieChowMassFlux>("rhie_chow_user_object")),
42 _adv_interp_method(getFVAdvectedInterpolationMethod(
43 getParam<InterpolationMethodName>("advected_interp_method_name"))),
44 _mass_face_flux(0.0),
45 _wall_boundary_names(getParam<std::vector<BoundaryName>>("walls"))
46{
47 if (_adv_interp_method.needsGradients())
48 _var.computeCellGradients(_adv_interp_method.gradientLimiter());
49}
FVAdvectedInterpolationMethod::gradientLimiter
virtual Moose::FV::GradientLimiterType gradientLimiter() const
FVAdvectedInterpolationMethod::needsGradients
virtual bool needsGradients() const
FVInterpolationMethodInterface::getFVAdvectedInterpolationMethod
const FVAdvectedInterpolationMethod & getFVAdvectedInterpolationMethod(const InterpolationMethodName &name) const
LinearFVFluxKernel::_var
MooseLinearVariableFV< Real > & _var
LinearFVTurbulentAdvection::_mass_flux_provider
const RhieChowMassFlux & _mass_flux_provider
The Rhie-Chow user object that provides us with the face velocity.
Definition LinearFVTurbulentAdvection.h:50
LinearFVTurbulentAdvection::_wall_boundary_names
const std::vector< BoundaryName > & _wall_boundary_names
Wall boundaries.
Definition LinearFVTurbulentAdvection.h:68
LinearFVTurbulentAdvection::_adv_interp_method
const FVAdvectedInterpolationMethod & _adv_interp_method
The interpolation method to use for the advected quantity.
Definition LinearFVTurbulentAdvection.h:54
LinearFVTurbulentAdvection::_mass_face_flux
Real _mass_face_flux
Container for the mass flux on the face which will be reused in the advection term's matrix and right...
Definition LinearFVTurbulentAdvection.h:65
LinearFVFluxKernel::getParam
const T & getParam(const std::string &name) const

Member Function Documentation

◆ addMatrixContribution()

void LinearFVTurbulentAdvection::addMatrixContribution ( )
overridevirtual

Reimplemented from LinearFVFluxKernel.

Definition at line 60 of file LinearFVTurbulentAdvection.C.

61{
62 // Coumputing bounding map
63 const Elem * elem = _current_face_info->elemPtr();
64 const auto bounded_elem = _wall_bounded.find(elem) != _wall_bounded.end();
65 const Elem * neighbor = _current_face_info->neighborPtr();
66 const auto bounded_neigh = _wall_bounded.find(neighbor) != _wall_bounded.end();
67
68 // If we are on an internal face, we populate the four entries in the system matrix
69 // which touch the face
70 if (_current_face_type == FaceInfo::VarFaceNeighbors::BOTH)
71 {
72 // The dof ids of the variable corresponding to the element and neighbor
73 _dof_indices(0) = _current_face_info->elemInfo()->dofIndices()[_sys_num][_var_num];
74 _dof_indices(1) = _current_face_info->neighborInfo()->dofIndices()[_sys_num][_var_num];
75
76 // Compute the entries which will go to the neighbor (offdiagonal) and element
77 // (diagonal).
78 const auto elem_matrix_contribution = computeElemMatrixContribution();
79 const auto neighbor_matrix_contribution = computeNeighborMatrixContribution();
80
81 // Populate matrix
82 if (hasBlocks(_current_face_info->elemInfo()->subdomain_id()) && !(bounded_elem))
83 {
84 _matrix_contribution(0, 0) = elem_matrix_contribution;
85 _matrix_contribution(0, 1) = neighbor_matrix_contribution;
86 }
87
88 if (hasBlocks(_current_face_info->neighborInfo()->subdomain_id()) && !(bounded_neigh))
89 {
90 _matrix_contribution(1, 0) = -elem_matrix_contribution;
91 _matrix_contribution(1, 1) = -neighbor_matrix_contribution;
92 }
93 (*_linear_system.matrix).add_matrix(_matrix_contribution, _dof_indices.get_values());
94 }
95 // We are at a block boundary where the variable is not defined on one of the adjacent cells.
96 // We check if we have a boundary condition here
97 else if (_current_face_type == FaceInfo::VarFaceNeighbors::ELEM ||
98 _current_face_type == FaceInfo::VarFaceNeighbors::NEIGHBOR)
99 {
100 mooseAssert(_current_face_info->boundaryIDs().size() == 1,
101 "We should only have one boundary on every face.");
102
103 LinearFVBoundaryCondition * bc_pointer =
104 _var.getBoundaryCondition(*_current_face_info->boundaryIDs().begin());
105
106 if (bc_pointer || _force_boundary_execution)
107 {
108 if (bc_pointer)
109 bc_pointer->setupFaceData(_current_face_info, _current_face_type);
110 const auto matrix_contribution = computeBoundaryMatrixContribution(*bc_pointer);
111
112 // We allow internal (for the mesh) boundaries too, so we have to check on which side we
113 // are on (assuming that this is a boundary for the variable)
114 if ((_current_face_type == FaceInfo::VarFaceNeighbors::ELEM) && !(bounded_elem))
115 {
116 const auto dof_id_elem = _current_face_info->elemInfo()->dofIndices()[_sys_num][_var_num];
117 (*_linear_system.matrix).add(dof_id_elem, dof_id_elem, matrix_contribution);
118 }
119 else if ((_current_face_type == FaceInfo::VarFaceNeighbors::NEIGHBOR) && !(bounded_neigh))
120 {
121 const auto dof_id_neighbor =
122 _current_face_info->neighborInfo()->dofIndices()[_sys_num][_var_num];
123 (*_linear_system.matrix).add(dof_id_neighbor, dof_id_neighbor, matrix_contribution);
124 }
125 }
126 }
127}
LinearFVFluxKernel::hasBlocks
bool hasBlocks(const SubdomainName &name) const
ElemInfo::subdomain_id
SubdomainID subdomain_id() const
ElemInfo::dofIndices
const std::vector< std::vector< dof_id_type > > & dofIndices() const
FaceInfo::boundaryIDs
const std::set< BoundaryID > & boundaryIDs() const
FaceInfo::neighborPtr
const Elem * neighborPtr() const
FaceInfo::elemInfo
const ElemInfo * elemInfo() const
FaceInfo::neighborInfo
const ElemInfo * neighborInfo() const
FaceInfo::elemPtr
const Elem * elemPtr() const
LinearFVBoundaryCondition::setupFaceData
void setupFaceData(const FaceInfo *face_info, const FaceInfo::VarFaceNeighbors face_type)
LinearFVFluxKernel::_current_face_type
FaceInfo::VarFaceNeighbors _current_face_type
LinearFVFluxKernel::_matrix_contribution
DenseMatrix< Real > _matrix_contribution
LinearFVFluxKernel::_dof_indices
DenseVector< dof_id_type > _dof_indices
LinearFVFluxKernel::_force_boundary_execution
const bool _force_boundary_execution
LinearFVFluxKernel::_current_face_info
const FaceInfo * _current_face_info
LinearFVFluxKernel::_sys_num
const unsigned int _sys_num
LinearFVFluxKernel::_var_num
const unsigned int _var_num
LinearFVTurbulentAdvection::computeNeighborMatrixContribution
virtual Real computeNeighborMatrixContribution() override
Definition LinearFVTurbulentAdvection.C:201
LinearFVTurbulentAdvection::computeBoundaryMatrixContribution
virtual Real computeBoundaryMatrixContribution(const LinearFVBoundaryCondition &bc) override
Definition LinearFVTurbulentAdvection.C:220
LinearFVTurbulentAdvection::computeElemMatrixContribution
virtual Real computeElemMatrixContribution() override
Definition LinearFVTurbulentAdvection.C:194
LinearFVTurbulentAdvection::_wall_bounded
std::unordered_set< const Elem * > _wall_bounded
List for wall bounded elements.
Definition LinearFVTurbulentAdvection.h:71
LinearFVFluxKernel::_linear_system
libMesh::LinearImplicitSystem & _linear_system
MooseLinearVariableFV::getBoundaryCondition
LinearFVBoundaryCondition * getBoundaryCondition(const BoundaryID bd_id) const
libMesh::DenseMatrix::get_values
std::vector< T > & get_values()
libMesh::LinearImplicitSystem::matrix
SparseMatrix< Number > * matrix

◆ addRightHandSideContribution()

void LinearFVTurbulentAdvection::addRightHandSideContribution ( )
overridevirtual

Reimplemented from LinearFVFluxKernel.

Definition at line 130 of file LinearFVTurbulentAdvection.C.

131{
132 // Coumputing bounding map
133 const Elem * elem = _current_face_info->elemPtr();
134 const auto bounded_elem = _wall_bounded.find(elem) != _wall_bounded.end();
135 const Elem * neighbor = _current_face_info->neighborPtr();
136 const auto bounded_neigh = _wall_bounded.find(neighbor) != _wall_bounded.end();
137
138 // If we are on an internal face, we populate the two entries in the right hand side
139 // which touch the face
140 if (_current_face_type == FaceInfo::VarFaceNeighbors::BOTH)
141 {
142 // The dof ids of the variable corresponding to the element and neighbor
143 _dof_indices(0) = _current_face_info->elemInfo()->dofIndices()[_sys_num][_var_num];
144 _dof_indices(1) = _current_face_info->neighborInfo()->dofIndices()[_sys_num][_var_num];
145
146 // Compute the entries which will go to the neighbor and element positions.
147 const auto elem_rhs_contribution = computeElemRightHandSideContribution();
148 const auto neighbor_rhs_contribution = computeNeighborRightHandSideContribution();
149
150 // Populate right hand side
151 if (hasBlocks(_current_face_info->elemInfo()->subdomain_id()) && !(bounded_elem))
152 _rhs_contribution(0) = elem_rhs_contribution;
153 if (hasBlocks(_current_face_info->neighborInfo()->subdomain_id()) && !(bounded_neigh))
154 _rhs_contribution(1) = neighbor_rhs_contribution;
155
156 (*_linear_system.rhs)
157 .add_vector(_rhs_contribution.get_values().data(), _dof_indices.get_values());
158 }
159 // We are at a block boundary where the variable is not defined on one of the adjacent cells.
160 // We check if we have a boundary condition here
161 else if (_current_face_type == FaceInfo::VarFaceNeighbors::ELEM ||
162 _current_face_type == FaceInfo::VarFaceNeighbors::NEIGHBOR)
163 {
164 mooseAssert(_current_face_info->boundaryIDs().size() == 1,
165 "We should only have one boundary on every face.");
166 LinearFVBoundaryCondition * bc_pointer =
167 _var.getBoundaryCondition(*_current_face_info->boundaryIDs().begin());
168
169 if (bc_pointer || _force_boundary_execution)
170 {
171 if (bc_pointer)
172 bc_pointer->setupFaceData(_current_face_info, _current_face_type);
173
174 const auto rhs_contribution = computeBoundaryRHSContribution(*bc_pointer);
175
176 // We allow internal (for the mesh) boundaries too, so we have to check on which side we
177 // are on (assuming that this is a boundary for the variable)
178 if ((_current_face_type == FaceInfo::VarFaceNeighbors::ELEM) && !(bounded_elem))
179 {
180 const auto dof_id_elem = _current_face_info->elemInfo()->dofIndices()[_sys_num][_var_num];
181 (*_linear_system.rhs).add(dof_id_elem, rhs_contribution);
182 }
183 else if ((_current_face_type == FaceInfo::VarFaceNeighbors::NEIGHBOR) && !(bounded_neigh))
184 {
185 const auto dof_id_neighbor =
186 _current_face_info->neighborInfo()->dofIndices()[_sys_num][_var_num];
187 (*_linear_system.rhs).add(dof_id_neighbor, rhs_contribution);
188 }
189 }
190 }
191}
LinearFVFluxKernel::_rhs_contribution
DenseVector< Real > _rhs_contribution
LinearFVTurbulentAdvection::computeNeighborRightHandSideContribution
virtual Real computeNeighborRightHandSideContribution() override
Definition LinearFVTurbulentAdvection.C:214
LinearFVTurbulentAdvection::computeElemRightHandSideContribution
virtual Real computeElemRightHandSideContribution() override
Definition LinearFVTurbulentAdvection.C:208
LinearFVTurbulentAdvection::computeBoundaryRHSContribution
virtual Real computeBoundaryRHSContribution(const LinearFVBoundaryCondition &bc) override
Definition LinearFVTurbulentAdvection.C:234
libMesh::LinearImplicitSystem::rhs
NumericVector< Number > * rhs

◆ computeBoundaryMatrixContribution()

Real LinearFVTurbulentAdvection::computeBoundaryMatrixContribution ( const LinearFVBoundaryCondition &  bc)
overridevirtual

Implements LinearFVFluxKernel.

Definition at line 220 of file LinearFVTurbulentAdvection.C.

221{
222 const auto * const adv_bc = static_cast<const LinearFVAdvectionDiffusionBC *>(&bc);
223 mooseAssert(adv_bc, "This should be a valid BC!");
224
225 const auto boundary_value_matrix_contrib = adv_bc->computeBoundaryValueMatrixContribution();
226
227 // We support internal boundaries too so we have to make sure the normal points always outward
228 const auto factor = (_current_face_type == FaceInfo::VarFaceNeighbors::ELEM) ? 1.0 : -1.0;
229
230 return boundary_value_matrix_contrib * factor * _mass_face_flux * _current_face_area;
231}

Referenced by addMatrixContribution().

◆ computeBoundaryRHSContribution()

Real LinearFVTurbulentAdvection::computeBoundaryRHSContribution ( const LinearFVBoundaryCondition &  bc)
overridevirtual

Implements LinearFVFluxKernel.

Definition at line 234 of file LinearFVTurbulentAdvection.C.

235{
236 const auto * const adv_bc = static_cast<const LinearFVAdvectionDiffusionBC *>(&bc);
237 mooseAssert(adv_bc, "This should be a valid BC!");
238
239 // We support internal boundaries too so we have to make sure the normal points always outward
240 const auto factor = (_current_face_type == FaceInfo::VarFaceNeighbors::ELEM ? 1.0 : -1.0);
241
242 const auto boundary_value_rhs_contrib = adv_bc->computeBoundaryValueRHSContribution();
243 return -boundary_value_rhs_contrib * factor * _mass_face_flux * _current_face_area;
244}

Referenced by addRightHandSideContribution().

◆ computeElemMatrixContribution()

Real LinearFVTurbulentAdvection::computeElemMatrixContribution ( )
overridevirtual

Implements LinearFVFluxKernel.

Definition at line 194 of file LinearFVTurbulentAdvection.C.

195{
196 const auto & coeffs = _adv_interp_result.weights_matrix;
197 return coeffs.first * _mass_face_flux * _current_face_area;
198}
LinearFVTurbulentAdvection::_adv_interp_result
FVAdvectedInterpolationMethod::AdvectedSystemContribution _adv_interp_result
Cached weights/correction for the current face (refreshed in setupFaceData)
Definition LinearFVTurbulentAdvection.h:57

Referenced by addMatrixContribution().

◆ computeElemRightHandSideContribution()

Real LinearFVTurbulentAdvection::computeElemRightHandSideContribution ( )
overridevirtual

Implements LinearFVFluxKernel.

Definition at line 208 of file LinearFVTurbulentAdvection.C.

209{
210 return _adv_interp_result.rhs_face_value * _mass_face_flux * _current_face_area;
211}

Referenced by addRightHandSideContribution().

◆ computeNeighborMatrixContribution()

Real LinearFVTurbulentAdvection::computeNeighborMatrixContribution ( )
overridevirtual

Implements LinearFVFluxKernel.

Definition at line 201 of file LinearFVTurbulentAdvection.C.

202{
203 const auto & coeffs = _adv_interp_result.weights_matrix;
204 return coeffs.second * _mass_face_flux * _current_face_area;
205}

Referenced by addMatrixContribution().

◆ computeNeighborRightHandSideContribution()

Real LinearFVTurbulentAdvection::computeNeighborRightHandSideContribution ( )
overridevirtual

Implements LinearFVFluxKernel.

Definition at line 214 of file LinearFVTurbulentAdvection.C.

215{
216 return -_adv_interp_result.rhs_face_value * _mass_face_flux * _current_face_area;
217}

Referenced by addRightHandSideContribution().

◆ initialSetup()

void LinearFVTurbulentAdvection::initialSetup ( )
overridevirtual

Reimplemented from LinearFVFluxKernel.

Definition at line 52 of file LinearFVTurbulentAdvection.C.

53{
54 LinearFVFluxKernel::initialSetup();
55 NS::getWallBoundedElements(
56 _wall_boundary_names, _fe_problem, _subproblem, blockIDs(), _wall_bounded);
57}
LinearFVFluxKernel::blockIDs
virtual const std::set< SubdomainID > & blockIDs() const
LinearFVFluxKernel::_fe_problem
FEProblemBase & _fe_problem
LinearFVFluxKernel::initialSetup
virtual void initialSetup()
LinearFVFluxKernel::_subproblem
SubProblem & _subproblem
NS::getWallBoundedElements
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.
Definition NavierStokesMethods.C:282

◆ setupFaceData()

void LinearFVTurbulentAdvection::setupFaceData ( const FaceInfo *  face_info)
overridevirtual

Reimplemented from LinearFVFluxKernel.

Definition at line 247 of file LinearFVTurbulentAdvection.C.

248{
249 LinearFVFluxKernel::setupFaceData(face_info);
250
251 // Caching the mass flux on the face which will be reused in the advection term's matrix and right
252 // hand side contributions
253 _mass_face_flux = _mass_flux_provider.getMassFlux(*face_info);
254
255 // Only internal faces need advected interpolation results; boundary contributions are handled
256 // through the linear FV boundary conditions.
257 if (_current_face_type != FaceInfo::VarFaceNeighbors::BOTH)
258 return;
259
260 const auto state = determineState();
261 const auto & elem_info = *_current_face_info->elemInfo();
262 const auto & neighbor_info = *_current_face_info->neighborInfo();
263
264 const Real elem_value = _var.getElemValue(elem_info, state);
265 const Real neighbor_value = _var.getElemValue(neighbor_info, state);
266 if (_adv_interp_method.needsGradients())
267 {
268 const auto limiter_type = _adv_interp_method.gradientLimiter();
269 _elem_grad_storage = _var.gradSln(elem_info, state, limiter_type);
270 _neighbor_grad_storage = _var.gradSln(neighbor_info, state, limiter_type);
271 }
272
273 _adv_interp_result = _adv_interp_method.advectedInterpolate(*_current_face_info,
274 elem_value,
275 neighbor_value,
276 &_elem_grad_storage,
277 &_neighbor_grad_storage,
278 _mass_face_flux);
279}
FVAdvectedInterpolationMethod::advectedInterpolate
virtual AdvectedSystemContribution advectedInterpolate(const FaceInfo &face, Real elem_value, Real neighbor_value, const VectorValue< Real > *elem_grad, const VectorValue< Real > *neighbor_grad, Real mass_flux) const=0
LinearFVFluxKernel::setupFaceData
virtual void setupFaceData(const FaceInfo *face_info)
LinearFVTurbulentAdvection::_elem_grad_storage
VectorValue< Real > _elem_grad_storage
Reusable gradient storage used when advected interpolation requires gradients.
Definition LinearFVTurbulentAdvection.h:60
MooseLinearVariableFV::getElemValue
Real getElemValue(const ElemInfo &elem_info, const StateArg &state) const
MooseLinearVariableFV::gradSln
VectorValue< Real > gradSln(const ElemInfo &elem_info, const StateArg &state) const
RhieChowMassFlux::getMassFlux
Real getMassFlux(const FaceInfo &fi) const
Get the face velocity times density (used in advection terms)
Definition RhieChowMassFlux.C:290
LinearFVFluxKernel::determineState
Moose::StateArg determineState() const
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ validParams()

InputParameters LinearFVTurbulentAdvection::validParams ( )
static

Definition at line 18 of file LinearFVTurbulentAdvection.C.

19{
20 InputParameters params = LinearFVFluxKernel::validParams();
21 params.addClassDescription("Represents the matrix and right hand side contributions of an "
22 "advection term for a turbulence variable.");
23
24 params.addRequiredParam<UserObjectName>(
25 "rhie_chow_user_object",
26 "The rhie-chow user-object which is used to determine the face velocity.");
27
28 params.addRequiredParam<InterpolationMethodName>(
29 "advected_interp_method_name",
30 "Name of the FVInterpolationMethod to use for the advected quantity.");
31
32 params.addParam<std::vector<BoundaryName>>(
33 "walls", {}, "Boundaries that correspond to solid walls.");
34
35 return params;
36}
InputParameters::addRequiredParam
void addRequiredParam(const std::string &name, const std::string &doc_string)
InputParameters::addParam
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
InputParameters::addClassDescription
void addClassDescription(const std::string &doc_string)
LinearFVFluxKernel::validParams
static InputParameters validParams()

Member Data Documentation

◆ _adv_interp_method

const FVAdvectedInterpolationMethod& LinearFVTurbulentAdvection::_adv_interp_method
private

The interpolation method to use for the advected quantity.

Definition at line 54 of file LinearFVTurbulentAdvection.h.

Referenced by LinearFVTurbulentAdvection(), and setupFaceData().

◆ _adv_interp_result

FVAdvectedInterpolationMethod::AdvectedSystemContribution LinearFVTurbulentAdvection::_adv_interp_result
private

Cached weights/correction for the current face (refreshed in setupFaceData)

Definition at line 57 of file LinearFVTurbulentAdvection.h.

Referenced by computeElemMatrixContribution(), computeElemRightHandSideContribution(), computeNeighborMatrixContribution(), computeNeighborRightHandSideContribution(), and setupFaceData().

◆ _elem_grad_storage

VectorValue<Real> LinearFVTurbulentAdvection::_elem_grad_storage
private

Reusable gradient storage used when advected interpolation requires gradients.

Definition at line 60 of file LinearFVTurbulentAdvection.h.

Referenced by setupFaceData().

◆ _mass_face_flux

Real LinearFVTurbulentAdvection::_mass_face_flux
private

Container for the mass flux on the face which will be reused in the advection term's matrix and right hand side contribution.

Definition at line 65 of file LinearFVTurbulentAdvection.h.

Referenced by computeBoundaryMatrixContribution(), computeBoundaryRHSContribution(), computeElemMatrixContribution(), computeElemRightHandSideContribution(), computeNeighborMatrixContribution(), computeNeighborRightHandSideContribution(), and setupFaceData().

◆ _mass_flux_provider

const RhieChowMassFlux& LinearFVTurbulentAdvection::_mass_flux_provider
protected

The Rhie-Chow user object that provides us with the face velocity.

Definition at line 50 of file LinearFVTurbulentAdvection.h.

Referenced by setupFaceData().

◆ _neighbor_grad_storage

VectorValue<Real> LinearFVTurbulentAdvection::_neighbor_grad_storage
private

Definition at line 61 of file LinearFVTurbulentAdvection.h.

Referenced by setupFaceData().

◆ _wall_boundary_names

const std::vector<BoundaryName>& LinearFVTurbulentAdvection::_wall_boundary_names
private

Wall boundaries.

Definition at line 68 of file LinearFVTurbulentAdvection.h.

Referenced by initialSetup().

◆ _wall_bounded

std::unordered_set<const Elem *> LinearFVTurbulentAdvection::_wall_bounded
private

List for wall bounded elements.

Definition at line 71 of file LinearFVTurbulentAdvection.h.

Referenced by addMatrixContribution(), addRightHandSideContribution(), and initialSetup().

The documentation for this class was generated from the following files:
Generated on Fri Aug 21 2026 13:43:12 for https://mooseframework.inl.gov by doxygen 1.9.8