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 | Static Protected Member Functions | Protected Attributes | Private Member Functions | Static Private Member Functions | Private Attributes | List of all members
WCNSLinearFVFlowPhysics Class Referencefinal

Creates all the objects needed to solve the Navier-Stokes equations with the SIMPLE algorithm using the linear finite volume discretization Currently does not implement: More...

#include <WCNSLinearFVFlowPhysics.h>

Inheritance diagram for WCNSLinearFVFlowPhysics:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 WCNSLinearFVFlowPhysics (const InputParameters &parameters)
 
virtual std::vector< UserObjectName > getSuppliedUserObjects () const override
 
bool hasFlowEquations () const
 Whether the physics is actually creating the flow equations.
 
bool addAxisymmetricViscousSourceEnabled () const
 Whether the cylindrical viscous source helper is enabled.
 
bool includeSymmetrizedViscousStress () const
 Whether to include the symmetrized contribution in the viscous stress.
 
bool includeIsotropicStress () const
 Whether to include the isotropic viscous stress contribution.
 
const std::vector< std::string > & getVelocityNames () const
 To interface with other Physics.
 
const NonlinearVariableName & getPressureName () const
 
const NonlinearVariableName & getFluidTemperatureName () const
 
MooseFunctorName getPorosityFunctorName (const bool smoothed) const
 
const MooseEnum & compressibility () const
 Return the compressibility of the flow equations selected.
 
bool porousMediumTreatment () const
 Return whether a porous medium treatment is applied.
 
RealVectorValue gravityVector () const
 Return the gravity vector.
 
const MooseFunctorName & densityName () const
 Return the name of the density functor.
 
const MooseFunctorName & dynamicViscosityName () const
 Return the name of the dynamic viscosity functor.
 
const MooseEnum & getVelocityFaceInterpolationMethod () const
 Get the face interpolation method for velocity.
 
const MooseEnum & getMomentumAdvectionFaceInterpolationMethod () const
 Get the face interpolation method for momentum in the advection term.
 
const MooseEnum & getMomentumFaceInterpolationMethod () const
 Get the face interpolation method for momentum (mostly used in the stress terms)
 
const std::vector< BoundaryName > & getInletBoundaries () const
 Get the inlet boundaries.
 
const std::vector< BoundaryName > & getOutletBoundaries () const
 Get the outlet boundaries.
 
const std::vector< BoundaryName > & getWallBoundaries () const
 Get the wall boundaries.
 
const std::vector< BoundaryName > & getHydraulicSeparators () const
 Get the hydraulic separator boundaries.
 
NS::MomentumInletTypes inletBoundaryType (const BoundaryName &boundary_name) const
 Get the type of the inlet BC.
 
const std::vector< Point > & getFluxInletDirections () const
 Get the inlet direction if using a flux inlet.
 
const std::vector< PostprocessorName > & getFluxInletPPs () const
 Get the inlet flux postprocessor if using a flux inlet.
 
const UserObjectName & rhieChowUOName () const
 Return the name of the Rhie Chow user object.
 
virtual void act () override final
 
void addBlocks (const std::vector< SubdomainName > &blocks)
 
void addBlocksById (const std::vector< SubdomainID > &block_ids)
 
const std::vector< SubdomainName > & blocks () const
 
bool checkBlockRestrictionIdentical (const std::string &object_name, const std::vector< SubdomainName > &blocks, const bool error_if_not_identical=true) const
 
bool hasBlocks (const std::vector< SubdomainName > &blocks) const
 
const T * getCoupledPhysics (const PhysicsName &phys_name, const bool allow_fail=false) const
 
const std::vector< T * > getCoupledPhysics (const bool allow_fail=false) const
 
unsigned int dimension () const
 
const ActionComponent & getActionComponent (const ComponentName &comp_name) const
 
void checkComponentType (const ActionComponent &component) const
 
virtual void addComponent (const ActionComponent &component)
 
const std::vector< VariableName > & solverVariableNames () const
 
const std::vector< VariableName > & auxVariableNames () const
 
void timedAct ()
 
MooseObjectName uniqueActionName () const
 
const std::string & specificTaskName () const
 
const std::set< std::string > & getAllTasks () const
 
void appendTask (const std::string &task)
 
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
 
PerfGraph & perfGraph ()
 
void assertParamDefined (const std::string &libmesh_dbg_var(param)) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 

Static Public Member Functions

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

Public Attributes

 usingCombinedWarningSolutionWarnings
 
const ConsoleStream _console
 

Static Public Attributes

static const std::string unique_action_name_param
 
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
 
static constexpr auto SYSTEM
 
static constexpr auto NAME
 

Protected Member Functions

virtual void initializePhysicsAdditional () override
 
virtual void actOnAdditionalTasks () override
 
virtual void addInitialConditions () override
 
virtual void addFVBCs () override
 
virtual void addMaterials () override
 
virtual void addPostprocessors () override
 
void addAxisymmetricViscousSource ()
 Adds the cylindrical source kernel for the radial momentum equation when requested and valid.
 
virtual void addAxisymmetricViscousSourceKernel (const std::vector< SubdomainName > &, unsigned int)
 Derived classes must override this hook to add the actual object that implements the axisymmetric viscous source term for their formulation once the helper has determined the relevant blocks and radial component.
 
void addPorousMediumSpeedMaterial ()
 Add material to define the local speed in porous medium flows.
 
void addNonPorousMediumSpeedMaterial ()
 Add material to define the local speed with no porous medium treatment.
 
void addFluidPropertiesFunctorMaterial ()
 Function which adds the general functor fluid properties functor material to define fluid functor material property.
 
VariableName getFlowVariableName (const std::string &default_name) const
 Convenience routine to be able to retrieve the actual variable names from their default names.
 
bool hasTurbulencePhysics () const
 Whether a turbulence Physics has been coupled in, to know which viscosity to pick on symmetry boundary conditions.
 
const WCNSFVTurbulencePhysicsBase * getCoupledTurbulencePhysics () const
 Find the turbulence physics.
 
std::vector< SubdomainName > getAxisymmetricRZBlocks () const
 Return the set of blocks restricted to an RZ coordinate system.
 
bool usingNavierStokesFVSyntax () const
 Detects if we are using the new Physics syntax or the old NavierStokesFV action.
 
InputParameters getAdditionalRMParams () const override
 Parameters to change or add relationship managers.
 
void addFVAdvectedInterpolationMethod (const MooseEnum &interpolation_method)
 Add the FVInterpolationMethod object for an advected interpolation method if absent.
 
void assertParamDefined (const std::string &param) const
 
bool isTransient () const
 
Factory & getFactory ()
 
Factory & getFactory () const
 
virtual FEProblemBase & getProblem ()
 
virtual const FEProblemBase & getProblem () const
 
void prepareCopyVariablesFromMesh () const
 
void copyVariablesFromMesh (const std::vector< VariableName > &variables_to_copy, bool are_nonlinear=true)
 
std::string prefix () const
 
void addUserObject (const std::string &uo_type, const std::string &uo_name, InputParameters &params)
 
void saveSolverVariableName (const VariableName &var_name)
 
void saveAuxVariableName (const VariableName &var_name)
 
void addSharedSolverVariable (const std::string &var_type, const VariableName &var_name, InputParameters params)
 
void addSharedSolverVariable (const VariableName &var_name, const std::string &family, const std::string &order, Real scaling=1.0)
 
void addSharedAuxVariable (const std::string &var_type, const VariableName &var_name, const InputParameters &params)
 
void addSharedAuxVariable (const VariableName &var_name, const std::string &family, const std::string &order)
 
bool variableExists (const VariableName &var_name, bool error_if_aux) const
 
bool solverVariableExists (const VariableName &var_name) const
 
const SolverSystemName & getSolverSystem (unsigned int variable_index) const
 
const SolverSystemName & getSolverSystem (const VariableName &variable_name) const
 
void addRequiredPhysicsTask (const std::string &task)
 
void assignBlocks (InputParameters &params, const std::vector< SubdomainName > &blocks) const
 
bool allMeshBlocks (const std::vector< SubdomainName > &blocks) const
 
bool allMeshBlocks (const std::set< SubdomainName > &blocks) const
 
std::set< SubdomainID > getSubdomainIDs (const std::set< SubdomainName > &blocks) const
 
std::vector< std::string > getSubdomainNamesAndIDs (const std::set< SubdomainID > &blocks) const
 
void addPetscPairsToPetscOptions (const std::vector< std::pair< MooseEnumItem, std::string > > &petsc_pair_options)
 
bool isVariableFV (const VariableName &var_name) const
 
bool isVariableScalar (const VariableName &var_name) const
 
bool shouldCreateVariable (const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool error_if_aux)
 
bool shouldCreateIC (const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool ic_is_default_ic, const bool error_if_already_defined) const
 
bool shouldCreateTimeDerivative (const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool error_if_already_defined) const
 
void reportPotentiallyMissedParameters (const std::vector< std::string > &param_names, const std::string &object_type, const std::string &object_name="") const
 
virtual void checkIntegrity () const
 
void associateWithParameter (const std::string &param_name, InputParameters &params) const
 
void associateWithParameter (const InputParameters &from_params, const std::string &param_name, InputParameters &params) const
 
const T & getMeshProperty (const std::string &data_name, const std::string &prefix)
 
const T & getMeshProperty (const std::string &data_name)
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name, const std::string &prefix) const
 
bool hasMeshProperty (const std::string &data_name) const
 
bool hasMeshProperty (const std::string &data_name) const
 
std::string meshPropertyName (const std::string &data_name) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level) const
 
PerfID registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const
 
std::string timedSectionName (const std::string &section_name) const
 
void flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const
 
InvalidSolutionID registerInvalidSolutionInternal (const std::string &message, const bool warning) const
 
void checkParamsBothSetOrNotSet (const std::string &param1, const std::string &param2) const
 
void checkAtMostOneParamSetByUser (const std::vector< std::string > &params) const
 
void checkSecondParamSetOnlyIfFirstOneTrue (const std::string &param1, const std::string &param2) const
 
void checkSecondParamSetOnlyIfFirstOneSet (const std::string &param1, const std::string &param2) const
 
void checkSecondParamNotSetIfFirstOneSet (const std::string &param1, const std::string &param2) const
 
void checkVectorParamsSameLength (const std::string &param1, const std::string &param2) const
 
void checkVectorParamAndMultiMooseEnumLength (const std::string &param1, const std::string &param2) const
 
void checkTwoDVectorParamsSameLength (const std::string &param1, const std::string &param2) const
 
void checkVectorParamsNoOverlap (const std::vector< std::string > &param_vecs) const
 
void checkTwoDVectorParamsNoRespectiveOverlap (const std::vector< std::string > &param_vecs) const
 
void checkTwoDVectorParamInnerSameLengthAsOneDVector (const std::string &param1, const std::string &param2) const
 
void checkTwoDVectorParamMultiMooseEnumSameLength (const std::string &param1, const std::string &param2, const bool error_for_param2) const
 
void checkVectorParamNotEmpty (const std::string &param1) const
 
void checkVectorParamsSameLengthIfSet (const std::string &param1, const std::string &param2, const bool ignore_empty_default_param2=false) const
 
void checkVectorParamLengthSameAsCombinedOthers (const std::string &param1, const std::string &param2, const std::string &param3) const
 
void checkBlockwiseConsistency (const std::string &block_param_name, const std::vector< std::string > &parameter_names) const
 
bool parameterConsistent (const InputParameters &other_param, const std::string &param_name) const
 
void warnInconsistent (const InputParameters &parameters, const std::string &param_name) const
 
void errorDependentParameter (const std::string &param1, const std::string &value_not_set, const std::vector< std::string > &dependent_params) const
 
void errorInconsistentDependentParameter (const std::string &param1, const std::string &value_set, const std::vector< std::string > &dependent_params) const
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 

Protected Attributes

const TagName _pressure_tag = "p_tag"
 Name of the vector to hold pressure momentum equation contributions.
 
const bool _has_flow_equations
 Boolean to keep track of whether the flow equations should be created.
 
const bool _add_rz_viscous_source
 Whether to automatically add the cylindrical viscous source term.
 
