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)

bool	hasFlowEquations () const
	Whether the physics is actually creating the flow equations. More...

const std::vector< std::string > &	getVelocityNames () const
	To interface with other Physics. More...

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. More...

bool	porousMediumTreatment () const
	Return whether a porous medium treatment is applied. More...

RealVectorValue	gravityVector () const
	Return the gravity vector. More...

const MooseFunctorName &	densityName () const
	Return the name of the density functor. More...

const MooseFunctorName &	dynamicViscosityName () const
	Return the name of the dynamic viscosity functor. More...

const MooseEnum &	getVelocityFaceInterpolationMethod () const
	Get the face interpolation method for velocity. More...

const MooseEnum &	getMomentumAdvectionFaceInterpolationMethod () const
	Get the face interpolation method for momentum in the advection term. More...

const MooseEnum &	getMomentumFaceInterpolationMethod () const
	Get the face interpolation method for momentum (mostly used in the stress terms) More...

const std::vector< BoundaryName > &	getInletBoundaries () const
	Get the inlet boundaries. More...

const std::vector< BoundaryName > &	getOutletBoundaries () const
	Get the outlet boundaries. More...

const std::vector< BoundaryName > &	getWallBoundaries () const
	Get the wall boundaries. More...

const std::vector< BoundaryName > &	getHydraulicSeparators () const
	Get the hydraulic separator boundaries. More...

const std::vector< Point > &	getFluxInletDirections () const
	Get the inlet direction if using a flux inlet. More...

const std::vector< PostprocessorName > &	getFluxInletPPs () const
	Get the inlet flux postprocessor if using a flux inlet. More...

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

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

virtual const std::string &	name () const

std::string	typeAndName () const

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

void	callMooseError (std::string msg, const bool with_prefix) const

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

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

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

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Public Attributes
const ConsoleStream	_console

Static Public Attributes
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	addPorousMediumSpeedMaterial ()
	Add material to define the local speed in porous medium flows. More...

void	addNonPorousMediumSpeedMaterial ()
	Add material to define the local speed with no porous medium treatment. More...

VariableName	getFlowVariableName (const std::string &default_name) const
	Convenience routine to be able to retrieve the actual variable names from their default names. More...

bool	hasTurbulencePhysics () const
	Whether a turbulence Physics has been coupled in, to know which viscosity to pick on symmetry boundary conditions. More...

const WCNSFVTurbulencePhysics *	getCoupledTurbulencePhysics () const
	Find the turbulence physics. More...

bool	usingNavierStokesFVSyntax () const
	Detects if we are using the new Physics syntax or the old NavierStokesFV action. More...

InputParameters	getAdditionalRMParams () const override
	Parameters to change or add relationship managers. More...

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

bool	addRelationshipManagers (Moose::RelationshipManagerType when_type, const InputParameters &moose_object_pars)

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	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. More...

const bool	_has_flow_equations
	Boolean to keep track of whether the flow equations should be created. More...

const MooseEnum	_compressibility
	Compressibility type, can be compressible, incompressible or weakly-compressible. More...

const bool	_solve_for_dynamic_pressure
	Whether we are solving for the total or dynamic pressure. More...

const bool	_porous_medium_treatment
	Whether to use the porous medium treatment. More...

const MooseFunctorName	_porosity_name
	Name of the porosity functor. More...

MooseFunctorName	_flow_porosity_functor_name
	Name of the porosity functor for the flow equations (if smoothed) More...

const std::vector< std::string >	_velocity_names
	Velocity names. More...

const NonlinearVariableName	_pressure_name
	Pressure name. More...

const NonlinearVariableName	_fluid_temperature_name
	Fluid temperature name. More...

const MooseFunctorName	_density_name
	Name of the density material property. More...

const MooseFunctorName	_density_gravity_name
	Name of the density material property used for gravity and Boussinesq terms. More...

const MooseFunctorName	_dynamic_viscosity_name
	Name of the dynamic viscosity material property. More...

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

const MooseEnum	_momentum_advection_interpolation
	The momentum face interpolation method for being advected. More...

const MooseEnum	_momentum_face_interpolation
	The momentum face interpolation method for stress terms. More...

const WCNSFVTurbulencePhysics *	_turbulence_physics
	Can be set to a coupled turbulence physics. More...

std::vector< std::vector< SubdomainName > >	_friction_blocks
	Subdomains where we want to have volumetric friction. More...

std::vector< std::vector< std::string > >	_friction_types
	The friction correlation types used for each block. More...

std::vector< std::vector< std::string > >	_friction_coeffs
	The coefficients used for each item if friction type. More...

const std::vector< BoundaryName >	_inlet_boundaries
	Boundaries with a flow inlet specified on them. More...

const std::vector< BoundaryName >	_outlet_boundaries
	Boundaries with a flow outlet specified on them. More...

const std::vector< BoundaryName >	_wall_boundaries
	Boundaries which define a wall (slip/noslip/etc.) More...

const std::vector< BoundaryName >	_hydraulic_separators
	Hydraulic separator boundaries. More...

std::map< BoundaryName, MooseEnum >	_momentum_inlet_types
	Momentum inlet boundary types. More...

std::map< BoundaryName, MooseEnum >	_momentum_outlet_types
	Momentum outlet boundary types. More...

std::map< BoundaryName, MooseEnum >	_momentum_wall_types
	Momentum wall boundary types. More...

std::vector< PostprocessorName >	_flux_inlet_pps
	Postprocessors describing the momentum inlet for each boundary. Indexing based on the number of flux boundaries. More...

std::vector< Point >	_flux_inlet_directions
	Direction of each flux inlet. Indexing based on the number of flux boundaries. More...

std::map< BoundaryName, std::vector< MooseFunctorName > >	_momentum_inlet_functors
	Functors describing the momentum inlet for each boundary. More...

std::map< BoundaryName, MooseFunctorName >	_pressure_functors
	Functors describing the outlet pressure on each boundary. More...

std::map< BoundaryName, std::vector< MooseFunctorName > >	_momentum_wall_functors
	Functors describing the momentum for each wall boundary. More...

bool	_define_variables
	Whether to define variables if they do not exist. More...

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

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

Private Member Functions
virtual void	addSolverVariables () override

virtual void	addFVKernels () override

virtual void	addUserObjects () override

void	addPressureCorrectionKernels ()
	Function adding kernels for the incompressible pressure correction equation. More...

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. More...

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. More...

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. More...

virtual void	addRhieChowUserObjects () override
	Function which adds the RhieChow interpolator user objects for weakly and incompressible formulations. More...

virtual void	addFunctorMaterials () override

virtual MooseFunctorName	getLinearFrictionCoefName () const override
	Get the name of the linear friction coefficient. Returns an empty string if no friction. More...

UserObjectName	rhieChowUOName () const override
	Return the name of the Rhie Chow user object. More...

unsigned short	getNumberAlgebraicGhostingLayersNeeded () const override
	Return the number of algebraic ghosting layers needed. More...

Private Attributes
const bool	_non_orthogonal_correction
	Whether to use the correction term for non-orthogonality. More...

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:

friction
other momentum sources and sinks
porous media
transients

Definition at line 24 of file WCNSLinearFVFlowPhysics.h.

Constructor & Destructor Documentation

◆ WCNSLinearFVFlowPhysics()

WCNSLinearFVFlowPhysics::WCNSLinearFVFlowPhysics ( const InputParameters & parameters )

Definition at line 63 of file WCNSLinearFVFlowPhysics.C.

   : WCNSFVFlowPhysicsBase(parameters),
     _non_orthogonal_correction(getParam<bool>("orthogonality_correction"))
 {
   if (_porous_medium_treatment)
     paramError("porous_medium_treatment", "Porous media unsupported");
   if (!_has_flow_equations)
     mooseError("Not supported");
 
   if (_hydraulic_separators.size())
     paramError("hydraulic_separator_sidesets",
                "Flow separators are not supported yet for linearFV!");
   if (getParam<bool>("pin_pressure"))
     paramError("pin_pressure",
                "Pressure pinning is implemented in the executioner for the linear finite volume "
                "segregated solves");
 }

