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 | Friends | List of all members
WCNSFVFlowPhysicsBase Class Referenceabstract

Base class for Physics which create the Navier Stokes flow equations. More...

#include <WCNSFVFlowPhysicsBase.h>

Inheritance diagram for WCNSFVFlowPhysicsBase:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 WCNSFVFlowPhysicsBase (const InputParameters &parameters)
 
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 MooseEnumcompressibility () 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 MooseEnumgetVelocityFaceInterpolationMethod () const
 Get the face interpolation method for velocity.
 
const MooseEnumgetMomentumAdvectionFaceInterpolationMethod () const
 Get the face interpolation method for momentum in the advection term.
 
const MooseEnumgetMomentumFaceInterpolationMethod () 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.
 
virtual MooseFunctorName getLinearFrictionCoefName () const =0
 Get the name of the linear friction coefficient. Returns an empty string if no friction.
 
const UserObjectName & rhieChowUOName () const
 Return the name of the Rhie Chow user object.
 
unsigned short getNumberAlgebraicGhostingLayersNeeded () const override
 Return the number of algebraic ghosting layers needed.
 
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 TgetCoupledPhysics (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 ActionComponentgetActionComponent (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
 
virtual std::vector< UserObjectName > getSuppliedUserObjects () const
 
void timedAct ()
 
MooseObjectName uniqueActionName () const
 
const std::string & specificTaskName () const
 
const std::set< std::string > & getAllTasks () const
 
void appendTask (const std::string &task)
 
MooseAppgetMooseApp () 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 InputParametersparameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const TgetParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const TqueryParam (const std::string &name) const
 
const TgetRenamedParam (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
 
PerfGraphperfGraph ()
 
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 addSolverVariables () override=0
 
virtual void addInitialConditions () override
 
virtual void addFVKernels () override=0
 
virtual void addFVBCs () override
 
virtual void addMaterials () override
 
virtual void addUserObjects () override=0
 
virtual void addPostprocessors () override
 
virtual void addMomentumTimeKernels ()=0
 Functions adding kernels for the flow momentum equations If the material properties are not constant, these can be used for weakly-compressible simulations (except the Boussinesq kernel) as well.
 
virtual void addMomentumPressureKernels ()=0
 
virtual void addMomentumGravityKernels ()=0
 
virtual void addMomentumFrictionKernels ()=0
 
virtual void addMomentumBoussinesqKernels ()=0
 
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.
 
virtual void addInletBC ()=0
 Functions adding boundary conditions for the flow simulation.
 
virtual void addOutletBC ()=0
 
virtual void addWallsBC ()=0
 
virtual void addSeparatorBC ()=0
 
virtual bool hasForchheimerFriction () const =0
 Return whether a Forchheimer friction model is in use.
 
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.
 
virtual void addRhieChowUserObjects ()=0
 Function which adds the RhieChow interpolator user objects for weakly and incompressible formulations.
 
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 WCNSFVTurbulencePhysicsBasegetCoupledTurbulencePhysics () 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
 
FactorygetFactory ()
 
FactorygetFactory () const
 
virtual FEProblemBasegetProblem ()
 
virtual const FEProblemBasegetProblem () 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)
 
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< SubdomainIDgetSubdomainIDs (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 TgetMeshProperty (const std::string &data_name, const std::string &prefix)
 
const TgetMeshProperty (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 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 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 addFVInterpolationMethods ()
 
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 addFunctorMaterials ()
 
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 RestartableDataValuegetMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const
 
void mooseErrorInternal (Args &&... args) const
 
const TforwardGetParam (const std::string &param_name) const
 
const InputParametersforwardParameters () 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

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
 

Friends

class WCNSFVTurbulencePhysics
 

Detailed Description

Base class for Physics which create the Navier Stokes flow equations.

Definition at line 29 of file WCNSFVFlowPhysicsBase.h.

Constructor & Destructor Documentation

◆ WCNSFVFlowPhysicsBase()

WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase ( const InputParameters parameters)

Definition at line 119 of file WCNSFVFlowPhysicsBase.C.

121 _has_flow_equations(getParam<bool>("add_flow_equations")),
122 _add_rz_viscous_source(getParam<bool>("add_rz_viscous_source")),
123 _compressibility(getParam<MooseEnum>("compressibility")),
124 _solve_for_dynamic_pressure(getParam<bool>("solve_for_dynamic_pressure")),
125 _porous_medium_treatment(getParam<bool>("porous_medium_treatment")),
126 _porosity_name(getParam<MooseFunctorName>("porosity")),
129 isParamValid("velocity_variable")
130 ? getParam<std::vector<std::string>>("velocity_variable")
132 ? std::vector<std::string>(NS::superficial_velocity_vector,
133 NS::superficial_velocity_vector + 3)
134 : std::vector<std::string>(NS::velocity_vector, NS::velocity_vector + 3))),
135 _pressure_name(isParamValid("pressure_variable")
136 ? getParam<NonlinearVariableName>("pressure_variable")
137 : NS::pressure),
138 _fluid_temperature_name(isParamValid("fluid_temperature_variable")
139 ? getParam<NonlinearVariableName>("fluid_temperature_variable")
140 : NS::T_fluid),
141 _density_name(getParam<MooseFunctorName>("density")),
142 _density_gravity_name(isParamValid("density_for_gravity_terms")
143 ? getParam<MooseFunctorName>("density_for_gravity_terms")
144 : getParam<MooseFunctorName>("density")),
145 _dynamic_viscosity_name(getParam<MooseFunctorName>("dynamic_viscosity")),
146 _include_symmetrized_viscous_stress(getParam<bool>("include_symmetrized_viscous_stress")),
147 _include_isotropic_viscous_stress(getParam<bool>("include_isotropic_viscous_stress")),
148 _rc_uo_name(isParamValid("rhie_chow_uo_name")
149 ? getParam<UserObjectName>("rhie_chow_uo_name")
150 : (_porous_medium_treatment ? "pins_rhie_chow_interpolator"
151 : "ins_rhie_chow_interpolator")),
152 _velocity_interpolation(getParam<MooseEnum>("velocity_interpolation")),
153 _momentum_advection_interpolation(getParam<MooseEnum>("momentum_advection_interpolation")),
154 _momentum_face_interpolation(getParam<MooseEnum>("momentum_face_interpolation")),
155 _friction_blocks(getParam<std::vector<std::vector<SubdomainName>>>("friction_blocks")),
156 _friction_types(getParam<std::vector<std::vector<std::string>>>("friction_types")),
157 _friction_coeffs(getParam<std::vector<std::vector<std::string>>>("friction_coeffs")),
158 _inlet_boundaries(getParam<std::vector<BoundaryName>>("inlet_boundaries")),
159 _outlet_boundaries(getParam<std::vector<BoundaryName>>("outlet_boundaries")),
160 _wall_boundaries(getParam<std::vector<BoundaryName>>("wall_boundaries")),
161 _hydraulic_separators(getParam<std::vector<BoundaryName>>("hydraulic_separator_sidesets")),
162 _flux_inlet_pps(getParam<std::vector<PostprocessorName>>("flux_inlet_pps")),
163 _flux_inlet_directions(getParam<std::vector<Point>>("flux_inlet_directions"))
164{
165 // Inlet boundary parameter checking
166 checkSecondParamSetOnlyIfFirstOneSet("flux_inlet_pps", "flux_inlet_directions");
167 if (_flux_inlet_directions.size())
168 checkVectorParamsSameLengthIfSet<PostprocessorName, Point>("flux_inlet_pps",
169 "flux_inlet_directions");
170
171 // Boussinesq parameters checks
172 checkSecondParamSetOnlyIfFirstOneTrue("boussinesq_approximation", "ref_temperature");
173 checkSecondParamSetOnlyIfFirstOneSet("gravity", "boussinesq_approximation");
174
175 // Dynamic pressure parameter checks
176 if (_compressibility != "incompressible" && _solve_for_dynamic_pressure)
177 paramError("compressibility",
178 "Solving for dynamic pressure is only implemented for incompressible flow");
179
180 // Boundary parameters checking
181 checkVectorParamAndMultiMooseEnumLength<BoundaryName>("inlet_boundaries", "momentum_inlet_types");
182 checkVectorParamAndMultiMooseEnumLength<BoundaryName>("outlet_boundaries",
183 "momentum_outlet_types");
184 checkVectorParamAndMultiMooseEnumLength<BoundaryName>("wall_boundaries", "momentum_wall_types");
186 std::vector<MooseFunctorName>,
187 PostprocessorName>(
188 "inlet_boundaries", "momentum_inlet_functors", "flux_inlet_pps");
189 checkVectorParamsNoOverlap<BoundaryName>(
190 {"inlet_boundaries", "outlet_boundaries", "wall_boundaries"});
191
192 // Porous media parameters
193 checkSecondParamSetOnlyIfFirstOneTrue("porous_medium_treatment", "porosity");
194
196 for (const auto & name : NS::velocity_vector)
197 {
198 const auto & it = std::find(_velocity_names.begin(), _velocity_names.end(), name);
199 if (it != _velocity_names.end())
200 paramError("velocity_variable",
201 "For porous medium simulations, functor name " + *it +
202 " is already reserved for the automatically-computed interstitial velocity. "
203 "Please choose another name for your external velocity variable!");
204 }
205
206 // Friction parameter checks
207 if (_friction_blocks.size())
208 checkVectorParamsSameLength<std::vector<SubdomainName>, std::vector<std::string>>(
209 "friction_blocks", "friction_types");
210 checkTwoDVectorParamsSameLength<std::string, std::string>("friction_types", "friction_coeffs");
211
212 // Create maps for boundary-restricted parameters
213 _momentum_inlet_types = Moose::createMapFromVectorAndMultiMooseEnum<BoundaryName>(
214 _inlet_boundaries, getParam<MultiMooseEnum>("momentum_inlet_types"));
215 _momentum_outlet_types = Moose::createMapFromVectorAndMultiMooseEnum<BoundaryName>(
216 _outlet_boundaries, getParam<MultiMooseEnum>("momentum_outlet_types"));
217 _momentum_wall_types = Moose::createMapFromVectorAndMultiMooseEnum<BoundaryName>(
218 _wall_boundaries, getParam<MultiMooseEnum>("momentum_wall_types"));
219 if (isParamSetByUser("momentum_inlet_functors"))
220 {
221 // Not all inlet boundary types require the specification of an inlet functor
222 std::vector<BoundaryName> inlet_boundaries_with_functors;
223 for (const auto & boundary : _inlet_boundaries)
224 if (libmesh_map_find(_momentum_inlet_types, boundary) == "fixed-velocity" ||
225 libmesh_map_find(_momentum_inlet_types, boundary) == "fixed-pressure")
226 inlet_boundaries_with_functors.push_back(boundary);
228 Moose::createMapFromVectors<BoundaryName, std::vector<MooseFunctorName>>(
229 inlet_boundaries_with_functors,
230 getParam<std::vector<std::vector<MooseFunctorName>>>("momentum_inlet_functors"));
231 }
232 if (isParamSetByUser("pressure_functors"))
233 {
234 // Not all outlet boundary types require the specification of an inlet functor
235 std::vector<BoundaryName> outlet_boundaries_with_functors;
236 for (const auto & boundary : _outlet_boundaries)
237 if (libmesh_map_find(_momentum_outlet_types, boundary) == "fixed-pressure-zero-gradient" ||
238 libmesh_map_find(_momentum_outlet_types, boundary) == "fixed-pressure")
239 outlet_boundaries_with_functors.push_back(boundary);
240 const auto & pressure_functors = getParam<std::vector<MooseFunctorName>>("pressure_functors");
241 if (outlet_boundaries_with_functors.size() != pressure_functors.size())
242 paramError("pressure_functors",
243 "Size (" + std::to_string(pressure_functors.size()) +
244 ") is not the same as the number of pressure outlet boundaries in "
245 "'fixed-pressure/fixed-pressure-zero-gradient' (size " +
246 std::to_string(outlet_boundaries_with_functors.size()) + ")");
247 _pressure_functors = Moose::createMapFromVectors<BoundaryName, MooseFunctorName>(
248 outlet_boundaries_with_functors, pressure_functors);
249 }
250
251 if (isParamSetByUser("momentum_wall_functors"))
252 {
253 // Not all wall boundary types require the specification of an inlet functor
254 std::vector<BoundaryName> wall_boundaries_with_functors;
255 for (const auto & boundary : _wall_boundaries)
256 if (libmesh_map_find(_momentum_wall_types, boundary) == "noslip")
257 wall_boundaries_with_functors.push_back(boundary);
258 const auto & momentum_wall_functors =
259 getParam<std::vector<std::vector<MooseFunctorName>>>("momentum_wall_functors");
260 if (wall_boundaries_with_functors.size() != momentum_wall_functors.size())
261 paramError("momentum_wall_functors",
262 "Size (" + std::to_string(momentum_wall_functors.size()) +
263 ") is not the same as the number of momentum_wall wall boundaries with "
264 "no-slip boundary conditions ' (size " +
265 std::to_string(wall_boundaries_with_functors.size()) + ")");
266
268 Moose::createMapFromVectors<BoundaryName, std::vector<MooseFunctorName>>(
269 wall_boundaries_with_functors, momentum_wall_functors);
270 }
271
272 addRequiredPhysicsTask("add_geometric_rm");
273 addRequiredPhysicsTask("add_variables_physics");
274 addRequiredPhysicsTask("add_ics_physics");
275 addRequiredPhysicsTask("add_materials_physics");
276 addRequiredPhysicsTask("add_user_object");
277 addRequiredPhysicsTask("add_postprocessor");
278 addRequiredPhysicsTask("add_corrector");
279 addRequiredPhysicsTask("get_turbulence_physics");
280}
void checkVectorParamLengthSameAsCombinedOthers(const std::string &param1, const std::string &param2, const std::string &param3) const
void checkSecondParamSetOnlyIfFirstOneTrue(const std::string &param1, const std::string &param2) const
void checkSecondParamSetOnlyIfFirstOneSet(const std::string &param1, const std::string &param2) const
const InputParameters & parameters() const
const std::string & name() const
void paramError(const std::string &param, Args... args) const
bool isParamSetByUser(const std::string &name) const
const T & getParam(const std::string &name) const
bool isParamValid(const std::string &name) const
Base class to hold common parameters and utilities between all the weakly compressible Navier Stokes-...
bool _define_variables
Whether to define variables if they do not exist.
void addRequiredPhysicsTask(const std::string &task)
std::map< BoundaryName, MooseEnum > _momentum_outlet_types
Momentum outlet boundary types.
const MooseEnum _velocity_interpolation
The velocity face interpolation method for advecting other quantities.
std::vector< std::vector< SubdomainName > > _friction_blocks
Subdomains where we want to have volumetric friction.
std::map< BoundaryName, std::vector< MooseFunctorName > > _momentum_wall_functors
Functors describing the momentum for each wall boundary.
std::map< BoundaryName, MooseEnum > _momentum_inlet_types
Momentum inlet boundary types.
std::map< BoundaryName, MooseFunctorName > _pressure_functors
Functors describing the outlet pressure on each boundary.
std::vector< PostprocessorName > _flux_inlet_pps
Postprocessors describing the momentum inlet for each boundary. Indexing based on the number of flux ...
const MooseEnum _momentum_advection_interpolation
The momentum face interpolation method for being advected.
std::vector< std::vector< std::string > > _friction_types
The friction correlation types used for each block.
const MooseEnum _compressibility
Compressibility type, can be compressible, incompressible or weakly-compressible.
const std::vector< std::string > _velocity_names
Velocity names.
const MooseEnum _momentum_face_interpolation
The momentum face interpolation method for stress terms.
const std::vector< BoundaryName > _wall_boundaries
Boundaries which define a wall (slip/noslip/etc.)
MooseFunctorName _flow_porosity_functor_name
Name of the porosity functor for the flow equations (if smoothed)
const NonlinearVariableName _pressure_name
Pressure name.
const bool _has_flow_equations
Boolean to keep track of whether the flow equations should be created.
const NonlinearVariableName _fluid_temperature_name
Fluid temperature name.
const bool _porous_medium_treatment
Whether to use the porous medium treatment.
std::map< BoundaryName, std::vector< MooseFunctorName > > _momentum_inlet_functors
Functors describing the momentum inlet for each boundary.
const std::vector< BoundaryName > _outlet_boundaries
Boundaries with a flow outlet specified on them.
std::vector< Point > _flux_inlet_directions
Direction of each flux inlet. Indexing based on the number of flux boundaries.
const MooseFunctorName _dynamic_viscosity_name
Name of the dynamic viscosity material property.
std::map< BoundaryName, MooseEnum > _momentum_wall_types
Momentum wall boundary types.
const std::vector< BoundaryName > _hydraulic_separators
Hydraulic separator boundaries.
const bool _solve_for_dynamic_pressure
Whether we are solving for the total or dynamic pressure.
UserObjectName _rc_uo_name
name of the Rhie Chow user object
const bool _add_rz_viscous_source
Whether to automatically add the cylindrical viscous source term.
const bool _include_isotropic_viscous_stress
Whether to include the isotropic viscous stress contribution.
const bool _include_symmetrized_viscous_stress
Whether to include the symmetrized viscous stress contribution.
const MooseFunctorName _density_name
Name of the density material property.
const MooseFunctorName _porosity_name
Name of the porosity functor.
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 MooseFunctorName _density_gravity_name
Name of the density material property used for gravity and Boussinesq terms.
const std::string velocity_vector[3]
Definition NS.h:50
if(subdm)

Member Function Documentation

◆ actOnAdditionalTasks()

void WCNSFVFlowPhysicsBase::actOnAdditionalTasks ( )
overrideprotectedvirtual

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

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.

Referenced by WCNSFVFlowPhysics::addFVKernels().

◆ addAxisymmetricViscousSourceEnabled()

bool WCNSFVFlowPhysicsBase::addAxisymmetricViscousSourceEnabled ( ) const
inline

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

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

◆ addFluidPropertiesFunctorMaterial()

void WCNSFVFlowPhysicsBase::addFluidPropertiesFunctorMaterial ( )
protected

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
Factory & getFactory()
void assignBlocks(InputParameters &params, const std::vector< SubdomainName > &blocks) const
std::string prefix() const
std::vector< SubdomainName > _blocks
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:147
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 addMaterials().

◆ 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 WCNSLinearFVFlowPhysics::addFVInterpolationMethods(), WCNSLinearFVScalarTransportPhysics::addFVInterpolationMethods(), and WCNSLinearFVTurbulencePhysics::addFVInterpolationMethods().

◆ addFVBCs()

void WCNSFVFlowPhysicsBase::addFVBCs ( )
overrideprotectedvirtual

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

◆ addFVKernels()

virtual void WCNSFVFlowPhysicsBase::addFVKernels ( )
overrideprotectedpure virtual

Reimplemented from PhysicsBase.

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

◆ addInitialConditions()

void WCNSFVFlowPhysicsBase::addInitialConditions ( )
overrideprotectedvirtual

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
bool shouldCreateIC(const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool ic_is_default_ic, const bool error_if_already_defined) const
unsigned int dimension() const
IntRange< T > make_range(T beg, T end)

◆ addInletBC()

virtual void WCNSFVFlowPhysicsBase::addInletBC ( )
protectedpure virtual

Functions adding boundary conditions for the flow simulation.

These are used for weakly-compressible simulations as well.

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

Referenced by addFVBCs().

◆ addMaterials()

void WCNSFVFlowPhysicsBase::addMaterials ( )
overrideprotectedvirtual

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

virtual void WCNSFVFlowPhysicsBase::addMomentumBoussinesqKernels ( )
protectedpure virtual

◆ addMomentumFrictionKernels()

virtual void WCNSFVFlowPhysicsBase::addMomentumFrictionKernels ( )
protectedpure virtual

◆ addMomentumGravityKernels()

virtual void WCNSFVFlowPhysicsBase::addMomentumGravityKernels ( )
protectedpure virtual

◆ addMomentumPressureKernels()

virtual void WCNSFVFlowPhysicsBase::addMomentumPressureKernels ( )
protectedpure virtual

◆ addMomentumTimeKernels()

virtual void WCNSFVFlowPhysicsBase::addMomentumTimeKernels ( )
protectedpure virtual

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

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

◆ addNonPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addNonPorousMediumSpeedMaterial ( )
protected

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

◆ addOutletBC()

virtual void WCNSFVFlowPhysicsBase::addOutletBC ( )
protectedpure virtual

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

Referenced by addFVBCs().

◆ addPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addPorousMediumSpeedMaterial ( )
protected

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}
static const std::string superficial_velocity_vector[3]
Definition NS.h:55

Referenced by addMaterials().

◆ addPostprocessors()

void WCNSFVFlowPhysicsBase::addPostprocessors ( )
overrideprotectedvirtual

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

◆ addRhieChowUserObjects()

virtual void WCNSFVFlowPhysicsBase::addRhieChowUserObjects ( )
protectedpure virtual

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

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

◆ addSeparatorBC()

virtual void WCNSFVFlowPhysicsBase::addSeparatorBC ( )
protectedpure virtual

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

Referenced by addFVBCs().

◆ addSolverVariables()

virtual void WCNSFVFlowPhysicsBase::addSolverVariables ( )
overrideprotectedpure virtual

Reimplemented from PhysicsBase.

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

◆ addUserObjects()

virtual void WCNSFVFlowPhysicsBase::addUserObjects ( )
overrideprotectedpure virtual

Reimplemented from PhysicsBase.

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

◆ addWallsBC()

virtual void WCNSFVFlowPhysicsBase::addWallsBC ( )
protectedpure virtual

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

Referenced by addFVBCs().

◆ compressibility()

const MooseEnum & WCNSFVFlowPhysicsBase::compressibility ( ) const
inline

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
inline

◆ dynamicViscosityName()

const MooseFunctorName & WCNSFVFlowPhysicsBase::dynamicViscosityName ( ) const
inline

◆ 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
protected

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}
bool allMeshBlocks(const std::vector< SubdomainName > &blocks) const
MeshBase & mesh
std::string stringify(const T &t)

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