const MooseEnum _compressibility
 Compressibility type, can be compressible, incompressible or weakly-compressible.
 
const bool _solve_for_dynamic_pressure
 Whether we are solving for the total or dynamic pressure.
 
const bool _porous_medium_treatment
 Whether to use the porous medium treatment.
 
const MooseFunctorName _porosity_name
 Name of the porosity functor.
 
MooseFunctorName _flow_porosity_functor_name
 Name of the porosity functor for the flow equations (if smoothed)
 
const std::vector< std::string > _velocity_names
 Velocity names.
 
const NonlinearVariableName _pressure_name
 Pressure name.
 
const NonlinearVariableName _fluid_temperature_name
 Fluid temperature name.
 
const MooseFunctorName _density_name
 Name of the density material property.
 
const MooseFunctorName _density_gravity_name
 Name of the density material property used for gravity and Boussinesq terms.
 
const MooseFunctorName _dynamic_viscosity_name
 Name of the dynamic viscosity material property.
 
const bool _include_symmetrized_viscous_stress
 Whether to include the symmetrized viscous stress contribution.
 
const bool _include_isotropic_viscous_stress
 Whether to include the isotropic viscous stress contribution.
 
UserObjectName _rc_uo_name
 name of the Rhie Chow user object
 
const MooseEnum _velocity_interpolation
 The velocity face interpolation method for advecting other quantities.
 
const MooseEnum _momentum_advection_interpolation
 The momentum face interpolation method for being advected.
 
const MooseEnum _momentum_face_interpolation
 The momentum face interpolation method for stress terms.
 
const WCNSFVTurbulencePhysicsBase * _turbulence_physics
 Can be set to a coupled turbulence physics.
 
std::vector< std::vector< SubdomainName > > _friction_blocks
 Subdomains where we want to have volumetric friction.
 
std::vector< std::vector< std::string > > _friction_types
 The friction correlation types used for each block.
 
std::vector< std::vector< std::string > > _friction_coeffs
 The coefficients used for each item if friction type.
 
const std::vector< BoundaryName > _inlet_boundaries
 Boundaries with a flow inlet specified on them.
 
const std::vector< BoundaryName > _outlet_boundaries
 Boundaries with a flow outlet specified on them.
 
const std::vector< BoundaryName > _wall_boundaries
 Boundaries which define a wall (slip/noslip/etc.)
 
const std::vector< BoundaryName > _hydraulic_separators
 Hydraulic separator boundaries.
 
std::map< BoundaryName, MooseEnum > _momentum_inlet_types
 Momentum inlet boundary types.
 
std::map< BoundaryName, MooseEnum > _momentum_outlet_types
 Momentum outlet boundary types.
 
std::map< BoundaryName, MooseEnum > _momentum_wall_types
 Momentum wall boundary types.
 
std::vector< PostprocessorName > _flux_inlet_pps
 Postprocessors describing the momentum inlet for each boundary. Indexing based on the number of flux boundaries.
 
std::vector< Point > _flux_inlet_directions
 Direction of each flux inlet. Indexing based on the number of flux boundaries.
 
std::map< BoundaryName, std::vector< MooseFunctorName > > _momentum_inlet_functors
 Functors describing the momentum inlet for each boundary.
 
std::map< BoundaryName, MooseFunctorName > _pressure_functors
 Functors describing the outlet pressure on each boundary.
 
std::map< BoundaryName, std::vector< MooseFunctorName > > _momentum_wall_functors
 Functors describing the momentum for each wall boundary.
 
bool _define_variables
 Whether to define variables if they do not exist.
 
std::vector< SolverSystemName > _system_names
 
std::vector< unsigned int > _system_numbers
 
const bool _verbose
 
const MooseEnum & _preconditioning
 
std::vector< SubdomainName > _blocks
 
std::string _registered_identifier
 
std::string _specific_task_name
 
std::set< std::string > _all_tasks
 
ActionWarehouse & _awh
 
const std::string & _current_task
 
std::shared_ptr< MooseMesh > & _mesh
 
std::shared_ptr< MooseMesh > & _displaced_mesh
 
std::shared_ptr< FEProblemBase > & _problem
 
PerfID _act_timer
 
MooseApp & _app
 
Factory & _factory
 
ActionFactory & _action_factory
 
const std::string & _type
 
const std::string & _name
 
const InputParameters & _pars
 
MooseApp & _pg_moose_app
 
const std::string _prefix
 
const Parallel::Communicator & _communicator
 

Private Member Functions

virtual void addSolverVariables () override
 
virtual void addFVInterpolationMethods () override
 
virtual void addFVKernels () override
 
virtual void addUserObjects () override
 
void addPressureCorrectionKernels ()
 Function adding kernels for the incompressible pressure correction equation.
 
void addMomentumTimeKernels () override
 Functions adding kernels for the incompressible momentum equation If the material properties are not constant, these can be used for weakly-compressible simulations (except the Boussinesq kernel) as well.
 
void addMomentumFluxKernels ()
 
virtual void addMomentumPressureKernels () override
 
virtual void addMomentumGravityKernels () override
 
virtual void addMomentumFrictionKernels () override
 
virtual void addMomentumBoussinesqKernels () override
 
virtual void addInletBC () override
 Functions adding boundary conditions for the flow simulation.
 
virtual void addOutletBC () override
 
virtual void addWallsBC () override
 
virtual void addSeparatorBC () override
 
virtual bool hasForchheimerFriction () const override
 Return whether a Forchheimer friction model is in use.
 
virtual void addRhieChowUserObjects () override
 Function which adds the RhieChow interpolator user objects for weakly and incompressible formulations.
 
virtual void addFunctorMaterials () override
 
virtual MooseFunctorName getLinearFrictionCoefName () const override
 Get the name of the linear friction coefficient. Returns an empty string if no friction.
 
unsigned short getNumberAlgebraicGhostingLayersNeeded () const override
 Return the number of algebraic ghosting layers needed.
 
virtual void addRelationshipManagers (Moose::RelationshipManagerType input_rm_type) override
 
bool addRelationshipManagers (Moose::RelationshipManagerType when_type, const InputParameters &moose_object_pars)
 
virtual void addRelationshipManagers (Moose::RelationshipManagerType when_type)
 
void initializePhysics ()
 
virtual void checkIntegrityEarly () const
 
virtual void addAuxiliaryVariables ()
 
virtual void addFEKernels ()
 
virtual void addNodalKernels ()
 
virtual void addDiracKernels ()
 
virtual void addDGKernels ()
 
virtual void addScalarKernels ()
 
virtual void addInterfaceKernels ()
 
virtual void addFVInterfaceKernels ()
 
virtual void addFEBCs ()
 
virtual void addNodalBCs ()
 
virtual void addPeriodicBCs ()
 
virtual void addFunctions ()
 
virtual void addAuxiliaryKernels ()
 
virtual void addCorrectors ()
 
virtual void addMultiApps ()
 
virtual void addTransfers ()
 
virtual void addVectorPostprocessors ()
 
virtual void addReporters ()
 
virtual void addOutputs ()
 
virtual void addPreconditioning ()
 
virtual void addExecutioner ()
 
virtual void addExecutors ()
 
void checkRequiredTasks () const
 
bool addRelationshipManager (Moose::RelationshipManagerType input_rm_type, const InputParameters &moose_object_pars, std::string rm_name, Moose::RelationshipManagerType rm_type, Moose::RelationshipManagerInputParameterCallback rm_input_parameter_func, Moose::RMSystemType sys_type=Moose::RMSystemType::NONE)
 
virtual std::string meshPropertyPrefix (const std::string &data_name) const
 
const RestartableDataValue & getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const
 
void mooseErrorInternal (Args &&... args) const
 
const T & forwardGetParam (const std::string &param_name) const
 
const InputParameters & forwardParameters () const
 
bool forwardIsParamSetByUser (const std::string &param_name) const
 
bool forwardIsParamValid (const std::string &param_name) const
 
void forwardParamError (Args &&... args) const
 
void forwardMooseError (Args &&... args) const
 
void forwardMooseWarning (Args &&... args) const
 
const std::string & forwardType () const
 
virtual const std::string & forwardName () const
 
const std::vector< SubdomainName > & forwardBlocks () const
 

Static Private Member Functions

static const hit::Node * getHitNode (const InputParameters &params)
 
static std::string messagePrefix (const InputParameters &params, const bool hit_prefix)
 

Private Attributes

const bool _non_orthogonal_correction
 Whether to use the correction term for non-orthogonality.
 
MooseEnum _is_transient
 
std::vector< VariableName > _solver_var_names
 
std::vector< VariableName > _aux_var_names
 
unsigned int _dim
 
std::set< std::string > _required_tasks
 
const ParallelParamObject & _parent
 
MooseApp & _meta_data_app
 
const MooseObject *const _meta_data_object
 
const MooseBase & _si_moose_base
 
const FEProblemBase * _si_problem
 
const PhysicsBase *const _customer_class
 

Detailed Description

Creates all the objects needed to solve the Navier-Stokes equations with the SIMPLE algorithm using the linear finite volume discretization Currently does not implement:

Definition at line 24 of file WCNSLinearFVFlowPhysics.h.

Constructor & Destructor Documentation

◆ WCNSLinearFVFlowPhysics()

WCNSLinearFVFlowPhysics::WCNSLinearFVFlowPhysics ( const InputParameters &  parameters)

Definition at line 80 of file WCNSLinearFVFlowPhysics.C.

82 _non_orthogonal_correction(getParam<bool>("orthogonality_correction"))
83{
84 addRequiredPhysicsTask("add_interpolation_method_physics");
85
87 paramError("porous_medium_treatment", "Porous media unsupported");
89 mooseError("Not supported");
90
91 if (_hydraulic_separators.size())
92 paramError("hydraulic_separator_sidesets",
93 "Flow separators are not supported yet for linearFV!");
94 if (getParam<bool>("pin_pressure"))
95 paramError("pin_pressure",
96 "Pressure pinning is implemented in the executioner for the linear finite volume "
97 "segregated solves");
98}
const InputParameters & parameters() const
void paramError(const std::string &param, Args... args) const
void mooseError(Args &&... args) const
void addRequiredPhysicsTask(const std::string &task)
Base class for Physics which create the Navier Stokes flow equations.
const bool _has_flow_equations
Boolean to keep track of whether the flow equations should be created.
const bool _porous_medium_treatment
Whether to use the porous medium treatment.
const std::vector< BoundaryName > _hydraulic_separators
Hydraulic separator boundaries.
const bool _non_orthogonal_correction
Whether to use the correction term for non-orthogonality.

Member Function Documentation

◆ actOnAdditionalTasks()

void WCNSFVFlowPhysicsBase::actOnAdditionalTasks ( )
overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 289 of file WCNSFVFlowPhysicsBase.C.

290{
291 // Turbulence physics would not be initialized before this task
292 if (_current_task == "get_turbulence_physics")
294}
const std::string & _current_task
const WCNSFVTurbulencePhysicsBase * getCoupledTurbulencePhysics() const
Find the turbulence physics.
const WCNSFVTurbulencePhysicsBase * _turbulence_physics
Can be set to a coupled turbulence physics.

◆ addAxisymmetricViscousSource()

void WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource ( )
protectedinherited

Adds the cylindrical source kernel for the radial momentum equation when requested and valid.

Definition at line 570 of file WCNSFVFlowPhysicsBase.C.

571{
573 return;
574
575 const auto rz_blocks = getAxisymmetricRZBlocks();
576 if (rz_blocks.empty())
577 return;
578
579 const auto radial_index = getProblem().mesh().getAxisymmetricRadialCoord();
580 addAxisymmetricViscousSourceKernel(rz_blocks, radial_index);
581}
virtual MooseMesh & mesh() override
unsigned int getAxisymmetricRadialCoord() const
virtual FEProblemBase & getProblem()
virtual void addAxisymmetricViscousSourceKernel(const std::vector< SubdomainName > &, unsigned int)
Derived classes must override this hook to add the actual object that implements the axisymmetric vis...
std::vector< SubdomainName > getAxisymmetricRZBlocks() const
Return the set of blocks restricted to an RZ coordinate system.
const bool _add_rz_viscous_source
Whether to automatically add the cylindrical viscous source term.