Member Function Documentation

◆ actOnAdditionalTasks()

void WCNSFVFlowPhysicsBase::actOnAdditionalTasks ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 252 of file WCNSFVFlowPhysicsBase.C.

 {
   // Turbulence physics would not be initialized before this task
   if (_current_task == "get_turbulence_physics")
     _turbulence_physics = getCoupledTurbulencePhysics();
 }

◆ addFunctorMaterials()

void WCNSLinearFVFlowPhysics::addFunctorMaterials ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 615 of file WCNSLinearFVFlowPhysics.C.

 {
   if (parameters().isParamValid("gravity"))
   {
     const auto gravity_vector = getParam<RealVectorValue>("gravity");
     const std::vector<std::string> comp_axis({"x", "y", "z"});
     for (const auto d : make_range(dimension()))
       if (gravity_vector(d) != 0)
       {
         // Add rho * g functor for each relevant direction
         // TODO: we could avoid using an AD functor material for non-AD density functor
         auto params = getFactory().getValidParams("ADParsedFunctorMaterial");
         assignBlocks(params, _blocks);
         params.set<std::string>("expression") =
             _density_gravity_name + " * " + std::to_string(gravity_vector(d));
         if (!MooseUtils::parsesToReal(_density_gravity_name))
           params.set<std::vector<std::string>>("functor_names") = {_density_gravity_name};
         params.set<std::string>("property_name") = "rho_g_" + comp_axis[d];
         // We don't output this helper material
         getProblem().addMaterial(
             "ADParsedFunctorMaterial", prefix() + "gravity_helper_" + comp_axis[d], params);
       }
   }
 }

◆ addFVBCs()

void WCNSFVFlowPhysicsBase::addFVBCs ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 260 of file WCNSFVFlowPhysicsBase.C.

 {
   addInletBC();
   addOutletBC();
   addWallsBC();
   addSeparatorBC();
 }

◆ addFVKernels()

void WCNSLinearFVFlowPhysics::addFVKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 162 of file WCNSLinearFVFlowPhysics.C.

 {
   if (!_has_flow_equations)
     return;
 
   // Pressure correction equation: divergence of momentum
   addPressureCorrectionKernels();
 
   // Momentum equation: time derivative
   if (isTransient())
     addMomentumTimeKernels();
 
   // Momentum equation: flux terms
   addMomentumFluxKernels();
 
   // Momentum equation: pressure term
   addMomentumPressureKernels();
 
   // Momentum equation: friction term
   if (_friction_types.size())
     addMomentumFrictionKernels();
 
   // Momentum equation: gravity source term
   addMomentumGravityKernels();
 
   // Momentum equation: boussinesq approximation
   if (getParam<bool>("boussinesq_approximation"))
     addMomentumBoussinesqKernels();
 }

◆ addInitialConditions()

void WCNSFVFlowPhysicsBase::addInitialConditions ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 310 of file WCNSFVFlowPhysicsBase.C.

 {
   if (!_define_variables && parameters().isParamSetByUser("initial_velocity") &&
       parameters().isParamSetByUser("velocity_variable") &&
       getParam<std::vector<FunctionName>>("initial_velocity").size() != 0)
     // TODO: Rework and remove this last statement once the NSFV action is removed
     paramError("initial_velocity",
                "Velocity is defined externally of WCNSFVFlowPhysicsBase, so should the inital "
                "conditions");
   if (!_define_variables && parameters().isParamSetByUser("initial_pressure") &&
       parameters().isParamSetByUser("pressure_variable"))
     paramError("initial_pressure",
                "Pressure is defined externally of WCNSFVFlowPhysicsBase, so should the inital "
                "condition");
 
   // Check dimension
   if (getParam<std::vector<FunctionName>>("initial_velocity").size() != dimension() &&
       getParam<std::vector<FunctionName>>("initial_velocity").size() != 3 &&
       getParam<std::vector<FunctionName>>("initial_velocity").size() != 0)
     // TODO: Rework and remove this last statement once the NSFV action is removed
     paramError("initial_velocity",
                "The number of velocity components in the " + type() + " initial condition is not " +
                    std::to_string(dimension()) + " or 3!");
 
   InputParameters params = getFactory().getValidParams("FunctionIC");
   assignBlocks(params, _blocks);
   auto vvalue = getParam<std::vector<FunctionName>>("initial_velocity");
 
   for (const auto d : make_range(dimension()))
   {
     params.set<VariableName>("variable") = _velocity_names[d];
     params.set<FunctionName>("function") = vvalue[d];
 
     if (shouldCreateIC(_velocity_names[d],
                        _blocks,
                        /*whether IC is a default*/ !isParamSetByUser("initial_velocity"),
                        /*error if already an IC*/ isParamSetByUser("initial_velocity")))
       getProblem().addInitialCondition("FunctionIC", prefix() + _velocity_names[d] + "_ic", params);
   }
 
   if (shouldCreateIC(_pressure_name,
                      _blocks,
                      /*whether IC is a default*/ !isParamSetByUser("initial_pressure"),
                      /*error if already an IC*/ isParamSetByUser("initial_pressure")))
   {
     params.set<VariableName>("variable") = _pressure_name;
     params.set<FunctionName>("function") = getParam<FunctionName>("initial_pressure");
 
     getProblem().addInitialCondition("FunctionIC", prefix() + _pressure_name + "_ic", params);
   }
 }

◆ 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 386 of file WCNSLinearFVFlowPhysics.C.

 {
   // Check the size of the BC parameters
   unsigned int num_velocity_functor_inlets = 0;
   for (const auto & [bdy, momentum_inlet_type] : _momentum_inlet_types)
     if (momentum_inlet_type == "fixed-velocity" || momentum_inlet_type == "fixed-pressure")
       num_velocity_functor_inlets++;
 
   if (num_velocity_functor_inlets != _momentum_inlet_functors.size())
     paramError("momentum_inlet_functors",
                "Size (" + std::to_string(_momentum_inlet_functors.size()) +
                    ") is not the same as the number of entries in the momentum_inlet_types "
                    "subvector for fixed-velocities/pressures functors (size " +
                    std::to_string(num_velocity_functor_inlets) + ")");
 
   unsigned int velocity_pressure_counter = 0;
   for (const auto & [inlet_bdy, momentum_inlet_type] : _momentum_inlet_types)
   {
     if (momentum_inlet_type == "fixed-velocity")
     {
       const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
       InputParameters params = getFactory().getValidParams(bc_type);
       params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
       if (_momentum_inlet_functors.size() < velocity_pressure_counter + 1)
         paramError("momentum_inlet_functors",
                    "More non-flux inlets than inlet functors (" +
                        std::to_string(_momentum_inlet_functors.size()) + ")");
 
       // Check that enough functors have been provided for the dimension of the problem
       const auto momentum_functors = libmesh_map_find(_momentum_inlet_functors, inlet_bdy);
       if (momentum_functors.size() < dimension())
         paramError("momentum_inlet_functors",
                    "Subvector for boundary '" + inlet_bdy + "' (size " +
                        std::to_string(momentum_functors.size()) +
                        ") is not the same size as the number of dimensions of the physics (" +
                        std::to_string(dimension()) + ")");
 
       for (const auto d : make_range(dimension()))
       {
         params.set<LinearVariableName>("variable") = _velocity_names[d];
         params.set<MooseFunctorName>("functor") = momentum_functors[d];
 
         getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + inlet_bdy, params);
       }
       ++velocity_pressure_counter;
 
       // Add the two term BC expansion for pressure if requested
       if (getParam<bool>("pressure_two_term_bc_expansion"))
       {
         const std::string bc_type = "LinearFVExtrapolatedPressureBC";
         InputParameters params = getFactory().getValidParams(bc_type);
         params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
         params.set<LinearVariableName>("variable") = _pressure_name;
         params.set<bool>("use_two_term_expansion") = true;
         getProblem().addLinearFVBC(bc_type,
                                    _pressure_name + "_extrapolation_inlet_" +
                                        Moose::stringify(inlet_bdy),
                                    params);
       }
     }
     else if (momentum_inlet_type == "fixed-pressure")
     {
       const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
       InputParameters params = getFactory().getValidParams(bc_type);
       params.set<LinearVariableName>("variable") = _pressure_name;
       if (_momentum_inlet_functors.size() < velocity_pressure_counter + 1)
         paramError("momentum_inlet_functors",
                    "More non-flux inlets than inlet functors (" +
                        std::to_string(_momentum_inlet_functors.size()) + ")");
 
       params.set<MooseFunctorName>("functor") =
           libmesh_map_find(_momentum_inlet_functors, inlet_bdy)[0];
       params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
 
       getProblem().addLinearFVBC(bc_type, _pressure_name + "_" + inlet_bdy, params);
       ++velocity_pressure_counter;
     }
     else
       mooseError("Unsupported inlet boundary condition type: ", momentum_inlet_type);
   }
 }