◆ getCoupledTurbulencePhysics()

const WCNSFVTurbulencePhysicsBase * WCNSFVFlowPhysicsBase::getCoupledTurbulencePhysics ( ) const
protected

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

◆ getFlowVariableName()

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

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}
void mooseError(Args &&... args) const
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
inline

◆ getFluxInletDirections()

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

Get the inlet direction if using a flux inlet.

Definition at line 89 of file WCNSFVFlowPhysicsBase.h.

89{ return _flux_inlet_directions; }

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

◆ getFluxInletPPs()

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

Get the inlet flux postprocessor if using a flux inlet.

Definition at line 91 of file WCNSFVFlowPhysicsBase.h.

91{ return _flux_inlet_pps; }

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

◆ getHydraulicSeparators()

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

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
inline

◆ getLinearFrictionCoefName()

virtual MooseFunctorName WCNSFVFlowPhysicsBase::getLinearFrictionCoefName ( ) const
pure virtual

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

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

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

◆ getMomentumAdvectionFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getMomentumAdvectionFaceInterpolationMethod ( ) const
inline

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

Definition at line 65 of file WCNSFVFlowPhysicsBase.h.

66 {
68 }

◆ getMomentumFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getMomentumFaceInterpolationMethod ( ) const
inline

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

Definition at line 70 of file WCNSFVFlowPhysicsBase.h.