Referenced by WCNSFVFlowPhysics::addFVKernels().

◆ addAxisymmetricViscousSourceEnabled()

bool WCNSFVFlowPhysicsBase::addAxisymmetricViscousSourceEnabled ( ) const
inlineinherited

Whether the cylindrical viscous source helper is enabled.

Definition at line 39 of file WCNSFVFlowPhysicsBase.h.

39{ return _add_rz_viscous_source; }

Referenced by WCNSFVTurbulencePhysics::addAxisymmetricTurbulentViscousSource().

◆ addAxisymmetricViscousSourceKernel()

virtual void WCNSFVFlowPhysicsBase::addAxisymmetricViscousSourceKernel ( const std::vector< SubdomainName > &  ,
unsigned int   
)
inlineprotectedvirtualinherited

Derived classes must override this hook to add the actual object that implements the axisymmetric viscous source term for their formulation once the helper has determined the relevant blocks and radial component.

Reimplemented in WCNSFVFlowPhysics.

Definition at line 127 of file WCNSFVFlowPhysicsBase.h.

129 {
130 }

Referenced by WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource().

◆ addFluidPropertiesFunctorMaterial()

void WCNSFVFlowPhysicsBase::addFluidPropertiesFunctorMaterial ( )
protectedinherited

Function which adds the general functor fluid properties functor material to define fluid functor material property.

Definition at line 354 of file WCNSFVFlowPhysicsBase.C.

355{
356 // Not very future-proof but it works
357 const bool use_ad = !dynamic_cast<WCNSLinearFVFlowPhysics *>(this);
358 const std::string class_name =
359 use_ad ? "GeneralFunctorFluidProps" : "NonADGeneralFunctorFluidProps";
360 InputParameters params = getFactory().getValidParams(class_name);
361 assignBlocks(params, _blocks);
362
363 params.set<MooseFunctorName>(NS::pressure) = _pressure_name;
364 params.set<MooseFunctorName>(NS::T_fluid) = _fluid_temperature_name;
365 params.set<MooseFunctorName>(NS::speed) = NS::speed;
366 params.applySpecificParameters(parameters(), {NS::fluid, NS::density, "mu_rampdown"});
368 params.set<bool>("force_define_density") = true;
370 {
371 params.set<MooseFunctorName>(NS::porosity) = "1";
372 params.set<MooseFunctorName>("characteristic_length") = "1";
373 }
374 else
375 // not implemented yet
376 paramInfo(
377 NS::fluid,
378 "Specifying the fluid properties user object does not define the GeneralFunctorFluidProps "
379 "when using the porous medium treatment. You have to define this object in the input");
380
381 // Dynamic pressure
382 params.set<bool>("solving_for_dynamic_pressure") = _solve_for_dynamic_pressure;
384 {
385 params.set<Point>("reference_pressure_point") = getParam<Point>("reference_pressure_point");
386 if (!isParamSetByUser("reference_pressure_point"))
387 paramWarning("reference_pressure_point",
388 "Default value of (0,0,0) used. If this point is outside the flow domain, the "
389 "simulation will error");
390 params.set<Real>("reference_pressure") = getParam<Real>("reference_pressure");
391 }
392 params.set<Point>("gravity") = getParam<RealVectorValue>("gravity");
393
395 getProblem().addFunctorMaterial(class_name, prefix() + "functor_fluidprops", params);
396}
virtual void addFunctorMaterial(const std::string &functor_material_name, const std::string &name, InputParameters &parameters)
InputParameters getValidParams(const std::string &name) const
void applySpecificParameters(const InputParameters &common, const std::vector< std::string > &include, bool allow_private=false)
T & set(const std::string &name, bool quiet_mode=false)
void paramInfo(const std::string &param, Args... args) const
void paramWarning(const std::string &param, Args... args) const
bool isParamSetByUser(const std::string &name) const
Factory & getFactory()
void assignBlocks(InputParameters &params, const std::vector< SubdomainName > &blocks) const
std::string prefix() const
std::vector< SubdomainName > _blocks
const NonlinearVariableName _pressure_name
Pressure name.
const NonlinearVariableName _fluid_temperature_name
Fluid temperature name.
const bool _solve_for_dynamic_pressure
Whether we are solving for the total or dynamic pressure.
const MooseFunctorName _density_name
Name of the density material property.
Creates all the objects needed to solve the Navier-Stokes equations with the SIMPLE algorithm using t...
bool parsesToReal(const std::string &input, Real *parsed_real)
static const std::string density
Definition NS.h:34
static const std::string T_fluid
Definition NS.h:110
static const std::string porosity
Definition NS.h:108
static const std::string speed
Definition NS.h:148
static const std::string fluid
Definition NS.h:88
static const std::string pressure
Definition NS.h:57
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

◆ addFunctorMaterials()

void WCNSLinearFVFlowPhysics::addFunctorMaterials ( )
overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 704 of file WCNSLinearFVFlowPhysics.C.

705{
706 if (parameters().isParamValid("gravity"))
707 {
708 const auto gravity_vector = getParam<RealVectorValue>("gravity");
709 const std::vector<std::string> comp_axis({"x", "y", "z"});
710 for (const auto d : make_range(dimension()))
711 if (gravity_vector(d) != 0)
712 {
713 // Add rho * g functor for each relevant direction
714 // TODO: we could avoid using an AD functor material for non-AD density functor
715 auto params = getFactory().getValidParams("ADParsedFunctorMaterial");
716 assignBlocks(params, _blocks);
717 params.set<std::string>("expression") =
718 _density_gravity_name + " * " + std::to_string(gravity_vector(d));
720 params.set<std::vector<std::string>>("functor_names") = {_density_gravity_name};
721 params.set<std::string>("property_name") = "rho_g_" + comp_axis[d];
722 // We don't output this helper material
724 "ADParsedFunctorMaterial", prefix() + "gravity_helper_" + comp_axis[d], params);
725 }
726 }
727}
virtual void addMaterial(const std::string &material_name, const std::string &name, InputParameters &parameters)
bool isParamValid(const std::string &name) const
unsigned int dimension() const
const MooseFunctorName _density_gravity_name
Name of the density material property used for gravity and Boussinesq terms.
if(subdm)
IntRange< T > make_range(T beg, T end)

◆ addFVAdvectedInterpolationMethod()

void NavierStokesPhysicsBase::addFVAdvectedInterpolationMethod ( const MooseEnum &  interpolation_method)
protectedinherited

Add the FVInterpolationMethod object for an advected interpolation method if absent.

Definition at line 43 of file NavierStokesPhysicsBase.C.

44{
45 const std::string method_name = interpolation_method;
46 if (getProblem().hasFVInterpolationMethod(method_name))
47 return;
48
49 const auto method_type = NS::fvAdvectedInterpolationMethodType(interpolation_method);
50
51 InputParameters params = getFactory().getValidParams(method_type);
52 getProblem().addFVInterpolationMethod(method_type, method_name, params);
53}
virtual void addFVInterpolationMethod(const std::string &method_type, const std::string &name, InputParameters &parameters)
std::string fvAdvectedInterpolationMethodType(const MooseEnum &interpolation_method)
Gets the FVInterpolationMethod object type for an advected interpolation method.
Definition NSFVUtils.C:67

Referenced by addFVInterpolationMethods(), WCNSLinearFVScalarTransportPhysics::addFVInterpolationMethods(), and WCNSLinearFVTurbulencePhysics::addFVInterpolationMethods().

◆ addFVBCs()

void WCNSFVFlowPhysicsBase::addFVBCs ( )
overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 297 of file WCNSFVFlowPhysicsBase.C.

298{
299 addInletBC();
300 addOutletBC();
301 addWallsBC();
303}
virtual void addInletBC()=0
Functions adding boundary conditions for the flow simulation.
virtual void addWallsBC()=0
virtual void addOutletBC()=0
virtual void addSeparatorBC()=0

◆ addFVInterpolationMethods()

void WCNSLinearFVFlowPhysics::addFVInterpolationMethods ( )
overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 101 of file WCNSLinearFVFlowPhysics.C.

102{
103 if (!_has_flow_equations || isParamValid("momentum_advection_interpolation_method_name"))
104 return;
105
106 addFVAdvectedInterpolationMethod(getParam<MooseEnum>("momentum_advection_interpolation"));
107}
void addFVAdvectedInterpolationMethod(const MooseEnum &interpolation_method)
Add the FVInterpolationMethod object for an advected interpolation method if absent.

◆ addFVKernels()

void WCNSLinearFVFlowPhysics::addFVKernels ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 186 of file WCNSLinearFVFlowPhysics.C.

187{
189 return;
190
191 // Pressure correction equation: divergence of momentum
193
194 // Momentum equation: time derivative
195 if (isTransient())
197
198 // Momentum equation: flux terms
200
201 // Momentum equation: pressure term
203
204 // Momentum equation: friction term
205 if (_friction_types.size())
207
208 // Momentum equation: gravity source term
210
211 // Momentum equation: boussinesq approximation
212 if (getParam<bool>("boussinesq_approximation"))
214}
bool isTransient() const
std::vector< std::vector< std::string > > _friction_types
The friction correlation types used for each block.
virtual void addMomentumPressureKernels() override
void addMomentumTimeKernels() override
Functions adding kernels for the incompressible momentum equation If the material properties are not ...
virtual void addMomentumBoussinesqKernels() override
virtual void addMomentumFrictionKernels() override
virtual void addMomentumGravityKernels() override
void addPressureCorrectionKernels()
Function adding kernels for the incompressible pressure correction equation.

◆ addInitialConditions()

void WCNSFVFlowPhysicsBase::addInitialConditions ( )
overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 399 of file WCNSFVFlowPhysicsBase.C.

400{
401 if (!_define_variables && parameters().isParamSetByUser("initial_velocity") &&
402 parameters().isParamSetByUser("velocity_variable") &&
403 getParam<std::vector<FunctionName>>("initial_velocity").size() != 0)
404 // TODO: Rework and remove this last statement once the NSFV action is removed
405 paramError("initial_velocity",
406 "Velocity is defined externally of WCNSFVFlowPhysicsBase, so should the inital "
407 "conditions");
408 if (!_define_variables && parameters().isParamSetByUser("initial_pressure") &&
409 parameters().isParamSetByUser("pressure_variable"))
410 paramError("initial_pressure",
411 "Pressure is defined externally of WCNSFVFlowPhysicsBase, so should the inital "
412 "condition");
413
414 // Check dimension
415 if (getParam<std::vector<FunctionName>>("initial_velocity").size() != dimension() &&
416 getParam<std::vector<FunctionName>>("initial_velocity").size() != 3 &&
417 getParam<std::vector<FunctionName>>("initial_velocity").size() != 0)
418 // TODO: Rework and remove this last statement once the NSFV action is removed
419 paramError("initial_velocity",
420 "The number of velocity components in the " + type() + " initial condition is not " +
421 std::to_string(dimension()) + " or 3!");
422
423 InputParameters params = getFactory().getValidParams("FVFunctionIC");
424 assignBlocks(params, _blocks);
425 auto vvalue = getParam<std::vector<FunctionName>>("initial_velocity");
426
427 for (const auto d : make_range(dimension()))
428 {
429 params.set<VariableName>("variable") = _velocity_names[d];
430 params.set<FunctionName>("function") = vvalue[d];
431
433 _blocks,
434 /*whether IC is a default*/ !isParamSetByUser("initial_velocity"),
435 /*error if already an IC*/ isParamSetByUser("initial_velocity")))
437 "FVFunctionIC", prefix() + _velocity_names[d] + "_ic", params);
438 }
439
441 _blocks,
442 /*whether IC is a default*/ !isParamSetByUser("initial_pressure"),
443 /*error if already an IC*/ isParamSetByUser("initial_pressure")))
444 {
445 params.set<VariableName>("variable") = _pressure_name;
446 params.set<FunctionName>("function") = getParam<FunctionName>("initial_pressure");
447
448 getProblem().addFVInitialCondition("FVFunctionIC", prefix() + _pressure_name + "_ic", params);
449 }
450}
virtual void addFVInitialCondition(const std::string &ic_name, const std::string &name, InputParameters &parameters)
const std::string & type() const
const T & getParam(const std::string &name) const
bool _define_variables
Whether to define variables if they do not exist.
bool shouldCreateIC(const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool ic_is_default_ic, const bool error_if_already_defined) const
const std::vector< std::string > _velocity_names
Velocity names.