◆ addMaterials()

void WCNSFVFlowPhysicsBase::addMaterials ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 269 of file WCNSFVFlowPhysicsBase.C.

 {
   if (hasForchheimerFriction() || _porous_medium_treatment)
     addPorousMediumSpeedMaterial();
   else
     addNonPorousMediumSpeedMaterial();
 }

◆ addMomentumBoussinesqKernels()

void WCNSLinearFVFlowPhysics::addMomentumBoussinesqKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 359 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   if (_compressibility == "weakly-compressible")
     paramError("boussinesq_approximation",
                "We cannot use boussinesq approximation while running in weakly-compressible mode!");
 
   std::string kernel_type = "LinearFVMomentumBoussinesq";
   std::string kernel_name = prefix() + "ins_momentum_boussinesq_";
 
   InputParameters params = getFactory().getValidParams(kernel_type);
   assignBlocks(params, _blocks);
   params.set<VariableName>(NS::T_fluid) = _fluid_temperature_name;
   params.set<MooseFunctorName>(NS::density) = _density_gravity_name;
   params.set<RealVectorValue>("gravity") = getParam<RealVectorValue>("gravity");
   params.set<Real>("ref_temperature") = getParam<Real>("ref_temperature");
   params.set<MooseFunctorName>("alpha_name") = getParam<MooseFunctorName>("thermal_expansion");
 
   for (const auto d : make_range(dimension()))
   {
     params.set<MooseEnum>("momentum_component") = NS::directions[d];
     params.set<LinearVariableName>("variable") = _velocity_names[d];
 
     getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
   }
 }

◆ addMomentumFluxKernels()

void WCNSLinearFVFlowPhysics::addMomentumFluxKernels ( )

private

Definition at line 240 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   const std::string u_names[3] = {"u", "v", "w"};
   std::string kernel_type = "LinearWCNSFVMomentumFlux";
   std::string kernel_name = prefix() + "ins_momentum_flux_";
 
   InputParameters params = getFactory().getValidParams(kernel_type);
   assignBlocks(params, _blocks);
   params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
   params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
   params.set<MooseEnum>("advected_interp_method") = _momentum_advection_interpolation;
   params.set<bool>("use_nonorthogonal_correction") = _non_orthogonal_correction;
   params.set<bool>("use_deviatoric_terms") = getParam<bool>("include_deviatoric_stress");
 
   for (unsigned int i = 0; i < dimension(); ++i)
     params.set<SolverVariableName>(u_names[i]) = _velocity_names[i];
 
   for (const auto d : make_range(dimension()))
   {
     params.set<LinearVariableName>("variable") = _velocity_names[d];
     params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
     getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
   }
 }

◆ addMomentumFrictionKernels()

void WCNSLinearFVFlowPhysics::addMomentumFrictionKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 285 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   unsigned int num_friction_blocks = _friction_blocks.size();
   unsigned int num_used_blocks = num_friction_blocks ? num_friction_blocks : 1;
 
   const std::string kernel_type = "LinearFVMomentumFriction";
   InputParameters params = getFactory().getValidParams(kernel_type);
 
   for (const auto block_i : make_range(num_used_blocks))
   {
     std::string block_name = "";
     if (num_friction_blocks)
     {
       params.set<std::vector<SubdomainName>>("block") = _friction_blocks[block_i];
       block_name = Moose::stringify(_friction_blocks[block_i]);
     }
     else
     {
       assignBlocks(params, _blocks);
       block_name = std::to_string(block_i);
     }
 
     for (const auto d : make_range(dimension()))
     {
       params.set<LinearVariableName>("variable") = _velocity_names[d];
       params.set<MooseEnum>("momentum_component") = NS::directions[d];
       for (unsigned int type_i = 0; type_i < _friction_types[block_i].size(); ++type_i)
       {
         const auto upper_name = MooseUtils::toUpper(_friction_types[block_i][type_i]);
         if (upper_name == "DARCY")
         {
           params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
           params.set<MooseFunctorName>("Darcy_name") = _friction_coeffs[block_i][type_i];
         }
         else
           paramError("momentum_friction_types",
                      "Friction type '",
                      _friction_types[block_i][type_i],
                      "' is not implemented");
       }
 
       getProblem().addLinearFVKernel(kernel_type,
                                      prefix() + "momentum_friction_" + block_name + "_" +
                                          NS::directions[d],
                                      params);
     }
   }
 }

◆ addMomentumGravityKernels()

void WCNSLinearFVFlowPhysics::addMomentumGravityKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 335 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   if (parameters().isParamValid("gravity") && !_solve_for_dynamic_pressure)
   {
     std::string kernel_type = "LinearFVSource";
     std::string kernel_name = prefix() + "ins_momentum_gravity_";
 
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     const auto gravity_vector = getParam<RealVectorValue>("gravity");
     const std::vector<std::string> comp_axis({"x", "y", "z"});
 
     for (const auto d : make_range(dimension()))
       if (gravity_vector(d) != 0)
       {
         params.set<MooseFunctorName>("source_density") = "rho_g_" + comp_axis[d];
         params.set<LinearVariableName>("variable") = _velocity_names[d];
 
         getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
       }
   }
 }

◆ addMomentumPressureKernels()

void WCNSLinearFVFlowPhysics::addMomentumPressureKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 267 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   std::string kernel_type = "LinearFVMomentumPressure";
   std::string kernel_name = prefix() + "ins_momentum_pressure_";
 
   InputParameters params = getFactory().getValidParams(kernel_type);
   assignBlocks(params, _blocks);
   params.set<VariableName>("pressure") = _pressure_name;
 
   for (const auto d : make_range(dimension()))
   {
     params.set<MooseEnum>("momentum_component") = NS::directions[d];
     params.set<LinearVariableName>("variable") = _velocity_names[d];
     getProblem().addLinearFVKernel(kernel_type, kernel_name + NS::directions[d], params);
   }
 }

◆ 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 222 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   std::string kernel_type = "LinearFVTimeDerivative";
   std::string kernel_name = prefix() + "ins_momentum_time";
 
   InputParameters params = getFactory().getValidParams(kernel_type);
   assignBlocks(params, _blocks);
   params.set<MooseFunctorName>("factor") = _density_name;
 
   for (const auto d : make_range(dimension()))
   {
     params.set<LinearVariableName>("variable") = _velocity_names[d];
     if (shouldCreateTimeDerivative(_velocity_names[d], _blocks, /*error if already defined*/ false))
       getProblem().addLinearFVKernel(kernel_type, kernel_name + "_" + NS::directions[d], params);
   }
 }

◆ addNonPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addNonPorousMediumSpeedMaterial ( )

protectedinherited

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