71 {
73 }

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

◆ getNumberAlgebraicGhostingLayersNeeded()

unsigned short WCNSFVFlowPhysicsBase::getNumberAlgebraicGhostingLayersNeeded ( ) const
overridevirtual

◆ getOutletBoundaries()

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

◆ getPorosityFunctorName()

MooseFunctorName WCNSFVFlowPhysicsBase::getPorosityFunctorName ( const bool  smoothed) const

◆ getPressureName()

const NonlinearVariableName & WCNSFVFlowPhysicsBase::getPressureName ( ) const
inline

◆ getVelocityFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getVelocityFaceInterpolationMethod ( ) const
inline

Get the face interpolation method for velocity.

Definition at line 63 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm().

◆ getVelocityNames()

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

◆ getWallBoundaries()

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

◆ gravityVector()

RealVectorValue WCNSFVFlowPhysicsBase::gravityVector ( ) const
inline

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
inline

◆ hasForchheimerFriction()

virtual bool WCNSFVFlowPhysicsBase::hasForchheimerFriction ( ) const
protectedpure virtual

Return whether a Forchheimer friction model is in use.

Implemented in WCNSFVFlowPhysics, and WCNSLinearFVFlowPhysics.

Referenced by addMaterials().

◆ hasTurbulencePhysics()