◆ addInletBC()

void WCNSLinearFVFlowPhysics::addInletBC ( )
overrideprivatevirtual

Functions adding boundary conditions for the flow simulation.

These are used for weakly-compressible simulations as well.

Implements WCNSFVFlowPhysicsBase.

Definition at line 421 of file WCNSLinearFVFlowPhysics.C.

422{
423 // Check the size of the BC parameters
424 unsigned int num_velocity_functor_inlets = 0;
425 for (const auto & [bdy, momentum_inlet_type] : _momentum_inlet_types)
426 if (momentum_inlet_type == "fixed-velocity" || momentum_inlet_type == "fixed-pressure")
427 num_velocity_functor_inlets++;
428
429 if (num_velocity_functor_inlets != _momentum_inlet_functors.size())
430 paramError("momentum_inlet_functors",
431 "Size (" + std::to_string(_momentum_inlet_functors.size()) +
432 ") is not the same as the number of entries in the momentum_inlet_types "
433 "subvector for fixed-velocities/pressures functors (size " +
434 std::to_string(num_velocity_functor_inlets) + ")");
435
436 unsigned int velocity_pressure_counter = 0;
437 for (const auto & [inlet_bdy, momentum_inlet_type] : _momentum_inlet_types)
438 {
439 if (momentum_inlet_type == "fixed-velocity")
440 {
441 const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
442 InputParameters params = getFactory().getValidParams(bc_type);
443 params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
444 if (_momentum_inlet_functors.size() < velocity_pressure_counter + 1)
445 paramError("momentum_inlet_functors",
446 "More non-flux inlets than inlet functors (" +
447 std::to_string(_momentum_inlet_functors.size()) + ")");
448
449 // Check that enough functors have been provided for the dimension of the problem
450 const auto momentum_functors = libmesh_map_find(_momentum_inlet_functors, inlet_bdy);
451 if (momentum_functors.size() < dimension())
452 paramError("momentum_inlet_functors",
453 "Subvector for boundary '" + inlet_bdy + "' (size " +
454 std::to_string(momentum_functors.size()) +
455 ") is not the same size as the number of dimensions of the physics (" +
456 std::to_string(dimension()) + ")");
457
458 for (const auto d : make_range(dimension()))
459 {
460 params.set<LinearVariableName>("variable") = _velocity_names[d];
461 params.set<MooseFunctorName>("functor") = momentum_functors[d];
462
463 getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + inlet_bdy, params);
464 }
465 ++velocity_pressure_counter;
466
467 // Add the two term BC expansion for pressure if requested
468 if (getParam<bool>("pressure_two_term_bc_expansion"))
469 {
470 const std::string bc_type = "LinearFVExtrapolatedPressureBC";
471 InputParameters params = getFactory().getValidParams(bc_type);
472 params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
473 params.set<LinearVariableName>("variable") = _pressure_name;
474 params.set<bool>("use_two_term_expansion") = true;
475 getProblem().addLinearFVBC(bc_type,
476 _pressure_name + "_extrapolation_inlet_" +
477 Moose::stringify(inlet_bdy),
478 params);
479 }
480 }
481 else if (momentum_inlet_type == "fixed-pressure")
482 {
483 const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
484 InputParameters params = getFactory().getValidParams(bc_type);
485 params.set<LinearVariableName>("variable") = _pressure_name;
486 if (_momentum_inlet_functors.size() < velocity_pressure_counter + 1)
487 paramError("momentum_inlet_functors",
488 "More non-flux inlets than inlet functors (" +
489 std::to_string(_momentum_inlet_functors.size()) + ")");
490
491 params.set<MooseFunctorName>("functor") =
492 libmesh_map_find(_momentum_inlet_functors, inlet_bdy)[0];
493 params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
494
495 getProblem().addLinearFVBC(bc_type, _pressure_name + "_" + inlet_bdy, params);
496 ++velocity_pressure_counter;
497 }
498 else
499 mooseError("Unsupported inlet boundary condition type: ", momentum_inlet_type);
500 }
501}
virtual void addLinearFVBC(const std::string &fv_bc_name, const std::string &name, InputParameters &parameters)
std::map< BoundaryName, MooseEnum > _momentum_inlet_types
Momentum inlet boundary types.
std::map< BoundaryName, std::vector< MooseFunctorName > > _momentum_inlet_functors
Functors describing the momentum inlet for each boundary.
std::string stringify(const T &t)

◆ addMaterials()

void WCNSFVFlowPhysicsBase::addMaterials ( )
overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 306 of file WCNSFVFlowPhysicsBase.C.

307{
310 else
312
315}
void addNonPorousMediumSpeedMaterial()
Add material to define the local speed with no porous medium treatment.
virtual bool hasForchheimerFriction() const =0
Return whether a Forchheimer friction model is in use.
void addFluidPropertiesFunctorMaterial()
Function which adds the general functor fluid properties functor material to define fluid functor mat...
void addPorousMediumSpeedMaterial()
Add material to define the local speed in porous medium flows.

◆ addMomentumBoussinesqKernels()

void WCNSLinearFVFlowPhysics::addMomentumBoussinesqKernels ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 394 of file WCNSLinearFVFlowPhysics.C.

395{
396 if (_compressibility == "weakly-compressible")
397 paramError("boussinesq_approximation",
398 "We cannot use boussinesq approximation while running in weakly-compressible mode!");
399
400 std::string kernel_type = "LinearFVMomentumBoussinesq";
401 std::string kernel_name = prefix() + "ins_momentum_boussinesq_";
402
403 InputParameters params = getFactory().getValidParams(kernel_type);
404 assignBlocks(params, _blocks);
405 params.set<VariableName>(NS::T_fluid) = _fluid_temperature_name;
406 params.set<MooseFunctorName>(NS::density) = _density_gravity_name;
407 params.set<RealVectorValue>("gravity") = getParam<RealVectorValue>("gravity");
408 params.set<Real>("ref_temperature") = getParam<Real>("ref_temperature");
409 params.set<MooseFunctorName>("alpha_name") = getParam<MooseFunctorName>("thermal_expansion");
410
411 for (const auto d : make_range(dimension()))
412 {
413 params.set<MooseEnum>("momentum_component") = NS::directions[d];
414 params.set<LinearVariableName>("variable") = _velocity_names[d];
415
416 getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
417 }
418}
virtual void addLinearFVKernel(const std::string &kernel_name, const std::string &name, InputParameters &parameters)
const MooseEnum _compressibility
Compressibility type, can be compressible, incompressible or weakly-compressible.
static const std::string directions[3]
Definition NS.h:23
VectorValue< Real > RealVectorValue

Referenced by addFVKernels().

◆ addMomentumFluxKernels()

void WCNSLinearFVFlowPhysics::addMomentumFluxKernels ( )
private

Definition at line 264 of file WCNSLinearFVFlowPhysics.C.

265{
266 const auto momentum_advection_method_name =
267 isParamValid("momentum_advection_interpolation_method_name")
268 ? getParam<InterpolationMethodName>("momentum_advection_interpolation_method_name")
269 : InterpolationMethodName(
270 std::string(getParam<MooseEnum>("momentum_advection_interpolation")));
271
272 const std::string u_names[3] = {"u", "v", "w"};
273 std::string kernel_type = "LinearWCNSFVMomentumFlux";
274 std::string kernel_name = prefix() + "ins_momentum_flux_";
275
276 InputParameters params = getFactory().getValidParams(kernel_type);
277 assignBlocks(params, _blocks);
279 params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
280 else
281 params.set<MooseFunctorName>(NS::mu) = NS::mu_eff;
282
283 params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
284 params.set<InterpolationMethodName>("advected_interp_method_name") =
285 momentum_advection_method_name;
286 params.set<bool>("use_nonorthogonal_correction") = _non_orthogonal_correction;
287 params.set<bool>("use_deviatoric_terms") = includeSymmetrizedViscousStress();
288
289 for (unsigned int i = 0; i < dimension(); ++i)
290 params.set<SolverVariableName>(u_names[i]) = _velocity_names[i];
291
292 for (const auto d : make_range(dimension()))
293 {
294 params.set<LinearVariableName>("variable") = _velocity_names[d];
295 params.set<MooseEnum>("momentum_component") = NS::directions[d];
296
297 getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
298 }
299}
const UserObjectName & rhieChowUOName() const
Return the name of the Rhie Chow user object.
const MooseFunctorName _dynamic_viscosity_name
Name of the dynamic viscosity material property.
bool includeSymmetrizedViscousStress() const
Whether to include the symmetrized contribution in the viscous stress.
static const std::string mu
Definition NS.h:127
static const std::string mu_eff
Definition NS.h:133

Referenced by addFVKernels().

◆ addMomentumFrictionKernels()

void WCNSLinearFVFlowPhysics::addMomentumFrictionKernels ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 320 of file WCNSLinearFVFlowPhysics.C.

321{
322 unsigned int num_friction_blocks = _friction_blocks.size();
323 unsigned int num_used_blocks = num_friction_blocks ? num_friction_blocks : 1;
324
325 const std::string kernel_type = "LinearFVMomentumFriction";
326 InputParameters params = getFactory().getValidParams(kernel_type);
327
328 for (const auto block_i : make_range(num_used_blocks))
329 {
330 std::string block_name = "";
331 if (num_friction_blocks)
332 {
333 params.set<std::vector<SubdomainName>>("block") = _friction_blocks[block_i];
334 block_name = Moose::stringify(_friction_blocks[block_i]);
335 }
336 else
337 {
338 assignBlocks(params, _blocks);
339 block_name = std::to_string(block_i);
340 }
341
342 for (const auto d : make_range(dimension()))
343 {
344 params.set<LinearVariableName>("variable") = _velocity_names[d];
345 params.set<MooseEnum>("momentum_component") = NS::directions[d];
346 for (unsigned int type_i = 0; type_i < _friction_types[block_i].size(); ++type_i)
347 {
348 const auto upper_name = MooseUtils::toUpper(_friction_types[block_i][type_i]);
349 if (upper_name == "DARCY")
350 {
351 params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
352 params.set<MooseFunctorName>("Darcy_name") = _friction_coeffs[block_i][type_i];
353 }
354 else
355 paramError("friction_types",
356 "Friction type '",
357 _friction_types[block_i][type_i],
358 "' is not implemented");
359 }
360
361 getProblem().addLinearFVKernel(kernel_type,
362 prefix() + "momentum_friction_" + block_name + "_" +
364 params);
365 }
366 }
367}
std::vector< std::vector< SubdomainName > > _friction_blocks
Subdomains where we want to have volumetric friction.
std::vector< std::vector< std::string > > _friction_coeffs
The coefficients used for each item if friction type.
std::string toUpper(std::string name)

Referenced by addFVKernels().

◆ addMomentumGravityKernels()

void WCNSLinearFVFlowPhysics::addMomentumGravityKernels ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 370 of file WCNSLinearFVFlowPhysics.C.

371{
373 {
374 std::string kernel_type = "LinearFVSource";
375 std::string kernel_name = prefix() + "ins_momentum_gravity_";
376
377 InputParameters params = getFactory().getValidParams(kernel_type);
378 assignBlocks(params, _blocks);
379 const auto gravity_vector = getParam<RealVectorValue>("gravity");
380 const std::vector<std::string> comp_axis({"x", "y", "z"});
381
382 for (const auto d : make_range(dimension()))
383 if (gravity_vector(d) != 0)
384 {
385 params.set<MooseFunctorName>("source_density") = "rho_g_" + comp_axis[d];
386 params.set<LinearVariableName>("variable") = _velocity_names[d];
387
388 getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
389 }
390 }
391}

Referenced by addFVKernels().

◆ addMomentumPressureKernels()