Definition at line 295 of file WCNSFVFlowPhysicsBase.C.

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

 {
   const std::string class_name = "ADVectorMagnitudeFunctorMaterial";
   InputParameters params = getFactory().getValidParams(class_name);
   assignBlocks(params, _blocks);
 
   const std::vector<std::string> param_names{"x_functor", "y_functor", "z_functor"};
   for (unsigned int dim_i = 0; dim_i < dimension(); ++dim_i)
     params.set<MooseFunctorName>(param_names[dim_i]) = _velocity_names[dim_i];
   params.set<MooseFunctorName>("vector_magnitude_name") = NS::speed;
 
   getProblem().addMaterial(class_name, prefix() + "ins_speed_material", params);
 }

◆ addOutletBC()

void WCNSLinearFVFlowPhysics::addOutletBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 469 of file WCNSLinearFVFlowPhysics.C.

 {
   // Check the BCs size
   unsigned int num_pressure_outlets = 0;
   for (const auto & [bdy, momentum_outlet_type] : _momentum_outlet_types)
     if (momentum_outlet_type == "fixed-pressure" ||
         momentum_outlet_type == "fixed-pressure-zero-gradient")
       num_pressure_outlets++;
 
   if (num_pressure_outlets != _pressure_functors.size())
     paramError("pressure_functors",
                "Size (" + std::to_string(_pressure_functors.size()) +
                    ") is not the same as the number of pressure outlet boundaries in "
                    "'fixed-pressure/fixed-pressure-zero-gradient' (size " +
                    std::to_string(num_pressure_outlets) + ")");
 
   const std::string u_names[3] = {"u", "v", "w"};
   for (const auto & [outlet_bdy, momentum_outlet_type] : _momentum_outlet_types)
   {
     // Zero tangeantial gradient condition on velocity
     if (momentum_outlet_type == "zero-gradient" || momentum_outlet_type == "fixed-pressure" ||
         momentum_outlet_type == "fixed-pressure-zero-gradient")
     {
       const std::string bc_type = "LinearFVAdvectionDiffusionOutflowBC";
       InputParameters params = getFactory().getValidParams(bc_type);
       params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
       params.set<bool>("use_two_term_expansion") = getParam<bool>("momentum_two_term_bc_expansion");
 
       for (const auto d : make_range(dimension()))
       {
         params.set<LinearVariableName>("variable") = _velocity_names[d];
         getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + outlet_bdy, params);
       }
     }
 
     // Fixed pressure condition, coming in the pressure correction equation
     if (momentum_outlet_type == "fixed-pressure" ||
         momentum_outlet_type == "fixed-pressure-zero-gradient")
     {
       const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
       InputParameters params = getFactory().getValidParams(bc_type);
       params.set<LinearVariableName>("variable") = _pressure_name;
       params.set<MooseFunctorName>("functor") = libmesh_map_find(_pressure_functors, outlet_bdy);
       params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
 
       getProblem().addLinearFVBC(bc_type, _pressure_name + "_" + outlet_bdy, params);
     }
   }
 }

◆ addPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addPorousMediumSpeedMaterial ( )

protectedinherited

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

Definition at line 278 of file WCNSFVFlowPhysicsBase.C.

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

 {
   InputParameters params = getFactory().getValidParams("PINSFVSpeedFunctorMaterial");
   assignBlocks(params, _blocks);
 
   for (unsigned int dim_i = 0; dim_i < dimension(); ++dim_i)
     params.set<MooseFunctorName>(NS::superficial_velocity_vector[dim_i]) = _velocity_names[dim_i];
   if (_porous_medium_treatment)
     params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
   else
     params.set<MooseFunctorName>(NS::porosity) = "1";
   params.set<bool>("define_interstitial_velocity_components") = _porous_medium_treatment;
 
   getProblem().addMaterial("PINSFVSpeedFunctorMaterial", prefix() + "pins_speed_material", params);
 }

◆ addPostprocessors()

void WCNSFVFlowPhysicsBase::addPostprocessors ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 370 of file WCNSFVFlowPhysicsBase.C.

 {
   const auto momentum_inlet_types = getParam<MultiMooseEnum>("momentum_inlet_types");
 
   for (unsigned int bc_ind = 0; bc_ind < momentum_inlet_types.size(); ++bc_ind)
     if (momentum_inlet_types[bc_ind] == "flux-mass" ||
         momentum_inlet_types[bc_ind] == "flux-velocity")
     {
       const std::string pp_type = "AreaPostprocessor";
       InputParameters params = getFactory().getValidParams(pp_type);
       params.set<std::vector<BoundaryName>>("boundary") = {_inlet_boundaries[bc_ind]};
       params.set<ExecFlagEnum>("execute_on") = EXEC_INITIAL;
 
       const auto name_pp = "area_pp_" + _inlet_boundaries[bc_ind];
       if (!getProblem().hasUserObject(name_pp))
         getProblem().addPostprocessor(pp_type, name_pp, params);
     }
 }

◆ addPressureCorrectionKernels()

void WCNSLinearFVFlowPhysics::addPressureCorrectionKernels ( )

private

Function adding kernels for the incompressible pressure correction equation.

Definition at line 193 of file WCNSLinearFVFlowPhysics.C.

Referenced by addFVKernels().

 {
   {
     std::string kernel_type = "LinearFVAnisotropicDiffusion";
     std::string kernel_name = prefix() + "p_diffusion";
 
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<LinearVariableName>("variable") = _pressure_name;
     params.set<MooseFunctorName>("diffusion_tensor") = "Ainv";
     params.set<bool>("use_nonorthogonal_correction") = _non_orthogonal_correction;
 
     getProblem().addLinearFVKernel(kernel_type, kernel_name, params);
   }
   {
     std::string kernel_type = "LinearFVDivergence";
     std::string kernel_name = prefix() + "HbyA_divergence";
 
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<LinearVariableName>("variable") = _pressure_name;
     params.set<MooseFunctorName>("face_flux") = "HbyA";
     params.set<bool>("force_boundary_execution") = true;
 
     getProblem().addLinearFVKernel(kernel_type, kernel_name, params);
   }
 }

◆ addRhieChowUserObjects()

void WCNSLinearFVFlowPhysics::addRhieChowUserObjects ( )

overrideprivatevirtual

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

Implements WCNSFVFlowPhysicsBase.

Definition at line 567 of file WCNSLinearFVFlowPhysics.C.

Referenced by addUserObjects().

 {
   mooseAssert(dimension(), "0-dimension not supported");
 
   // First make sure that we only add this object once
   // Potential cases:
   // - there is a flow physics, and an advection one (UO should be added by one)
   // - there is only an advection physics (UO should be created)
   // - there are two advection physics on different blocks with set velocities (first one picks)
   // Counting RC UOs defined on the same blocks seems to be the most fool proof option
   std::vector<UserObject *> objs;
   getProblem()
       .theWarehouse()
       .query()
       .condition<AttribSystem>("UserObject")
       .condition<AttribThread>(0)
       .queryInto(objs);
   unsigned int num_rc_uo = 0;
   for (const auto & obj : objs)
     if (dynamic_cast<RhieChowMassFlux *>(obj))
     {
       const auto rc_obj = dynamic_cast<RhieChowMassFlux *>(obj);
       if (rc_obj->blocks() == _blocks)
         num_rc_uo++;
       // one of the RC user object is defined everywhere
       else if (rc_obj->blocks().size() == 0 || _blocks.size() == 0)
         num_rc_uo++;
     }
 
   if (num_rc_uo)
     return;
 
   const std::string u_names[3] = {"u", "v", "w"};
   const auto object_type = "RhieChowMassFlux";
 
   auto params = getFactory().getValidParams(object_type);
   assignBlocks(params, _blocks);
   for (unsigned int d = 0; d < dimension(); ++d)
     params.set<VariableName>(u_names[d]) = _velocity_names[d];
 
   params.set<VariableName>("pressure") = _pressure_name;
   params.set<std::string>("p_diffusion_kernel") = prefix() + "p_diffusion";
   params.set<MooseFunctorName>(NS::density) = _density_name;
 
   getProblem().addUserObject(object_type, rhieChowUOName(), params);
 }

◆ addSeparatorBC()

virtual void WCNSLinearFVFlowPhysics::addSeparatorBC ( )

inlineoverrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 57 of file WCNSLinearFVFlowPhysics.h.

57 {}

◆ addSolverVariables()

void WCNSLinearFVFlowPhysics::addSolverVariables ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 102 of file WCNSLinearFVFlowPhysics.C.

 {
   if (!_has_flow_equations)
     return;
 
   for (const auto d : make_range(dimension()))
     saveSolverVariableName(_velocity_names[d]);
   saveSolverVariableName(_pressure_name);
 
   const std::vector<std::string> v_short = {"u", "v", "w"};
 
   // Check number of variables
   if (_velocity_names.size() != dimension() && _velocity_names.size() != 3)
     paramError("velocity_variable",
                "The number of velocity variable names supplied to the NSFVAction is not " +
                    Moose::stringify(dimension()) + " (mesh dimension)" +
                    ((dimension() == 3) ? "" : " or 3!") + "\nVelocity variables " +
                    Moose::stringify(_velocity_names));
 
   // Velocities
   for (const auto d : make_range(dimension()))
   {
     if (!shouldCreateVariable(_velocity_names[d], _blocks, /*error if aux*/ true))
       reportPotentiallyMissedParameters({"system_names"}, "MooseLinearVariableFVReal");
     else if (_define_variables)
     {
       std::string variable_type = "MooseLinearVariableFVReal";
 
       auto params = getFactory().getValidParams(variable_type);
       assignBlocks(params, _blocks);
       params.set<SolverSystemName>("solver_sys") = getSolverSystem(_velocity_names[d]);
 
       getProblem().addVariable(variable_type, _velocity_names[d], params);
     }
     else
       paramError("velocity_variable",
                  "Variable (" + _velocity_names[d] +
                      ") supplied to the WCNSLinearFVFlowPhysics does not exist!");
   }
 
   // Pressure
   if (!shouldCreateVariable(_pressure_name, _blocks, /*error if aux*/ true))
     reportPotentiallyMissedParameters({"system_names"}, "MooseLinearVariableFVReal");
   else if (_define_variables)
   {
     const auto pressure_type = "MooseLinearVariableFVReal";
 
     auto params = getFactory().getValidParams(pressure_type);
     assignBlocks(params, _blocks);
     params.set<SolverSystemName>("solver_sys") = getSolverSystem(_pressure_name);
 
     getProblem().addVariable(pressure_type, _pressure_name, params);
   }
   else
     paramError("pressure_variable",
                "Variable (" + _pressure_name +
                    ") supplied to the WCNSLinearFVFlowPhysics does not exist!");
 }

◆ addUserObjects()

void WCNSLinearFVFlowPhysics::addUserObjects ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 560 of file WCNSLinearFVFlowPhysics.C.

 {
   // Rhie Chow user object for interpolation velocities
   addRhieChowUserObjects();
 }

◆ addWallsBC()

void WCNSLinearFVFlowPhysics::addWallsBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 520 of file WCNSLinearFVFlowPhysics.C.

 {
   const std::string u_names[3] = {"u", "v", "w"};
 
   for (const auto & [boundary_name, momentum_wall_type] : _momentum_wall_types)
   {
     if (momentum_wall_type == "noslip")
     {
       const std::string bc_type = "LinearFVAdvectionDiffusionFunctorDirichletBC";
       InputParameters params = getFactory().getValidParams(bc_type);
       params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
 
       for (const auto d : make_range(dimension()))
       {
         params.set<LinearVariableName>("variable") = _velocity_names[d];
         if (_momentum_wall_functors.count(boundary_name) == 0)
           params.set<MooseFunctorName>("functor") = "0";
         else
           params.set<MooseFunctorName>("functor") = _momentum_wall_functors[boundary_name][d];
 
         getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
       }
     }
     else
       mooseError("Unsupported wall boundary condition type: " + std::string(momentum_wall_type));
   }
 
   if (getParam<bool>("pressure_two_term_bc_expansion"))
   {
     const std::string bc_type = "LinearFVExtrapolatedPressureBC";
     InputParameters params = getFactory().getValidParams(bc_type);
     params.set<std::vector<BoundaryName>>("boundary") = _wall_boundaries;
     params.set<LinearVariableName>("variable") = _pressure_name;
     params.set<bool>("use_two_term_expansion") = true;
     getProblem().addLinearFVBC(
         bc_type, _pressure_name + "_extrapolation_" + Moose::stringify(_wall_boundaries), params);
   }
 }

◆ compressibility()

const MooseEnum& WCNSFVFlowPhysicsBase::compressibility ( ) const

inlineinherited

Return the compressibility of the flow equations selected.

Definition at line 47 of file WCNSFVFlowPhysicsBase.h.

47 { return _compressibility; }

WCNSFVFlowPhysicsBase::_compressibility

const MooseEnum _compressibility

Compressibility type, can be compressible, incompressible or weakly-compressible. ...

Definition: WCNSFVFlowPhysicsBase.h:149

◆ densityName()

const MooseFunctorName& WCNSFVFlowPhysicsBase::densityName ( ) const

inlineinherited

Return the name of the density functor.

Definition at line 53 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVTurbulencePhysics::addAuxiliaryKernels(), WCNSFVTurbulencePhysics::addFVBCs(), WCNSFVTurbulencePhysics::addInitialConditions(), WCNSFVTurbulencePhysics::addKEpsilonAdvection(), WCNSFVTurbulencePhysics::addKEpsilonSink(), WCNSFVTurbulencePhysics::addMaterials(), PNSFVSolidHeatTransferPhysics::checkFluidAndSolidHeatTransferPhysicsParameters(), WCNSFVTwoPhaseMixturePhysics::WCNSFVTwoPhaseMixturePhysics(), and WCNSLinearFVTwoPhaseMixturePhysics::WCNSLinearFVTwoPhaseMixturePhysics().

53 { return _density_name; }

WCNSFVFlowPhysicsBase::_density_name

const MooseFunctorName _density_name

Name of the density material property.

Definition: WCNSFVFlowPhysicsBase.h:168

◆ dynamicViscosityName()

const MooseFunctorName& WCNSFVFlowPhysicsBase::dynamicViscosityName ( ) const

inlineinherited

Return the name of the dynamic viscosity functor.

Definition at line 55 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVTurbulencePhysics::addAuxiliaryKernels(), WCNSFVTurbulencePhysics::addFlowTurbulenceKernels(), WCNSFVTurbulencePhysics::addFVBCs(), WCNSFVTurbulencePhysics::addKEpsilonDiffusion(), WCNSFVTurbulencePhysics::addKEpsilonSink(), and WCNSFVTurbulencePhysics::addMaterials().

55 { return _dynamic_viscosity_name; }

WCNSFVFlowPhysicsBase::_dynamic_viscosity_name

const MooseFunctorName _dynamic_viscosity_name

Name of the dynamic viscosity material property.

Definition: WCNSFVFlowPhysicsBase.h:172

◆ getAdditionalRMParams()

InputParameters NavierStokesPhysicsBase::getAdditionalRMParams ( ) const

overrideprotectedvirtualinherited

Parameters to change or add relationship managers.

Reimplemented from PhysicsBase.

Definition at line 42 of file NavierStokesPhysicsBase.C.

 {
   unsigned short necessary_layers = getParam<unsigned short>("ghost_layers");
   necessary_layers = std::max(necessary_layers, getNumberAlgebraicGhostingLayersNeeded());
 
   // Just an object that has a ghost_layers parameter
   const std::string kernel_type = "INSFVMixingLengthReynoldsStress";
   InputParameters params = getFactory().getValidParams(kernel_type);
   params.template set<unsigned short>("ghost_layers") = necessary_layers;
 
   return params;
 }

◆ getCoupledTurbulencePhysics()

const WCNSFVTurbulencePhysics * WCNSFVFlowPhysicsBase::getCoupledTurbulencePhysics ( ) const

protectedinherited

Find the turbulence physics.

Definition at line 416 of file WCNSFVFlowPhysicsBase.C.