bool WCNSFVFlowPhysicsBase::hasTurbulencePhysics ( ) const
inlineprotected

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
inline

Whether to include the isotropic viscous stress contribution.

Definition at line 43 of file WCNSFVFlowPhysicsBase.h.

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

◆ includeSymmetrizedViscousStress()

bool WCNSFVFlowPhysicsBase::includeSymmetrizedViscousStress ( ) const
inline

◆ initializePhysicsAdditional()

void WCNSFVFlowPhysicsBase::initializePhysicsAdditional ( )
overrideprotectedvirtual

Reimplemented from PhysicsBase.

Reimplemented in WCNSLinearFVFlowPhysics.

Definition at line 283 of file WCNSFVFlowPhysicsBase.C.

284{
285 getProblem().needFV();
286}
virtual void needFV() override

Referenced by WCNSLinearFVFlowPhysics::initializePhysicsAdditional().

◆ inletBoundaryType()

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

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:216

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC().

◆ porousMediumTreatment()

bool WCNSFVFlowPhysicsBase::porousMediumTreatment ( ) const
inline

Return whether a porous medium treatment is applied.

Definition at line 55 of file WCNSFVFlowPhysicsBase.h.

◆ rhieChowUOName()

const UserObjectName & WCNSFVFlowPhysicsBase::rhieChowUOName ( ) const

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}

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(), WCNSLinearFVFlowPhysics::addMomentumFluxKernels(), WCNSFVFlowPhysics::addMomentumFrictionKernels(), WCNSFVFlowPhysics::addMomentumGravityKernels(), WCNSFVFlowPhysics::addMomentumPressureKernels(), WCNSFVFlowPhysics::addMomentumTimeKernels(), WCNSFVFlowPhysics::addMomentumViscousDissipationKernels(), WCNSFVFlowPhysics::addOutletBC(), WCNSFVTwoPhaseMixturePhysics::addPhaseDriftFluxTerm(), WCNSLinearFVTwoPhaseMixturePhysics::addPhaseDriftFluxTerm(), WCNSLinearFVFlowPhysics::addRhieChowUserObjects(), WCNSFVScalarTransportPhysics::addScalarAdvectionKernels(), WCNSLinearFVScalarTransportPhysics::addScalarAdvectionKernels(), WCNSFVFlowPhysics::addSeparatorBC(), WCNSFVFlowPhysics::addWallsBC(), WCNSFVFlowPhysics::getSuppliedUserObjects(), and WCNSLinearFVFlowPhysics::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 WCNSFVFlowPhysicsBase::validParams ( )
static