void WCNSLinearFVFlowPhysics::addMomentumPressureKernels ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 302 of file WCNSLinearFVFlowPhysics.C.

303{
304 std::string kernel_type = "LinearFVMomentumPressure";
305 std::string kernel_name = prefix() + "ins_momentum_pressure_";
306
307 InputParameters params = getFactory().getValidParams(kernel_type);
308 assignBlocks(params, _blocks);
309 params.set<VariableName>("pressure") = _pressure_name;
310
311 for (const auto d : make_range(dimension()))
312 {
313 params.set<MooseEnum>("momentum_component") = NS::directions[d];
314 params.set<LinearVariableName>("variable") = _velocity_names[d];
315 getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
316 }
317}

Referenced by addFVKernels().

◆ addMomentumTimeKernels()

void WCNSLinearFVFlowPhysics::addMomentumTimeKernels ( )
overrideprivatevirtual

Functions adding kernels for the incompressible momentum equation If the material properties are not constant, these can be used for weakly-compressible simulations (except the Boussinesq kernel) as well.

Implements WCNSFVFlowPhysicsBase.

Definition at line 246 of file WCNSLinearFVFlowPhysics.C.

247{
248 std::string kernel_type = "LinearFVTimeDerivative";
249 std::string kernel_name = prefix() + "ins_momentum_time";
250
251 InputParameters params = getFactory().getValidParams(kernel_type);
252 assignBlocks(params, _blocks);
253 params.set<MooseFunctorName>("factor") = _density_name;
254
255 for (const auto d : make_range(dimension()))
256 {
257 params.set<LinearVariableName>("variable") = _velocity_names[d];
258 if (shouldCreateTimeDerivative(_velocity_names[d], _blocks, /*error if already defined*/ false))
259 getProblem().addLinearFVKernel(kernel_type, kernel_name + "_" + NS::directions[d], params);
260 }
261}
bool shouldCreateTimeDerivative(const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool error_if_already_defined) const

Referenced by addFVKernels().

◆ addNonPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addNonPorousMediumSpeedMaterial ( )
protectedinherited

Add material to define the local speed with no porous medium treatment.

Definition at line 336 of file WCNSFVFlowPhysicsBase.C.

337{
338 // Not very future-proof but it works
339 const bool use_ad = !dynamic_cast<WCNSLinearFVFlowPhysics *>(this);
340 const std::string class_name =
341 use_ad ? "ADVectorMagnitudeFunctorMaterial" : "VectorMagnitudeFunctorMaterial";
342 InputParameters params = getFactory().getValidParams(class_name);
343 assignBlocks(params, _blocks);
344
345 const std::vector<std::string> param_names{"x_functor", "y_functor", "z_functor"};
346 for (unsigned int dim_i = 0; dim_i < dimension(); ++dim_i)
347 params.set<MooseFunctorName>(param_names[dim_i]) = _velocity_names[dim_i];
348 params.set<MooseFunctorName>("vector_magnitude_name") = NS::speed;
349
350 getProblem().addFunctorMaterial(class_name, prefix() + "ins_speed_material", params);
351}

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

◆ addOutletBC()

void WCNSLinearFVFlowPhysics::addOutletBC ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 504 of file WCNSLinearFVFlowPhysics.C.

505{
506 // Check the BCs size
507 unsigned int num_pressure_outlets = 0;
508 for (const auto & [bdy, momentum_outlet_type] : _momentum_outlet_types)
509 if (momentum_outlet_type == "fixed-pressure" ||
510 momentum_outlet_type == "fixed-pressure-zero-gradient")
511 num_pressure_outlets++;
512
513 if (num_pressure_outlets != _pressure_functors.size())
514 paramError("pressure_functors",
515 "Size (" + std::to_string(_pressure_functors.size()) +
516 ") is not the same as the number of pressure outlet boundaries in "
517 "'fixed-pressure/fixed-pressure-zero-gradient' (size " +
518 std::to_string(num_pressure_outlets) + ")");
519
520 const std::string u_names[3] = {"u", "v", "w"};
521 for (const auto & [outlet_bdy, momentum_outlet_type] : _momentum_outlet_types)
522 {
523 // Zero tangeantial gradient condition on velocity
524 if (momentum_outlet_type == "zero-gradient" || momentum_outlet_type == "fixed-pressure" ||
525 momentum_outlet_type == "fixed-pressure-zero-gradient")
526 {
527 const std::string bc_type = "LinearFVAdvectionDiffusionOutflowBC";
528 InputParameters params = getFactory().getValidParams(bc_type);
529 params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
530 params.set<bool>("use_two_term_expansion") = getParam<bool>("momentum_two_term_bc_expansion");
531
532 for (const auto d : make_range(dimension()))
533 {
534 params.set<LinearVariableName>("variable") = _velocity_names[d];
535 getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + outlet_bdy, params);
536 }
537 }
538
539 // Fixed pressure condition, coming in the pressure correction equation
540 if (momentum_outlet_type == "fixed-pressure" ||
541 momentum_outlet_type == "fixed-pressure-zero-gradient")
542 {
543 const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
544 InputParameters params = getFactory().getValidParams(bc_type);
545 params.set<LinearVariableName>("variable") = _pressure_name;
546 params.set<MooseFunctorName>("functor") = libmesh_map_find(_pressure_functors, outlet_bdy);
547 params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
548
549 getProblem().addLinearFVBC(bc_type, _pressure_name + "_" + outlet_bdy, params);
550 }
551 }
552}
std::map< BoundaryName, MooseEnum > _momentum_outlet_types
Momentum outlet boundary types.
std::map< BoundaryName, MooseFunctorName > _pressure_functors
Functors describing the outlet pressure on each boundary.

◆ addPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addPorousMediumSpeedMaterial ( )
protectedinherited

Add material to define the local speed in porous medium flows.

Definition at line 318 of file WCNSFVFlowPhysicsBase.C.

319{
320 InputParameters params = getFactory().getValidParams("PINSFVSpeedFunctorMaterial");
321 assignBlocks(params, _blocks);
322
323 for (unsigned int dim_i = 0; dim_i < dimension(); ++dim_i)
324 params.set<MooseFunctorName>(NS::superficial_velocity_vector[dim_i]) = _velocity_names[dim_i];
326 params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
327 else
328 params.set<MooseFunctorName>(NS::porosity) = "1";
329 params.set<bool>("define_interstitial_velocity_components") = _porous_medium_treatment;
330
332 "PINSFVSpeedFunctorMaterial", prefix() + "pins_speed_material", params);
333}
MooseFunctorName _flow_porosity_functor_name
Name of the porosity functor for the flow equations (if smoothed)
static const std::string superficial_velocity_vector[3]
Definition NS.h:55

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

◆ addPostprocessors()

void WCNSFVFlowPhysicsBase::addPostprocessors ( )
overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 460 of file WCNSFVFlowPhysicsBase.C.

461{
462 const auto momentum_inlet_types = getParam<MultiMooseEnum>("momentum_inlet_types");
463
464 for (unsigned int bc_ind = 0; bc_ind < momentum_inlet_types.size(); ++bc_ind)
465 if (momentum_inlet_types[bc_ind] == "flux-mass" ||
466 momentum_inlet_types[bc_ind] == "flux-velocity")
467 {
468 const std::string pp_type = "AreaPostprocessor";
469 InputParameters params = getFactory().getValidParams(pp_type);
470 params.set<std::vector<BoundaryName>>("boundary") = {_inlet_boundaries[bc_ind]};
471 params.set<ExecFlagEnum>("execute_on") = EXEC_INITIAL;
472
473 const auto name_pp = "area_pp_" + _inlet_boundaries[bc_ind];
474 if (!getProblem().hasUserObject(name_pp))
475 getProblem().addPostprocessor(pp_type, name_pp, params);
476 }
477}
virtual void addPostprocessor(const std::string &pp_name, const std::string &name, InputParameters &parameters)
bool hasUserObject(const std::string &name) const
const std::vector< BoundaryName > _inlet_boundaries
Boundaries with a flow inlet specified on them.

◆ addPressureCorrectionKernels()

void WCNSLinearFVFlowPhysics::addPressureCorrectionKernels ( )
private

Function adding kernels for the incompressible pressure correction equation.

Definition at line 217 of file WCNSLinearFVFlowPhysics.C.

218{
219 {
220 std::string kernel_type = "LinearFVPressureCorrectionDiffusion";
221 std::string kernel_name = prefix() + "p_diffusion";
222
223 InputParameters params = getFactory().getValidParams(kernel_type);
224 assignBlocks(params, _blocks);
225 params.set<LinearVariableName>("variable") = _pressure_name;
226 params.set<MooseFunctorName>("diffusion_tensor") = "Ainv";
227 params.set<bool>("use_nonorthogonal_correction") = _non_orthogonal_correction;
228
229 getProblem().addLinearFVKernel(kernel_type, kernel_name, params);
230 }
231 {
232 std::string kernel_type = "LinearFVDivergence";
233 std::string kernel_name = prefix() + "HbyA_divergence";
234
235 InputParameters params = getFactory().getValidParams(kernel_type);
236 assignBlocks(params, _blocks);
237 params.set<LinearVariableName>("variable") = _pressure_name;
238 params.set<MooseFunctorName>("face_flux") = "HbyA";
239 params.set<bool>("force_boundary_execution") = true;
240
241 getProblem().addLinearFVKernel(kernel_type, kernel_name, params);
242 }
243}

Referenced by addFVKernels().

◆ addRhieChowUserObjects()

void WCNSLinearFVFlowPhysics::addRhieChowUserObjects ( )
overrideprivatevirtual

Function which adds the RhieChow interpolator user objects for weakly and incompressible formulations.

Implements WCNSFVFlowPhysicsBase.

Definition at line 646 of file WCNSLinearFVFlowPhysics.C.

647{
648 mooseAssert(dimension(), "0-dimension not supported");
649
650 // First make sure that we only add this object once
651 // Potential cases:
652 // - there is a flow physics, and an advection one (UO should be added by one)
653 // - there is only an advection physics (UO should be created)
654 // - there are two advection physics on different blocks with set velocities (first one picks)
655 // Counting RC UOs defined on the same blocks seems to be the most fool proof option
656 std::vector<UserObject *> objs;
657 getProblem()
658 .theWarehouse()
659 .query()
660 .condition<AttribSystem>("UserObject")
661 .condition<AttribThread>(0)
662 .queryInto(objs);
663 unsigned int num_rc_uo = 0;
664 for (const auto & obj : objs)
665 if (dynamic_cast<RhieChowMassFlux *>(obj))
666 {
667 const auto rc_obj = dynamic_cast<RhieChowMassFlux *>(obj);
668 if (rc_obj->blocks() == _blocks)
669 num_rc_uo++;
670 // one of the RC user object is defined everywhere
671 else if (rc_obj->blocks().size() == 0 || _blocks.size() == 0)
672 num_rc_uo++;
673 }
674
675 if (num_rc_uo)
676 return;
677
678 const std::string u_names[3] = {"u", "v", "w"};
679 const auto object_type = "RhieChowMassFlux";
680
681 auto params = getFactory().getValidParams(object_type);
682 assignBlocks(params, _blocks);
683 for (unsigned int d = 0; d < dimension(); ++d)
684 params.set<VariableName>(u_names[d]) = _velocity_names[d];
685
686 params.set<VariableName>("pressure") = _pressure_name;
687 params.set<std::string>("p_diffusion_kernel") = prefix() + "p_diffusion";
688 params.set<MooseFunctorName>(NS::density) = _density_name;
689 params.set<MooseEnum>("pressure_projection_method") =
690 getParam<MooseEnum>("pressure_projection_method");
691 params.set<MooseEnum>("pressure_diffusion_interpolation") =
692 getParam<MooseEnum>("pressure_diffusion_interpolation");
693
694 addUserObject(object_type, rhieChowUOName(), params);
695}
TheWarehouse & theWarehouse() const
void addUserObject(const std::string &uo_type, const std::string &uo_name, InputParameters &params)
User object responsible for determining the face fluxes using the Rhie-Chow interpolation in a segreg...
Query query()

Referenced by addUserObjects().