Referenced by WCNSFVFlowPhysicsBase::actOnAdditionalTasks().

 {
   // User passed it, just use that
   if (isParamValid("coupled_turbulence_physics"))
     return getCoupledPhysics<WCNSFVTurbulencePhysics>(
         getParam<PhysicsName>("coupled_flow_physics"));
   // Look for any physics of the right type, and check the block restriction
   else
   {
     const auto all_turbulence_physics = getCoupledPhysics<const WCNSFVTurbulencePhysics>(true);
     for (const auto physics : all_turbulence_physics)
       if (checkBlockRestrictionIdentical(
               physics->name(), physics->blocks(), /*error_if_not_identical=*/false))
         return physics;
   }
   // Did not find one
   return nullptr;
 }

◆ 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 390 of file WCNSFVFlowPhysicsBase.C.

 {
   if (short_name == NS::pressure)
     return getPressureName();
   else if (short_name == NS::velocity_x && dimension() > 0)
     return getVelocityNames()[0];
   else if (short_name == NS::velocity_y && dimension() > 1)
     return getVelocityNames()[1];
   else if (short_name == NS::velocity_z && dimension() > 2)
     return getVelocityNames()[2];
   else if (short_name == NS::temperature)
     return getFluidTemperatureName();
   else
     mooseError("Short Variable name '", short_name, "' not recognized.");
 }

◆ getFluidTemperatureName()

const NonlinearVariableName& WCNSFVFlowPhysicsBase::getFluidTemperatureName ( ) const

inlineinherited

Definition at line 42 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVTurbulencePhysics::addFluidEnergyTurbulenceKernels(), and WCNSFVFlowPhysicsBase::getFlowVariableName().

42 { return _fluid_temperature_name; }

WCNSFVFlowPhysicsBase::_fluid_temperature_name

const NonlinearVariableName _fluid_temperature_name

Fluid temperature name.

Definition: WCNSFVFlowPhysicsBase.h:165

◆ getFluxInletDirections()

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

inlineinherited

Get the inlet direction if using a flux inlet.

Definition at line 77 of file WCNSFVFlowPhysicsBase.h.

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

77 { return _flux_inlet_directions; }

WCNSFVFlowPhysicsBase::_flux_inlet_directions

std::vector< Point > _flux_inlet_directions

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

Definition: WCNSFVFlowPhysicsBase.h:210

◆ getFluxInletPPs()

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

inlineinherited

Get the inlet flux postprocessor if using a flux inlet.

Definition at line 79 of file WCNSFVFlowPhysicsBase.h.

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

79 { return _flux_inlet_pps; }

WCNSFVFlowPhysicsBase::_flux_inlet_pps

std::vector< PostprocessorName > _flux_inlet_pps

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

Definition: WCNSFVFlowPhysicsBase.h:208

◆ getHydraulicSeparators()

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

inlineinherited

Get the hydraulic separator boundaries.

Definition at line 75 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergySeparatorBC().

75 { return _hydraulic_separators; }

WCNSFVFlowPhysicsBase::_hydraulic_separators

const std::vector< BoundaryName > _hydraulic_separators

Hydraulic separator boundaries.

Definition: WCNSFVFlowPhysicsBase.h:198

◆ getInletBoundaries()

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

inlineinherited

Get the inlet boundaries.

Definition at line 69 of file WCNSFVFlowPhysicsBase.h.

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

69 { return _inlet_boundaries; }

WCNSFVFlowPhysicsBase::_inlet_boundaries

const std::vector< BoundaryName > _inlet_boundaries

Boundaries with a flow inlet specified on them.

Definition: WCNSFVFlowPhysicsBase.h:192

◆ 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 64 of file WCNSLinearFVFlowPhysics.h.

   {
     mooseError("Not implemented");
   }

◆ getMomentumAdvectionFaceInterpolationMethod()

const MooseEnum& WCNSFVFlowPhysicsBase::getMomentumAdvectionFaceInterpolationMethod ( ) const

inlineinherited

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

Definition at line 59 of file WCNSFVFlowPhysicsBase.h.

   {
     return _momentum_advection_interpolation;
   }

◆ getMomentumFaceInterpolationMethod()

const MooseEnum& WCNSFVFlowPhysicsBase::getMomentumFaceInterpolationMethod ( ) const

inlineinherited

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

Definition at line 64 of file WCNSFVFlowPhysicsBase.h.

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

   {
     return _momentum_face_interpolation;
   }

◆ getNumberAlgebraicGhostingLayersNeeded()

unsigned short WCNSLinearFVFlowPhysics::getNumberAlgebraicGhostingLayersNeeded ( ) const

overrideprivatevirtual

Return the number of algebraic ghosting layers needed.

Reimplemented from WCNSFVFlowPhysicsBase.

Definition at line 648 of file WCNSLinearFVFlowPhysics.C.

 {
   return 1;
 }

◆ getOutletBoundaries()

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

inlineinherited

Get the outlet boundaries.

Definition at line 71 of file WCNSFVFlowPhysicsBase.h.

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

71 { return _outlet_boundaries; }

WCNSFVFlowPhysicsBase::_outlet_boundaries

const std::vector< BoundaryName > _outlet_boundaries

Boundaries with a flow outlet specified on them.

Definition: WCNSFVFlowPhysicsBase.h:194

◆ getPorosityFunctorName()

MooseFunctorName WCNSFVFlowPhysicsBase::getPorosityFunctorName ( const bool smoothed ) const

inherited

Definition at line 407 of file WCNSFVFlowPhysicsBase.C.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyHeatConductionKernels(), WCNSFVFluidHeatTransferPhysics::addEnergyTimeKernels(), WCNSFVScalarTransportPhysics::addScalarAdvectionKernels(), and WCNSFVCoupledAdvectionPhysicsHelper::getPorosityFunctorName().

 {
   if (smoothed)
     return _flow_porosity_functor_name;
   else
     return _porosity_name;
 }

◆ getPressureName()

const NonlinearVariableName& WCNSFVFlowPhysicsBase::getPressureName ( ) const

inlineinherited

Definition at line 41 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFluidHeatTransferPhysicsBase::addMaterials(), and WCNSFVFlowPhysicsBase::getFlowVariableName().

41 { return _pressure_name; }

WCNSFVFlowPhysicsBase::_pressure_name

const NonlinearVariableName _pressure_name

Pressure name.

Definition: WCNSFVFlowPhysicsBase.h:163

◆ getVelocityFaceInterpolationMethod()

const MooseEnum& WCNSFVFlowPhysicsBase::getVelocityFaceInterpolationMethod ( ) const

inlineinherited

Get the face interpolation method for velocity.

Definition at line 57 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm().

57 { return _velocity_interpolation; }

WCNSFVFlowPhysicsBase::_velocity_interpolation

const MooseEnum _velocity_interpolation

The velocity face interpolation method for advecting other quantities.

Definition: WCNSFVFlowPhysicsBase.h:175

◆ getVelocityNames()

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

inlineinherited

To interface with other Physics.

Definition at line 40 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm(), WCNSLinearFVTwoPhaseMixturePhysics::addMaterials(), WCNSFVTwoPhaseMixturePhysics::addMaterials(), WCNSLinearFVTwoPhaseMixturePhysics::addPhaseDriftFluxTerm(), WCNSFVTwoPhaseMixturePhysics::addPhaseDriftFluxTerm(), and WCNSFVFlowPhysicsBase::getFlowVariableName().

40 { return _velocity_names; }

WCNSFVFlowPhysicsBase::_velocity_names

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

Velocity names.

Definition: WCNSFVFlowPhysicsBase.h:161

◆ getWallBoundaries()

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

inlineinherited

Get the wall boundaries.

Definition at line 73 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyWallBC(), and WCNSLinearFVFluidHeatTransferPhysics::addEnergyWallBC().

73 { return _wall_boundaries; }

WCNSFVFlowPhysicsBase::_wall_boundaries

const std::vector< BoundaryName > _wall_boundaries

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

Definition: WCNSFVFlowPhysicsBase.h:196

◆ gravityVector()