Definition at line 19 of file WCNSFVFlowPhysicsBase.C.

20{
22 params.addClassDescription("Base class for Physics defining the Navier Stokes flow equations");
23
25 params.addParam<bool>("add_flow_equations",
26 true,
27 "Whether to add the flow equations. This parameter is not necessary when "
28 "using the Physics syntax");
29
30 // We pull in parameters from various flow objects. This helps make sure the parameters are
31 // spelled the same way and match the evolution of other objects.
32 // If we remove these objects, or change their parameters, these parameters should be updated
33 // Downstream actions must either implement all these options, or redefine the parameter with
34 // a restricted MooseEnum, or place an error in the constructor for unsupported configurations
35 // We mostly pull the boundary parameters from NSFV Action
36
38 params.addParam<bool>("include_deviatoric_stress",
39 false,
40 "Whether to include the symmetrized viscous stress contribution "
41 "(grad(u)+grad(u)^T).");
42 params.deprecateParam(
43 "include_deviatoric_stress", "include_symmetrized_viscous_stress", "12/31/2025");
44 params.addParam<bool>("include_isotropic_viscous_stress",
45 false,
46 "Whether to add the isotropic -(2/3) mu div(u) I contribution to the "
47 "viscous stress (requires complete expansion of the velocity gradient).");
48 params.addParam<bool>("add_rz_viscous_source",
49 true,
50 "When true, automatically adds INSFVMomentumViscousSourceRZ to the radial "
51 "momentum equation on RZ blocks. Disable if this term should be omitted.");
52
53 // Momentum boundary conditions are important for advection problems as well
55
56 // Convenient for sharing the executioner block with a non-Physics syntax
57 // or for defining the RhieChow user object outside the Physics
58 params.addParam<UserObjectName>(
59 "rhie_chow_uo_name",
60 "Name of the Rhie Chow user object. Defaults to 'pins_rhie_chow_interpolator' for porous "
61 "flow and 'ins_rhie_chow_interpolator' for non-porous.");
62
63 // Specify the weakly compressible boundary flux information. They are used for specifying in flux
64 // boundary conditions for advection physics in WCNSFV
66 params.addParam<std::vector<std::vector<MooseFunctorName>>>(
67 "momentum_wall_functors",
68 {},
69 "Functors for each component of the velocity value on walls. This is only necessary for the "
70 "fixed-velocity momentum wall types.");
71
72 // Most downstream physics implementations are valid for porous media too
73 // If yours is not, please remember to disable the 'porous_medium_treatment' parameter
74 params.transferParam<bool>(NSFVBase::validParams(), "porous_medium_treatment");
75 params.transferParam<MooseFunctorName>(NSFVBase::validParams(), "porosity");
76
77 // Initialization parameters
78 params.transferParam<std::vector<FunctionName>>(NSFVBase::validParams(), "initial_velocity");
79 params.transferParam<FunctionName>(NSFVBase::validParams(), "initial_pressure");
80
81 // Spatial discretization scheme
82 // Specify the numerical schemes for interpolations of velocity and pressure
83 params.transferParam<MooseEnum>(NSFVBase::validParams(), "velocity_interpolation");
84 params.transferParam<MooseEnum>(NSFVBase::validParams(), "momentum_advection_interpolation");
85 params.transferParam<MooseEnum>(NSFVBase::validParams(), "momentum_face_interpolation");
86 params.transferParam<bool>(NSFVBase::validParams(), "momentum_two_term_bc_expansion");
87 params.transferParam<bool>(NSFVBase::validParams(), "pressure_two_term_bc_expansion");
88 MooseEnum coeff_interp_method("average harmonic", "harmonic");
89 params.addParam<MooseEnum>("mu_interp_method",
90 coeff_interp_method,
91 "Switch that can select face interpolation method for the viscosity.");
92
93 // Fluid properties
94 params.addParam<UserObjectName>(NS::fluid, "Fluid properties userobject");
95 params.addParam<FunctionName>(
96 "mu_rampdown", 1, "A function describing a ramp down of viscosity over time");
97 params.addParamNamesToGroup(NS::fluid + " mu_rampdown", "Material properties");
98
99 // Parameter groups
100 params.addParamNamesToGroup(
101 "velocity_variable pressure_variable initial_pressure initial_velocity", "Variables");
102 params.addParamNamesToGroup("density dynamic_viscosity", "Material properties");
103 params.addParamNamesToGroup("inlet_boundaries momentum_inlet_types momentum_inlet_functors",
104 "Inlet boundary conditions");
105 params.addParamNamesToGroup("outlet_boundaries momentum_outlet_types pressure_functors",
106 "Outlet boundary conditions");
107 params.addParamNamesToGroup("wall_boundaries momentum_wall_types momentum_wall_functors",
108 "Wall boundary conditions");
109 params.addParamNamesToGroup(
110 "include_symmetrized_viscous_stress include_isotropic_viscous_stress velocity_interpolation "
111 "momentum_advection_interpolation momentum_two_term_bc_expansion "
112 "pressure_two_term_bc_expansion mu_interp_method momentum_face_interpolation",
113 "Numerical scheme");
114 params.addParamNamesToGroup("thermal_expansion", "Gravity treatment");
115
116 return params;
117}
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 commonMomentumEquationParams()
Definition NSFVBase.C:55
static InputParameters commonMomentumBoundaryTypesParams()
Definition NSFVBase.C:144
static InputParameters validParams()
Definition NSFVBase.C:371
static InputParameters commonMomentumBoundaryFluxesParams()
Definition NSFVBase.C:168
static InputParameters commonNavierStokesFlowParams()
Definition NSFVBase.C:17
static InputParameters validParams()