◆ addSeparatorBC()

virtual void WCNSLinearFVFlowPhysics::addSeparatorBC ( )
inlineoverrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 60 of file WCNSLinearFVFlowPhysics.h.

60{}

◆ addSolverVariables()

void WCNSLinearFVFlowPhysics::addSolverVariables ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 130 of file WCNSLinearFVFlowPhysics.C.

131{
133 return;
134
135 for (const auto d : make_range(dimension()))
138
139 const std::vector<std::string> v_short = {"u", "v", "w"};
140
141 // Check number of variables
142 if (_velocity_names.size() != dimension() && _velocity_names.size() != 3)
143 paramError("velocity_variable",
144 "The number of velocity variable names supplied to the NSFVAction is not " +
145 Moose::stringify(dimension()) + " (mesh dimension)" +
146 ((dimension() == 3) ? "" : " or 3!") + "\nVelocity variables " +
148
149 // Velocities
150 for (const auto d : make_range(dimension()))
151 {
152 if (!shouldCreateVariable(_velocity_names[d], _blocks, /*error if aux*/ true))
153 reportPotentiallyMissedParameters({"system_names"}, "MooseLinearVariableFVReal");
154 else if (_define_variables)
155 {
156 std::string variable_type = "MooseLinearVariableFVReal";
157
158 auto params = getFactory().getValidParams(variable_type);
159
160 addSharedSolverVariable(variable_type, _velocity_names[d], params);
161 }
162 else
163 paramError("velocity_variable",
164 "Variable (" + _velocity_names[d] +
165 ") supplied to the WCNSLinearFVFlowPhysics does not exist!");
166 }
167
168 // Pressure
169 if (!shouldCreateVariable(_pressure_name, _blocks, /*error if aux*/ true))
170 reportPotentiallyMissedParameters({"system_names"}, "MooseLinearVariableFVReal");
171 else if (_define_variables)
172 {
173 const auto pressure_type = "MooseLinearVariableFVReal";
174
175 auto params = getFactory().getValidParams(pressure_type);
176
177 addSharedSolverVariable(pressure_type, _pressure_name, params);
178 }
179 else
180 paramError("pressure_variable",
181 "Variable (" + _pressure_name +
182 ") supplied to the WCNSLinearFVFlowPhysics does not exist!");
183}
void reportPotentiallyMissedParameters(const std::vector< std::string > &param_names, const std::string &object_type, const std::string &object_name="") const
void saveSolverVariableName(const VariableName &var_name)
void addSharedSolverVariable(const std::string &var_type, const VariableName &var_name, InputParameters params)
bool shouldCreateVariable(const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool error_if_aux)

◆ addUserObjects()

void WCNSLinearFVFlowPhysics::addUserObjects ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 638 of file WCNSLinearFVFlowPhysics.C.

639{
640 mooseAssert(!_porous_medium_treatment, "Not implemented");
641 // Rhie Chow user object for interpolation velocities
643}
virtual void addRhieChowUserObjects() override
Function which adds the RhieChow interpolator user objects for weakly and incompressible formulations...

◆ addWallsBC()

void WCNSLinearFVFlowPhysics::addWallsBC ( )
overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 555 of file WCNSLinearFVFlowPhysics.C.

556{
557 const std::string u_names[3] = {"u", "v", "w"};
558 bool has_symmetry_bc = false;
559
560 for (const auto & [boundary_name, momentum_wall_type] : _momentum_wall_types)
561 {
562 if (momentum_wall_type == "noslip")
563 {
564 const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
565 InputParameters params = getFactory().getValidParams(bc_type);
566 params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
567
568 for (const auto d : make_range(dimension()))
569 {
570 params.set<LinearVariableName>("variable") = _velocity_names[d];
571 if (_momentum_wall_functors.count(boundary_name) == 0)
572 params.set<MooseFunctorName>("functor") = "0";
573 else
574 params.set<MooseFunctorName>("functor") = _momentum_wall_functors[boundary_name][d];
575
576 getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
577 }
578 }
579 else if (momentum_wall_type == "symmetry")
580 {
581 has_symmetry_bc = true;
582 {
583 const std::string bc_type = "LinearFVVelocitySymmetryBC";
584 InputParameters params = getFactory().getValidParams(bc_type);
585 params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
586 for (unsigned int d = 0; d < dimension(); ++d)
587 params.set<SolverVariableName>(u_names[d]) = _velocity_names[d];
588
589 for (const auto d : make_range(dimension()))
590 {
591 params.set<LinearVariableName>("variable") = _velocity_names[d];
592 params.set<MooseEnum>("momentum_component") = NS::directions[d];
593
594 getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
595 }
596 }
597 {
598 const std::string bc_type = "LinearFVPressureSymmetryBC";
599 InputParameters params = getFactory().getValidParams(bc_type);
600 params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
601 params.set<LinearVariableName>("variable") = _pressure_name;
602 params.set<MooseFunctorName>("HbyA_flux") = "HbyA";
603 getProblem().addLinearFVBC(bc_type, _pressure_name + "_" + boundary_name, params);
604 }
605 }
606 else
607 mooseError("Unsupported wall boundary condition type: " + std::string(momentum_wall_type));
608 }
609
610 if (getParam<bool>("pressure_two_term_bc_expansion"))
611 {
612 if (!has_symmetry_bc)
613 {
614 const std::string bc_type = "LinearFVExtrapolatedPressureBC";
615 InputParameters params = getFactory().getValidParams(bc_type);
616 params.set<std::vector<BoundaryName>>("boundary") = _wall_boundaries;
617 params.set<LinearVariableName>("variable") = _pressure_name;
618 params.set<bool>("use_two_term_expansion") = true;
620 bc_type, _pressure_name + "_extrapolation_" + Moose::stringify(_wall_boundaries), params);
621 }
622 else
623 for (const auto & [boundary_name, momentum_wall_type] : _momentum_wall_types)
624 if (momentum_wall_type != "symmetry")
625 {
626 const std::string bc_type = "LinearFVExtrapolatedPressureBC";
627 InputParameters params = getFactory().getValidParams(bc_type);
628 params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
629 params.set<LinearVariableName>("variable") = _pressure_name;
630 params.set<bool>("use_two_term_expansion") = true;
632 bc_type, _pressure_name + "_extrapolation_" + boundary_name, params);
633 }
634 }
635}
std::map< BoundaryName, std::vector< MooseFunctorName > > _momentum_wall_functors
Functors describing the momentum for each wall boundary.
const std::vector< BoundaryName > _wall_boundaries
Boundaries which define a wall (slip/noslip/etc.)
std::map< BoundaryName, MooseEnum > _momentum_wall_types
Momentum wall boundary types.

◆ compressibility()

const MooseEnum & WCNSFVFlowPhysicsBase::compressibility ( ) const
inlineinherited

Return the compressibility of the flow equations selected.

Definition at line 53 of file WCNSFVFlowPhysicsBase.h.

53{ return _compressibility; }

◆ densityName()

const MooseFunctorName & WCNSFVFlowPhysicsBase::densityName ( ) const
inlineinherited

◆ dynamicViscosityName()

const MooseFunctorName & WCNSFVFlowPhysicsBase::dynamicViscosityName ( ) const
inlineinherited

◆ getAdditionalRMParams()

InputParameters NavierStokesPhysicsBase::getAdditionalRMParams ( ) const
overrideprotectedvirtualinherited

Parameters to change or add relationship managers.

Reimplemented from PhysicsBase.

Definition at line 56 of file NavierStokesPhysicsBase.C.

57{
58 unsigned short necessary_layers = getParam<unsigned short>("ghost_layers");
59 necessary_layers = std::max(necessary_layers, getNumberAlgebraicGhostingLayersNeeded());
60
61 // Just an object that has a ghost_layers parameter
62 const std::string kernel_type = "INSFVMixingLengthReynoldsStress";
63 InputParameters params = getFactory().getValidParams(kernel_type);
64 params.template set<unsigned short>("ghost_layers") = necessary_layers;
65
66 return params;
67}
virtual unsigned short getNumberAlgebraicGhostingLayersNeeded() const =0
Return the number of ghosting layers needed.

◆ getAxisymmetricRZBlocks()

std::vector< SubdomainName > WCNSFVFlowPhysicsBase::getAxisymmetricRZBlocks ( ) const
protectedinherited

Return the set of blocks restricted to an RZ coordinate system.

Definition at line 539 of file WCNSFVFlowPhysicsBase.C.

540{
541 std::vector<SubdomainName> rz_blocks;
542 const auto & mesh = getProblem().mesh();
543
544 const bool use_all_blocks =
545 _blocks.empty() || allMeshBlocks(_blocks) ||
546 std::find(_blocks.begin(), _blocks.end(), "ANY_BLOCK_ID") != _blocks.end();
547
548 std::vector<SubdomainID> block_ids;
549 if (use_all_blocks)
550 {
551 const auto & mesh_blocks = mesh.meshSubdomains();
552 block_ids.insert(block_ids.end(), mesh_blocks.begin(), mesh_blocks.end());
553 }
554 else
555 block_ids = mesh.getSubdomainIDs(_blocks);
556
557 for (const auto subdomain_id : block_ids)
558 if (mesh.getCoordSystem(subdomain_id) == Moose::COORD_RZ)
559 {
560 auto name = mesh.getSubdomainName(subdomain_id);
561 if (name.empty())
562 name = Moose::stringify(subdomain_id);
563 rz_blocks.push_back(name);
564 }
565
566 return rz_blocks;
567}
const std::string & name() const
bool allMeshBlocks(const std::vector< SubdomainName > &blocks) const
MeshBase & mesh

Referenced by WCNSFVTurbulencePhysics::addAxisymmetricTurbulentViscousSource(), and WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource().

◆ getCoupledTurbulencePhysics()

const WCNSFVTurbulencePhysicsBase * WCNSFVFlowPhysicsBase::getCoupledTurbulencePhysics ( ) const
protectedinherited

Find the turbulence physics.

Definition at line 506 of file WCNSFVFlowPhysicsBase.C.

507{
508 // User passed it, just use that
509 if (isParamValid("coupled_turbulence_physics"))
510 return getCoupledPhysics<WCNSFVTurbulencePhysicsBase>(
511 getParam<PhysicsName>("coupled_turbulence_physics"));
512 // Look for any physics of the right type, and check the block restriction
513 else
514 {
515 const auto all_turbulence_physics = getCoupledPhysics<const WCNSFVTurbulencePhysicsBase>(true);
516 for (const auto physics : all_turbulence_physics)
517 {
519 physics->name(), physics->blocks(), /*error_if_not_identical=*/false))
520 return physics;
521 else if (_verbose)
522 mooseInfoRepeated("Detected Turbulence Physics '" + physics->name() +
523 "' with an incompatible block restriction. It will thus not be coupled "
524 "to this flow equations physics");
525 }
526 }
527 // Did not find one
528 return nullptr;
529}
void mooseInfoRepeated(Args &&... args)
bool checkBlockRestrictionIdentical(const std::string &object_name, const std::vector< SubdomainName > &blocks, const bool error_if_not_identical=true) const
const bool _verbose

Referenced by WCNSFVFlowPhysicsBase::actOnAdditionalTasks().

◆ getFlowVariableName()

VariableName WCNSFVFlowPhysicsBase::getFlowVariableName ( const std::string &  default_name) const
protectedinherited

Convenience routine to be able to retrieve the actual variable names from their default names.

Definition at line 480 of file WCNSFVFlowPhysicsBase.C.

481{
482 if (short_name == NS::pressure)
483 return getPressureName();
484 else if (short_name == NS::velocity_x && dimension() > 0)
485 return getVelocityNames()[0];
486 else if (short_name == NS::velocity_y && dimension() > 1)
487 return getVelocityNames()[1];
488 else if (short_name == NS::velocity_z && dimension() > 2)
489 return getVelocityNames()[2];
490 else if (short_name == NS::temperature)
492 else
493 mooseError("Short Variable name '", short_name, "' not recognized.");
494}
const NonlinearVariableName & getPressureName() const
const std::vector< std::string > & getVelocityNames() const
To interface with other Physics.
const NonlinearVariableName & getFluidTemperatureName() const
static const std::string velocity_y
Definition NS.h:48
static const std::string temperature
Definition NS.h:60
static const std::string velocity_z
Definition NS.h:49
static const std::string velocity_x
Definition NS.h:47