RealVectorValue WCNSFVFlowPhysicsBase::gravityVector ( ) const

inlineinherited

Return the gravity vector.

Definition at line 51 of file WCNSFVFlowPhysicsBase.h.

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

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

◆ hasFlowEquations()

bool WCNSFVFlowPhysicsBase::hasFlowEquations ( ) const

inlineinherited

Whether the physics is actually creating the flow equations.

Definition at line 37 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSLinearFVTwoPhaseMixturePhysics::addFVKernels(), WCNSFVTwoPhaseMixturePhysics::addFVKernels(), and WCNSFVTurbulencePhysics::retrieveCoupledPhysics().

37 { return _has_flow_equations; }

WCNSFVFlowPhysicsBase::_has_flow_equations

const bool _has_flow_equations

Boolean to keep track of whether the flow equations should be created.

Definition: WCNSFVFlowPhysicsBase.h:146

◆ hasForchheimerFriction()

virtual bool WCNSLinearFVFlowPhysics::hasForchheimerFriction ( ) const

inlineoverrideprivatevirtual

Return whether a Forchheimer friction model is in use.

Implements WCNSFVFlowPhysicsBase.

Definition at line 59 of file WCNSLinearFVFlowPhysics.h.

59 { 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 131 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addWallsBC().

   {
     if (_turbulence_physics)
       return _turbulence_physics->hasTurbulenceModel();
     else
       return false;
   }

◆ initializePhysicsAdditional()

void WCNSLinearFVFlowPhysics::initializePhysicsAdditional ( )

overrideprotectedvirtual

Reimplemented from WCNSFVFlowPhysicsBase.

Definition at line 82 of file WCNSLinearFVFlowPhysics.C.

 {
   WCNSFVFlowPhysicsBase::initializePhysicsAdditional();
   // TODO Add support for multi-system by either:
   // - creating the problem in the Physics or,
   // - checking that the right systems are being created
   getProblem().needSolutionState(2, Moose::SolutionIterationType::Nonlinear);
   // TODO Ban all other nonlinear Physics for now
 
   // Fix the default system names if using a different dimension
   if (!isParamSetByUser("system_name"))
   {
     if (dimension() == 1)
       _system_names = {"u_system", "pressure_system"};
     else if (dimension() == 2)
       _system_names = {"u_system", "v_system", "pressure_system"};
   }
 }

◆ porousMediumTreatment()

bool WCNSFVFlowPhysicsBase::porousMediumTreatment ( ) const

inlineinherited

Return whether a porous medium treatment is applied.

Definition at line 49 of file WCNSFVFlowPhysicsBase.h.

49 { return _porous_medium_treatment; }

WCNSFVFlowPhysicsBase::_porous_medium_treatment

const bool _porous_medium_treatment

Whether to use the porous medium treatment.

Definition: WCNSFVFlowPhysicsBase.h:154

◆ rhieChowUOName()

UserObjectName WCNSLinearFVFlowPhysics::rhieChowUOName ( ) const

overrideprivatevirtual

Return the name of the Rhie Chow user object.

Implements WCNSFVFlowPhysicsBase.

Definition at line 641 of file WCNSLinearFVFlowPhysics.C.

Referenced by addMomentumFluxKernels(), and addRhieChowUserObjects().

 {
   mooseAssert(!_porous_medium_treatment, "Not implemented");
   return "ins_rhie_chow_interpolator";
 }

◆ usingNavierStokesFVSyntax()

bool NavierStokesPhysicsBase::usingNavierStokesFVSyntax ( ) const

inlineprotectedinherited

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

Definition at line 32 of file NavierStokesPhysicsBase.h.

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

◆ validParams()

InputParameters WCNSLinearFVFlowPhysics::validParams ( )

static

Definition at line 25 of file WCNSLinearFVFlowPhysics.C.

 {
   InputParameters params = WCNSFVFlowPhysicsBase::validParams();
   params.addClassDescription(
       "Define the Navier Stokes weakly-compressible equations with the linear "
       "solver implementation of the SIMPLE scheme");
 
   params.addParam<bool>(
       "orthogonality_correction", false, "Whether to use orthogonality correction");
   params.set<unsigned short>("ghost_layers") = 1;
 
   // This will be adapted based on the dimension
   params.set<std::vector<SolverSystemName>>("system_names") = {
       "u_system", "v_system", "w_system", "pressure_system"};
 
   // Implemented in the executioner
   params.suppressParameter<MooseEnum>("pinned_pressure_type");
   params.suppressParameter<Point>("pinned_pressure_point");
   params.suppressParameter<PostprocessorName>("pinned_pressure_value");
 
   // Not supported
   params.suppressParameter<bool>("add_flow_equations");
   params.set<bool>("porous_medium_treatment") = false;
   params.suppressParameter<bool>("porous_medium_treatment");
   params.set<MooseFunctorName>("porosity") = "1";
   params.suppressParameter<MooseFunctorName>("porosity");
   params.suppressParameter<MooseEnum>("mu_interp_method");
   // Not needed
   params.suppressParameter<bool>("add_flow_equations");
   params.suppressParameter<MooseEnum>("preconditioning");
 
   // No other options so far
   params.set<MooseEnum>("velocity_interpolation") = "rc";
   params.suppressParameter<MooseEnum>("velocity_interpolation");
 
   return params;
 }

Member Data Documentation

◆ _compressibility

const MooseEnum WCNSFVFlowPhysicsBase::_compressibility

protectedinherited

Compressibility type, can be compressible, incompressible or weakly-compressible.

Definition at line 149 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addFVKernels(), addMomentumBoussinesqKernels(), WCNSFVFlowPhysics::addMomentumBoussinesqKernels(), WCNSFVFlowPhysics::addMomentumTimeKernels(), WCNSFVFlowPhysicsBase::compressibility(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _define_variables

bool NavierStokesPhysicsBase::_define_variables

protectedinherited

Whether to define variables if they do not exist.

Definition at line 44 of file NavierStokesPhysicsBase.h.

Referenced by WCNSLinearFVFluidHeatTransferPhysics::addAuxiliaryVariables(), WCNSFVTurbulencePhysics::addAuxiliaryVariables(), WCNSFVTurbulencePhysics::addInitialConditions(), WCNSFVFluidHeatTransferPhysicsBase::addInitialConditions(), WCNSFVScalarTransportPhysicsBase::addInitialConditions(), WCNSFVFlowPhysicsBase::addInitialConditions(), WCNSFVFluidHeatTransferPhysics::addSolverVariables(), WCNSLinearFVFluidHeatTransferPhysics::addSolverVariables(), WCNSFVFlowPhysics::addSolverVariables(), addSolverVariables(), WCNSFVTurbulencePhysics::addSolverVariables(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _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 170 of file WCNSFVFlowPhysicsBase.h.

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

◆ _density_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_name

protectedinherited

Name of the density material property.

Definition at line 168 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addInletBC(), WCNSFVFlowPhysics::addMassKernels(), WCNSFVFlowPhysics::addMassTimeKernels(), WCNSFVFlowPhysics::addMomentumAdvectionKernels(), WCNSFVFlowPhysics::addMomentumFrictionKernels(), addMomentumTimeKernels(), WCNSFVFlowPhysics::addMomentumTimeKernels(), WCNSFVFlowPhysics::addOutletBC(), addRhieChowUserObjects(), WCNSFVFlowPhysics::addSolverVariables(), WCNSFVFlowPhysics::addWallsBC(), and WCNSFVFlowPhysicsBase::densityName().

◆ _dynamic_viscosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_dynamic_viscosity_name

protectedinherited

Name of the dynamic viscosity material property.

Definition at line 172 of file WCNSFVFlowPhysicsBase.h.

Referenced by addMomentumFluxKernels(), addMomentumFrictionKernels(), WCNSFVFlowPhysics::addMomentumFrictionKernels(), WCNSFVFlowPhysics::addMomentumViscousDissipationKernels(), WCNSFVFlowPhysics::addWallsBC(), and WCNSFVFlowPhysicsBase::dynamicViscosityName().

◆ _flow_porosity_functor_name

MooseFunctorName WCNSFVFlowPhysicsBase::_flow_porosity_functor_name

protectedinherited

Name of the porosity functor for the flow equations (if smoothed)

Definition at line 158 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addMassTimeKernels(), WCNSFVFlowPhysics::addMomentumAdvectionKernels(), WCNSFVFlowPhysics::addMomentumBoussinesqKernels(), WCNSFVFlowPhysics::addMomentumGravityKernels(), WCNSFVFlowPhysics::addMomentumPressureKernels(), WCNSFVFlowPhysics::addMomentumViscousDissipationKernels(), WCNSFVFlowPhysics::addOutletBC(), WCNSFVFlowPhysicsBase::addPorousMediumSpeedMaterial(), WCNSFVFlowPhysicsBase::getPorosityFunctorName(), and WCNSFVFlowPhysics::WCNSFVFlowPhysics().

◆ _fluid_temperature_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_fluid_temperature_name

protectedinherited

Fluid temperature name.

Definition at line 165 of file WCNSFVFlowPhysicsBase.h.

Referenced by addMomentumBoussinesqKernels(), WCNSFVFlowPhysics::addMomentumBoussinesqKernels(), and WCNSFVFlowPhysicsBase::getFluidTemperatureName().

◆ _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 210 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 208 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 185 of file WCNSFVFlowPhysicsBase.h.

Referenced by addMomentumFrictionKernels(), WCNSFVFlowPhysics::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 189 of file WCNSFVFlowPhysicsBase.h.

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

◆ _friction_types

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

protectedinherited

The friction correlation types used for each block.

Definition at line 187 of file WCNSFVFlowPhysicsBase.h.

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

◆ _has_flow_equations

const bool WCNSFVFlowPhysicsBase::_has_flow_equations

protectedinherited

Boolean to keep track of whether the flow equations should be created.

Definition at line 146 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addCorrectors(), WCNSFVFlowPhysics::addFVKernels(), addFVKernels(), WCNSFVFlowPhysics::addRhieChowUserObjects(), WCNSFVFlowPhysics::addSolverVariables(), addSolverVariables(), WCNSFVFlowPhysicsBase::hasFlowEquations(), and WCNSLinearFVFlowPhysics().

◆ _hydraulic_separators

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

protectedinherited

Hydraulic separator boundaries.

Definition at line 198 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addSeparatorBC(), WCNSFVFlowPhysicsBase::getHydraulicSeparators(), and WCNSLinearFVFlowPhysics().

◆ _inlet_boundaries

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

protectedinherited

Boundaries with a flow inlet specified on them.

Definition at line 192 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_advection_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_advection_interpolation

protectedinherited

The momentum face interpolation method for being advected.

Definition at line 177 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_face_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_face_interpolation

protectedinherited

The momentum face interpolation method for stress terms.

Definition at line 179 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 213 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_inlet_types

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

protectedinherited

Momentum inlet boundary types.

Definition at line 201 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_outlet_types

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

protectedinherited

Momentum outlet boundary types.

Definition at line 203 of file WCNSFVFlowPhysicsBase.h.

Referenced by addOutletBC(), WCNSFVFlowPhysics::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 217 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_wall_types

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

protectedinherited

Momentum wall boundary types.

Definition at line 205 of file WCNSFVFlowPhysicsBase.h.

Referenced by addWallsBC(), WCNSFVFlowPhysics::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 74 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 194 of file WCNSFVFlowPhysicsBase.h.

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

◆ _porosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_porosity_name

protectedinherited

Name of the porosity functor.

Definition at line 156 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addInletBC(), WCNSFVFlowPhysics::addRhieChowUserObjects(), WCNSFVFlowPhysics::addSolverVariables(), WCNSFVFlowPhysicsBase::getPorosityFunctorName(), and WCNSFVFlowPhysics::WCNSFVFlowPhysics().

◆ _porous_medium_treatment

const bool WCNSFVFlowPhysicsBase::_porous_medium_treatment

protectedinherited

Whether to use the porous medium treatment.

Definition at line 154 of file WCNSFVFlowPhysicsBase.h.

◆ _pressure_functors

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

protectedinherited

Functors describing the outlet pressure on each boundary.

Definition at line 215 of file WCNSFVFlowPhysicsBase.h.

Referenced by addOutletBC(), WCNSFVFlowPhysics::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 143 of file WCNSFVFlowPhysicsBase.h.

◆ _solve_for_dynamic_pressure

const bool WCNSFVFlowPhysicsBase::_solve_for_dynamic_pressure

protectedinherited

Whether we are solving for the total or dynamic pressure.

Definition at line 151 of file WCNSFVFlowPhysicsBase.h.

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

◆ _turbulence_physics

const WCNSFVTurbulencePhysics* WCNSFVFlowPhysicsBase::_turbulence_physics

protectedinherited

Can be set to a coupled turbulence physics.

Definition at line 182 of file WCNSFVFlowPhysicsBase.h.

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

◆ _velocity_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_velocity_interpolation

protectedinherited

The velocity face interpolation method for advecting other quantities.

Definition at line 175 of file WCNSFVFlowPhysicsBase.h.

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

◆ _velocity_names

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

protectedinherited

Velocity names.

Definition at line 161 of file WCNSFVFlowPhysicsBase.h.

◆ _wall_boundaries

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

protectedinherited

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

Definition at line 196 of file WCNSFVFlowPhysicsBase.h.

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

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

navier_stokes/include/physics/WCNSLinearFVFlowPhysics.h
navier_stokes/src/physics/WCNSLinearFVFlowPhysics.C

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

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ WCNSLinearFVFlowPhysics()

Member Function Documentation

◆ actOnAdditionalTasks()

◆ addFunctorMaterials()

◆ addFVBCs()

◆ addFVKernels()

◆ addInitialConditions()

◆ addInletBC()

◆ addMaterials()

◆ addMomentumBoussinesqKernels()

◆ addMomentumFluxKernels()

◆ addMomentumFrictionKernels()

◆ addMomentumGravityKernels()

◆ addMomentumPressureKernels()

◆ addMomentumTimeKernels()

◆ addNonPorousMediumSpeedMaterial()

◆ addOutletBC()

◆ addPorousMediumSpeedMaterial()

◆ addPostprocessors()

◆ addPressureCorrectionKernels()

◆ addRhieChowUserObjects()

◆ addSeparatorBC()

◆ addSolverVariables()

◆ addUserObjects()

◆ addWallsBC()

◆ compressibility()

◆ densityName()

◆ dynamicViscosityName()

◆ getAdditionalRMParams()

◆ getCoupledTurbulencePhysics()

◆ getFlowVariableName()

◆ getFluidTemperatureName()

◆ getFluxInletDirections()

◆ getFluxInletPPs()

◆ getHydraulicSeparators()

◆ getInletBoundaries()

◆ getLinearFrictionCoefName()

◆ getMomentumAdvectionFaceInterpolationMethod()

◆ getMomentumFaceInterpolationMethod()

◆ getNumberAlgebraicGhostingLayersNeeded()

◆ getOutletBoundaries()

◆ getPorosityFunctorName()

◆ getPressureName()

◆ getVelocityFaceInterpolationMethod()

◆ getVelocityNames()

◆ getWallBoundaries()

◆ gravityVector()

◆ hasFlowEquations()

◆ hasForchheimerFriction()

◆ hasTurbulencePhysics()

◆ initializePhysicsAdditional()

◆ porousMediumTreatment()

◆ rhieChowUOName()

◆ usingNavierStokesFVSyntax()

◆ validParams()

Member Data Documentation

◆ _compressibility

◆ _define_variables

◆ _density_gravity_name

◆ _density_name

◆ _dynamic_viscosity_name

◆ _flow_porosity_functor_name

◆ _fluid_temperature_name

◆ _flux_inlet_directions

◆ _flux_inlet_pps

◆ _friction_blocks

◆ _friction_coeffs

◆ _friction_types