Referenced by WCNSFVFlowPhysics::validParams(), and WCNSLinearFVFlowPhysics::validParams().

Friends And Related Symbol Documentation

◆ WCNSFVTurbulencePhysics

friend class WCNSFVTurbulencePhysics
friend

Definition at line 255 of file WCNSFVFlowPhysicsBase.h.

Member Data Documentation

◆ _add_rz_viscous_source

const bool WCNSFVFlowPhysicsBase::_add_rz_viscous_source
protected

Whether to automatically add the cylindrical viscous source term.

Definition at line 176 of file WCNSFVFlowPhysicsBase.h.

Referenced by addAxisymmetricViscousSource(), and addAxisymmetricViscousSourceEnabled().

◆ _compressibility

const MooseEnum WCNSFVFlowPhysicsBase::_compressibility
protected

◆ _define_variables

bool NavierStokesPhysicsBase::_define_variables
protectedinherited

◆ _density_gravity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_gravity_name
protected

◆ _density_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_name
protected

◆ _dynamic_viscosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_dynamic_viscosity_name
protected

◆ _flow_porosity_functor_name

MooseFunctorName WCNSFVFlowPhysicsBase::_flow_porosity_functor_name
protected

◆ _fluid_temperature_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_fluid_temperature_name
protected

◆ _flux_inlet_directions

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

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(), getFluxInletDirections(), and WCNSFVFlowPhysicsBase().