◆ getFluidTemperatureName()

const NonlinearVariableName & WCNSFVFlowPhysicsBase::getFluidTemperatureName ( ) const
inlineinherited

◆ getFluxInletDirections()

const std::vector< Point > & WCNSFVFlowPhysicsBase::getFluxInletDirections ( ) const
inlineinherited

Get the inlet direction if using a flux inlet.

Definition at line 89 of file WCNSFVFlowPhysicsBase.h.

89{ return _flux_inlet_directions; }
std::vector< Point > _flux_inlet_directions
Direction of each flux inlet. Indexing based on the number of flux boundaries.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC(), and WCNSFVScalarTransportPhysics::addScalarInletBC().

◆ getFluxInletPPs()

const std::vector< PostprocessorName > & WCNSFVFlowPhysicsBase::getFluxInletPPs ( ) const
inlineinherited

Get the inlet flux postprocessor if using a flux inlet.

Definition at line 91 of file WCNSFVFlowPhysicsBase.h.

91{ return _flux_inlet_pps; }
std::vector< PostprocessorName > _flux_inlet_pps
Postprocessors describing the momentum inlet for each boundary. Indexing based on the number of flux ...

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC(), and WCNSFVScalarTransportPhysics::addScalarInletBC().

◆ getHydraulicSeparators()

const std::vector< BoundaryName > & WCNSFVFlowPhysicsBase::getHydraulicSeparators ( ) const
inlineinherited

Get the hydraulic separator boundaries.

Definition at line 81 of file WCNSFVFlowPhysicsBase.h.

81{ return _hydraulic_separators; }

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergySeparatorBC().

◆ getInletBoundaries()

const std::vector< BoundaryName > & WCNSFVFlowPhysicsBase::getInletBoundaries ( ) const
inlineinherited

◆ getLinearFrictionCoefName()

virtual MooseFunctorName WCNSLinearFVFlowPhysics::getLinearFrictionCoefName ( ) const
inlineoverrideprivatevirtual

Get the name of the linear friction coefficient. Returns an empty string if no friction.

Implements WCNSFVFlowPhysicsBase.

Definition at line 67 of file WCNSLinearFVFlowPhysics.h.

68 {
69 mooseError("Not implemented");
70 }

◆ getMomentumAdvectionFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getMomentumAdvectionFaceInterpolationMethod ( ) const
inlineinherited

Get the face interpolation method for momentum in the advection term.

Definition at line 65 of file WCNSFVFlowPhysicsBase.h.

66 {
68 }
const MooseEnum _momentum_advection_interpolation
The momentum face interpolation method for being advected.

◆ getMomentumFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getMomentumFaceInterpolationMethod ( ) const
inlineinherited

Get the face interpolation method for momentum (mostly used in the stress terms)

Definition at line 70 of file WCNSFVFlowPhysicsBase.h.

71 {
73 }
const MooseEnum _momentum_face_interpolation
The momentum face interpolation method for stress terms.

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm(), and WCNSFVTurbulencePhysics::addFlowTurbulenceKernels().

◆ getNumberAlgebraicGhostingLayersNeeded()

unsigned short WCNSLinearFVFlowPhysics::getNumberAlgebraicGhostingLayersNeeded ( ) const
overrideprivatevirtual

Return the number of algebraic ghosting layers needed.

Reimplemented from WCNSFVFlowPhysicsBase.

Definition at line 730 of file WCNSLinearFVFlowPhysics.C.

731{
732 return 1;
733}

◆ getOutletBoundaries()

const std::vector< BoundaryName > & WCNSFVFlowPhysicsBase::getOutletBoundaries ( ) const
inlineinherited

Get the outlet boundaries.

Definition at line 77 of file WCNSFVFlowPhysicsBase.h.

77{ return _outlet_boundaries; }
const std::vector< BoundaryName > _outlet_boundaries
Boundaries with a flow outlet specified on them.

Referenced by WCNSLinearFVFluidHeatTransferPhysics::addEnergyOutletBC(), and WCNSLinearFVScalarTransportPhysics::addScalarOutletBC().

◆ getPorosityFunctorName()

MooseFunctorName WCNSFVFlowPhysicsBase::getPorosityFunctorName ( const bool  smoothed) const
inherited

◆ getPressureName()

const NonlinearVariableName & WCNSFVFlowPhysicsBase::getPressureName ( ) const
inlineinherited

◆ getSuppliedUserObjects()

std::vector< UserObjectName > WCNSLinearFVFlowPhysics::getSuppliedUserObjects ( ) const
overridevirtual

Reimplemented from PhysicsBase.

Definition at line 698 of file WCNSLinearFVFlowPhysics.C.

699{
700 return {rhieChowUOName()};
701}

◆ getVelocityFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getVelocityFaceInterpolationMethod ( ) const
inlineinherited

Get the face interpolation method for velocity.

Definition at line 63 of file WCNSFVFlowPhysicsBase.h.

const MooseEnum _velocity_interpolation
The velocity face interpolation method for advecting other quantities.

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm().

◆ getVelocityNames()

const std::vector< std::string > & WCNSFVFlowPhysicsBase::getVelocityNames ( ) const
inlineinherited

◆ getWallBoundaries()

const std::vector< BoundaryName > & WCNSFVFlowPhysicsBase::getWallBoundaries ( ) const
inlineinherited

◆ gravityVector()

RealVectorValue WCNSFVFlowPhysicsBase::gravityVector ( ) const
inlineinherited

Return the gravity vector.

Definition at line 57 of file WCNSFVFlowPhysicsBase.h.

57{ return getParam<RealVectorValue>("gravity"); }

Referenced by WCNSFVTwoPhaseMixturePhysics::addFunctorMaterials(), WCNSLinearFVTwoPhaseMixturePhysics::addMaterials(), and WCNSLinearFVTwoPhaseMixturePhysics::checkIntegrity().

◆ hasFlowEquations()

bool WCNSFVFlowPhysicsBase::hasFlowEquations ( ) const
inlineinherited

◆ hasForchheimerFriction()

virtual bool WCNSLinearFVFlowPhysics::hasForchheimerFriction ( ) const
inlineoverrideprivatevirtual

Return whether a Forchheimer friction model is in use.

Implements WCNSFVFlowPhysicsBase.

Definition at line 62 of file WCNSLinearFVFlowPhysics.h.

62{ return false; };

◆ hasTurbulencePhysics()

bool WCNSFVFlowPhysicsBase::hasTurbulencePhysics ( ) const
inlineprotectedinherited

Whether a turbulence Physics has been coupled in, to know which viscosity to pick on symmetry boundary conditions.

Definition at line 156 of file WCNSFVFlowPhysicsBase.h.

157 {
160 else
161 return false;
162 }
bool hasTurbulenceModel() const
Whether a turbulence model is in use.

Referenced by WCNSFVFlowPhysics::addWallsBC().

◆ includeIsotropicStress()

bool WCNSFVFlowPhysicsBase::includeIsotropicStress ( ) const
inlineinherited

Whether to include the isotropic viscous stress contribution.

Definition at line 43 of file WCNSFVFlowPhysicsBase.h.

const bool _include_isotropic_viscous_stress
Whether to include the isotropic viscous stress contribution.

Referenced by WCNSFVTurbulencePhysics::addFlowTurbulenceKernels(), and WCNSFVFlowPhysics::addMomentumViscousDissipationKernels().

◆ includeSymmetrizedViscousStress()

bool WCNSFVFlowPhysicsBase::includeSymmetrizedViscousStress ( ) const
inlineinherited

Whether to include the symmetrized contribution in the viscous stress.

Definition at line 41 of file WCNSFVFlowPhysicsBase.h.

const bool _include_symmetrized_viscous_stress
Whether to include the symmetrized viscous stress contribution.

Referenced by WCNSFVTurbulencePhysics::addAxisymmetricTurbulentViscousSource(), WCNSFVFlowPhysics::addAxisymmetricViscousSourceKernel(), addMomentumFluxKernels(), and WCNSFVFlowPhysics::addMomentumViscousDissipationKernels().

◆ initializePhysicsAdditional()

void WCNSLinearFVFlowPhysics::initializePhysicsAdditional ( )
overrideprotectedvirtual

Reimplemented from WCNSFVFlowPhysicsBase.

Definition at line 110 of file WCNSLinearFVFlowPhysics.C.

111{
113 // TODO Add support for multi-system by either:
114 // - creating the problem in the Physics or,
115 // - checking that the right systems are being created
116 getProblem().needSolutionState(2, Moose::SolutionIterationType::Nonlinear);
117 // TODO Ban all other nonlinear Physics for now
118
119 // Fix the default system names if using a different dimension
120 if (!isParamSetByUser("system_name"))
121 {
122 if (dimension() == 1)
123 _system_names = {"u_system", "pressure_system"};
124 else if (dimension() == 2)
125 _system_names = {"u_system", "v_system", "pressure_system"};
126 }
127}
void needSolutionState(unsigned int oldest_needed, Moose::SolutionIterationType iteration_type)
std::vector< SolverSystemName > _system_names
virtual void initializePhysicsAdditional() override

◆ inletBoundaryType()

NS::MomentumInletTypes WCNSFVFlowPhysicsBase::inletBoundaryType ( const BoundaryName &  boundary_name) const
inlineinherited

Get the type of the inlet BC.

Definition at line 83 of file WCNSFVFlowPhysicsBase.h.

84 {
86 static_cast<int>(libmesh_map_find(_momentum_inlet_types, boundary_name)));
87 }
MomentumInletTypes
Definition NS.h:217

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC().

◆ porousMediumTreatment()

bool WCNSFVFlowPhysicsBase::porousMediumTreatment ( ) const
inlineinherited

Return whether a porous medium treatment is applied.

Definition at line 55 of file WCNSFVFlowPhysicsBase.h.

◆ rhieChowUOName()

const UserObjectName & WCNSFVFlowPhysicsBase::rhieChowUOName ( ) const
inherited

Return the name of the Rhie Chow user object.

Definition at line 532 of file WCNSFVFlowPhysicsBase.C.

533{
534 mooseAssert(!_rc_uo_name.empty(), "The Rhie-Chow user-object name should be set!");
535 return _rc_uo_name;
536}
UserObjectName _rc_uo_name
name of the Rhie Chow user object

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm(), WCNSFVTurbulencePhysics::addAxisymmetricTurbulentViscousSource(), WCNSFVFlowPhysics::addAxisymmetricViscousSourceKernel(), WCNSFVFluidHeatTransferPhysics::addEnergyAdvectionKernels(), WCNSLinearFVFluidHeatTransferPhysics::addEnergyAdvectionKernels(), WCNSFVTurbulencePhysics::addFlowTurbulenceKernels(), WCNSFVFlowPhysics::addInletBC(), WCNSFVTurbulencePhysics::addKEpsilonAdvection(), WCNSLinearFVTurbulencePhysics::addKEpsilonAdvection(), WCNSFVFlowPhysics::addMassKernels(), WCNSFVFlowPhysics::addMomentumAdvectionKernels(), WCNSFVFlowPhysics::addMomentumBoussinesqKernels(), addMomentumFluxKernels(), WCNSFVFlowPhysics::addMomentumFrictionKernels(), WCNSFVFlowPhysics::addMomentumGravityKernels(), WCNSFVFlowPhysics::addMomentumPressureKernels(), WCNSFVFlowPhysics::addMomentumTimeKernels(), WCNSFVFlowPhysics::addMomentumViscousDissipationKernels(), WCNSFVFlowPhysics::addOutletBC(), WCNSFVTwoPhaseMixturePhysics::addPhaseDriftFluxTerm(), WCNSLinearFVTwoPhaseMixturePhysics::addPhaseDriftFluxTerm(), addRhieChowUserObjects(), WCNSFVScalarTransportPhysics::addScalarAdvectionKernels(), WCNSLinearFVScalarTransportPhysics::addScalarAdvectionKernels(), WCNSFVFlowPhysics::addSeparatorBC(), WCNSFVFlowPhysics::addWallsBC(), WCNSFVFlowPhysics::getSuppliedUserObjects(), and getSuppliedUserObjects().