◆ _flux_inlet_pps

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

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

◆ _friction_blocks

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

Subdomains where we want to have volumetric friction.

Definition at line 221 of file WCNSFVFlowPhysicsBase.h.

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

◆ _friction_coeffs

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

◆ _friction_types

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

◆ _has_flow_equations

const bool WCNSFVFlowPhysicsBase::_has_flow_equations
protected

◆ _hydraulic_separators

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

◆ _include_isotropic_viscous_stress

const bool WCNSFVFlowPhysicsBase::_include_isotropic_viscous_stress
protected

Whether to include the isotropic viscous stress contribution.

Definition at line 206 of file WCNSFVFlowPhysicsBase.h.

Referenced by includeIsotropicStress().

◆ _include_symmetrized_viscous_stress

const bool WCNSFVFlowPhysicsBase::_include_symmetrized_viscous_stress
protected

Whether to include the symmetrized viscous stress contribution.

Definition at line 204 of file WCNSFVFlowPhysicsBase.h.

Referenced by includeSymmetrizedViscousStress().

◆ _inlet_boundaries

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

Boundaries with a flow inlet specified on them.

Definition at line 228 of file WCNSFVFlowPhysicsBase.h.

Referenced by addPostprocessors(), getInletBoundaries(), and WCNSFVFlowPhysicsBase().