◆ usingNavierStokesFVSyntax()

bool NavierStokesPhysicsBase::usingNavierStokesFVSyntax ( ) const
inlineprotectedinherited

Detects if we are using the new Physics syntax or the old NavierStokesFV action.

Definition at line 33 of file NavierStokesPhysicsBase.h.

34 {
35 return (parameters().get<std::string>("registered_identifier") == "Modules/NavierStokesFV");
36 }

◆ validParams()

InputParameters WCNSLinearFVFlowPhysics::validParams ( )
static

Definition at line 27 of file WCNSLinearFVFlowPhysics.C.

28{
31 "Define the Navier Stokes weakly-compressible equations with the linear "
32 "solver implementation of the SIMPLE scheme");
33 params.set<MooseEnum>("momentum_advection_interpolation") = NS::fvAdvectedInterpolationMethods();
34 params.addParam<InterpolationMethodName>(
35 "momentum_advection_interpolation_method_name",
36 "Name of an externally defined FVInterpolationMethod to use for momentum advection. When "
37 "provided, this overrides 'momentum_advection_interpolation'.");
38
39 params.addParam<bool>(
40 "orthogonality_correction", false, "Whether to use orthogonality correction");
41 params.renameParam("orthogonality_correction", "use_nonorthogonal_correction", "");
42 params.addParam<MooseEnum>(
43 "pressure_diffusion_interpolation",
45 "The face interpolation method for Ainv in the pressure correction diffusion term.");
46 params.addParamNamesToGroup("pressure_diffusion_interpolation", "Numerical scheme");
47 params.set<unsigned short>("ghost_layers") = 1;
48
49 // This will be adapted based on the dimension
50 params.set<std::vector<SolverSystemName>>("system_names") = {
51 "u_system", "v_system", "w_system", "pressure_system"};
52
53 // Implemented in the executioner
54 params.suppressParameter<MooseEnum>("pinned_pressure_type");
55 params.suppressParameter<Point>("pinned_pressure_point");
56 params.suppressParameter<PostprocessorName>("pinned_pressure_value");
57
58 // Not supported
59 params.suppressParameter<bool>("add_flow_equations");
60 params.set<bool>("porous_medium_treatment") = false;
61 params.suppressParameter<bool>("porous_medium_treatment");
62 params.set<MooseFunctorName>("porosity") = "1";
63 params.suppressParameter<MooseFunctorName>("porosity");
64 params.suppressParameter<MooseEnum>("mu_interp_method");
65 // Not needed
66 params.suppressParameter<bool>("add_flow_equations");
67 params.suppressParameter<MooseEnum>("preconditioning");
68
69 // No other options so far
70 params.set<MooseEnum>("velocity_interpolation") = "rc";
71 params.suppressParameter<MooseEnum>("velocity_interpolation");
72
73 // Rhie-Chow
74 params.transferParam<MooseEnum>(RhieChowMassFlux::validParams(), "pressure_projection_method");
75 params.addParamNamesToGroup("momentum_advection_interpolation_method_name", "Numerical scheme");
76
77 return params;
78}
void suppressParameter(const std::string &name)
void addParamNamesToGroup(const std::string &space_delim_names, const std::string group_name)
void renameParam(const std::string &old_name, const std::string &new_name, const std::string &new_docstring)
void transferParam(const InputParameters &source_param, const std::string &name, const std::string &new_name="", const std::string &new_description="")
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
void addClassDescription(const std::string &doc_string)
static InputParameters validParams()
static InputParameters validParams()
MooseEnum fvFaceInterpolationMethods()
Enum of the interpolation methods supported by FVInterpolationMethod objects.
Definition NSFVUtils.C:85
MooseEnum fvAdvectedInterpolationMethods()
Enum of the advected interpolation methods supported by FVInterpolationMethod objects.
Definition NSFVUtils.C:61

Member Data Documentation

◆ _add_rz_viscous_source

const bool WCNSFVFlowPhysicsBase::_add_rz_viscous_source
protectedinherited

Whether to automatically add the cylindrical viscous source term.

Definition at line 176 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource(), and WCNSFVFlowPhysicsBase::addAxisymmetricViscousSourceEnabled().

◆ _compressibility

const MooseEnum WCNSFVFlowPhysicsBase::_compressibility
protectedinherited

◆ _define_variables

bool NavierStokesPhysicsBase::_define_variables
protectedinherited

◆ _density_gravity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_gravity_name
protectedinherited

Name of the density material property used for gravity and Boussinesq terms.

Definition at line 200 of file WCNSFVFlowPhysicsBase.h.

Referenced by addFunctorMaterials(), WCNSFVFlowPhysics::addMomentumBoussinesqKernels(), addMomentumBoussinesqKernels(), and WCNSFVFlowPhysics::addMomentumGravityKernels().

◆ _density_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_name
protectedinherited

◆ _dynamic_viscosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_dynamic_viscosity_name
protectedinherited

◆ _flow_porosity_functor_name

MooseFunctorName WCNSFVFlowPhysicsBase::_flow_porosity_functor_name
protectedinherited

◆ _fluid_temperature_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_fluid_temperature_name
protectedinherited

◆ _flux_inlet_directions

std::vector<Point> WCNSFVFlowPhysicsBase::_flux_inlet_directions
protectedinherited

Direction of each flux inlet. Indexing based on the number of flux boundaries.

Definition at line 246 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addInletBC(), WCNSFVFlowPhysicsBase::getFluxInletDirections(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _flux_inlet_pps

std::vector<PostprocessorName> WCNSFVFlowPhysicsBase::_flux_inlet_pps
protectedinherited

Postprocessors describing the momentum inlet for each boundary. Indexing based on the number of flux boundaries.

Definition at line 244 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addInletBC(), and WCNSFVFlowPhysicsBase::getFluxInletPPs().

◆ _friction_blocks

std::vector<std::vector<SubdomainName> > WCNSFVFlowPhysicsBase::_friction_blocks
protectedinherited

Subdomains where we want to have volumetric friction.

Definition at line 221 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addMomentumFrictionKernels(), addMomentumFrictionKernels(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _friction_coeffs

std::vector<std::vector<std::string> > WCNSFVFlowPhysicsBase::_friction_coeffs
protectedinherited

The coefficients used for each item if friction type.

Definition at line 225 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addMomentumFrictionKernels(), addMomentumFrictionKernels(), and WCNSFVFlowPhysics::getLinearFrictionCoefName().

◆ _friction_types

std::vector<std::vector<std::string> > WCNSFVFlowPhysicsBase::_friction_types
protectedinherited

◆ _has_flow_equations

const bool WCNSFVFlowPhysicsBase::_has_flow_equations
protectedinherited

◆ _hydraulic_separators

const std::vector<BoundaryName> WCNSFVFlowPhysicsBase::_hydraulic_separators
protectedinherited

◆ _include_isotropic_viscous_stress

const bool WCNSFVFlowPhysicsBase::_include_isotropic_viscous_stress
protectedinherited

Whether to include the isotropic viscous stress contribution.

Definition at line 206 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::includeIsotropicStress().

◆ _include_symmetrized_viscous_stress

const bool WCNSFVFlowPhysicsBase::_include_symmetrized_viscous_stress
protectedinherited

Whether to include the symmetrized viscous stress contribution.

Definition at line 204 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::includeSymmetrizedViscousStress().

◆ _inlet_boundaries

const std::vector<BoundaryName> WCNSFVFlowPhysicsBase::_inlet_boundaries
protectedinherited

◆ _momentum_advection_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_advection_interpolation
protectedinherited

The momentum face interpolation method for being advected.

Definition at line 213 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addMomentumAdvectionKernels(), and WCNSFVFlowPhysicsBase::getMomentumAdvectionFaceInterpolationMethod().

◆ _momentum_face_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_face_interpolation
protectedinherited

The momentum face interpolation method for stress terms.

Definition at line 215 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::getMomentumFaceInterpolationMethod().

◆ _momentum_inlet_functors

std::map<BoundaryName, std::vector<MooseFunctorName> > WCNSFVFlowPhysicsBase::_momentum_inlet_functors
protectedinherited

Functors describing the momentum inlet for each boundary.

Definition at line 249 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addInletBC(), addInletBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _momentum_inlet_types

std::map<BoundaryName, MooseEnum> WCNSFVFlowPhysicsBase::_momentum_inlet_types
protectedinherited

◆ _momentum_outlet_types

std::map<BoundaryName, MooseEnum> WCNSFVFlowPhysicsBase::_momentum_outlet_types
protectedinherited

Momentum outlet boundary types.

Definition at line 239 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addOutletBC(), addOutletBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _momentum_wall_functors

std::map<BoundaryName, std::vector<MooseFunctorName> > WCNSFVFlowPhysicsBase::_momentum_wall_functors
protectedinherited

Functors describing the momentum for each wall boundary.

Definition at line 253 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addWallsBC(), addWallsBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _momentum_wall_types

std::map<BoundaryName, MooseEnum> WCNSFVFlowPhysicsBase::_momentum_wall_types
protectedinherited

Momentum wall boundary types.

Definition at line 241 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addWallsBC(), addWallsBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _non_orthogonal_correction

const bool WCNSLinearFVFlowPhysics::_non_orthogonal_correction
private

Whether to use the correction term for non-orthogonality.

Definition at line 75 of file WCNSLinearFVFlowPhysics.h.

Referenced by addMomentumFluxKernels(), and addPressureCorrectionKernels().

◆ _outlet_boundaries

const std::vector<BoundaryName> WCNSFVFlowPhysicsBase::_outlet_boundaries
protectedinherited

Boundaries with a flow outlet specified on them.

Definition at line 230 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::getOutletBoundaries(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _porosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_porosity_name
protectedinherited

◆ _porous_medium_treatment

const bool WCNSFVFlowPhysicsBase::_porous_medium_treatment
protectedinherited

◆ _pressure_functors

std::map<BoundaryName, MooseFunctorName> WCNSFVFlowPhysicsBase::_pressure_functors
protectedinherited

Functors describing the outlet pressure on each boundary.

Definition at line 251 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addOutletBC(), addOutletBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _pressure_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_pressure_name
protectedinherited

◆ _pressure_tag

const TagName WCNSFVFlowPhysicsBase::_pressure_tag = "p_tag"
protectedinherited

Name of the vector to hold pressure momentum equation contributions.

Definition at line 171 of file WCNSFVFlowPhysicsBase.h.

◆ _rc_uo_name

UserObjectName WCNSFVFlowPhysicsBase::_rc_uo_name
protectedinherited

name of the Rhie Chow user object

Definition at line 209 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::rhieChowUOName().

◆ _solve_for_dynamic_pressure

const bool WCNSFVFlowPhysicsBase::_solve_for_dynamic_pressure
protectedinherited

◆ _turbulence_physics

const WCNSFVTurbulencePhysicsBase* WCNSFVFlowPhysicsBase::_turbulence_physics
protectedinherited

Can be set to a coupled turbulence physics.

Definition at line 218 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::actOnAdditionalTasks(), addMomentumFluxKernels(), and WCNSFVFlowPhysicsBase::hasTurbulencePhysics().

◆ _velocity_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_velocity_interpolation
protectedinherited

The velocity face interpolation method for advecting other quantities.

Definition at line 211 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addMassKernels(), WCNSFVFlowPhysics::addMomentumAdvectionKernels(), and WCNSFVFlowPhysicsBase::getVelocityFaceInterpolationMethod().

◆ _velocity_names

const std::vector<std::string> WCNSFVFlowPhysicsBase::_velocity_names
protectedinherited

◆ _wall_boundaries

const std::vector<BoundaryName> WCNSFVFlowPhysicsBase::_wall_boundaries
protectedinherited

Boundaries which define a wall (slip/noslip/etc.)

Definition at line 232 of file WCNSFVFlowPhysicsBase.h.

Referenced by addWallsBC(), WCNSFVFlowPhysicsBase::getWallBoundaries(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().


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