◆ _momentum_advection_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_advection_interpolation
protected

The momentum face interpolation method for being advected.

Definition at line 213 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_face_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_face_interpolation
protected

The momentum face interpolation method for stress terms.

Definition at line 215 of file WCNSFVFlowPhysicsBase.h.

Referenced by getMomentumFaceInterpolationMethod().

◆ _momentum_inlet_functors

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

Functors describing the momentum inlet for each boundary.

Definition at line 249 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_inlet_types

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

◆ _momentum_outlet_types

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

Momentum outlet boundary types.

Definition at line 239 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_wall_functors

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

Functors describing the momentum for each wall boundary.

Definition at line 253 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_wall_types

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

Momentum wall boundary types.

Definition at line 241 of file WCNSFVFlowPhysicsBase.h.

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

◆ _outlet_boundaries

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

Boundaries with a flow outlet specified on them.

Definition at line 230 of file WCNSFVFlowPhysicsBase.h.

Referenced by getOutletBoundaries(), and WCNSFVFlowPhysicsBase().

◆ _porosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_porosity_name
protected

◆ _porous_medium_treatment

const bool WCNSFVFlowPhysicsBase::_porous_medium_treatment
protected

◆ _pressure_functors

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

Functors describing the outlet pressure on each boundary.

Definition at line 251 of file WCNSFVFlowPhysicsBase.h.

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

◆ _pressure_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_pressure_name
protected

◆ _pressure_tag

const TagName WCNSFVFlowPhysicsBase::_pressure_tag = "p_tag"
protected

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
protected

name of the Rhie Chow user object

Definition at line 209 of file WCNSFVFlowPhysicsBase.h.

Referenced by rhieChowUOName().

◆ _solve_for_dynamic_pressure

const bool WCNSFVFlowPhysicsBase::_solve_for_dynamic_pressure
protected

◆ _turbulence_physics

const WCNSFVTurbulencePhysicsBase* WCNSFVFlowPhysicsBase::_turbulence_physics
protected

Can be set to a coupled turbulence physics.

Definition at line 218 of file WCNSFVFlowPhysicsBase.h.

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

◆ _velocity_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_velocity_interpolation
protected

The velocity face interpolation method for advecting other quantities.

Definition at line 211 of file WCNSFVFlowPhysicsBase.h.

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

◆ _velocity_names

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

◆ _wall_boundaries

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

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

Definition at line 232 of file WCNSFVFlowPhysicsBase.h.

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


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