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

#include <WCNSLinearFVFlowPhysics.h>

Inheritance diagram for WCNSLinearFVFlowPhysics:

[legend]

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	WCNSLinearFVFlowPhysics (const InputParameters &parameters)

virtual std::vector< UserObjectName >	getSuppliedUserObjects () const override

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

bool	addAxisymmetricViscousSourceEnabled () const
	Whether the cylindrical viscous source helper is enabled.

bool	includeSymmetrizedViscousStress () const
	Whether to include the symmetrized contribution in the viscous stress.

bool	includeIsotropicStress () const
	Whether to include the isotropic viscous stress contribution.

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

const NonlinearVariableName &	getPressureName () const

const NonlinearVariableName &	getFluidTemperatureName () const

MooseFunctorName	getPorosityFunctorName (const bool smoothed) const

const MooseEnum &	compressibility () const
	Return the compressibility of the flow equations selected.

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

RealVectorValue	gravityVector () const
	Return the gravity vector.

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

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

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

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

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

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

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

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

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

NS::MomentumInletTypes	inletBoundaryType (const BoundaryName &boundary_name) const
	Get the type of the inlet BC.

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

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

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

virtual void	act () override final

void	addBlocks (const std::vector< SubdomainName > &blocks)

void	addBlocksById (const std::vector< SubdomainID > &block_ids)

const std::vector< SubdomainName > &	blocks () const

bool	checkBlockRestrictionIdentical (const std::string &object_name, const std::vector< SubdomainName > &blocks, const bool error_if_not_identical=true) const

bool	hasBlocks (const std::vector< SubdomainName > &blocks) const

const T *	getCoupledPhysics (const PhysicsName &phys_name, const bool allow_fail=false) const

const std::vector< T * >	getCoupledPhysics (const bool allow_fail=false) const

unsigned int	dimension () const

const ActionComponent &	getActionComponent (const ComponentName &comp_name) const

void	checkComponentType (const ActionComponent &component) const

virtual void	addComponent (const ActionComponent &component)

const std::vector< VariableName > &	solverVariableNames () const

const std::vector< VariableName > &	auxVariableNames () const

void	timedAct ()

MooseObjectName	uniqueActionName () const

const std::string &	specificTaskName () const

const std::set< std::string > &	getAllTasks () const

void	appendTask (const std::string &task)

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

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

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () const

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

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

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

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

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

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

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

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

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

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

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

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

std::string	messagePrefix (const bool hit_prefix=true) const

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

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

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

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

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

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

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

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

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

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

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

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

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) const

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

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

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

PerfGraph &	perfGraph ()

void	assertParamDefined (const std::string &libmesh_dbg_var(param)) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

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

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::string	unique_action_name_param

static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Member Functions
virtual void	initializePhysicsAdditional () override

virtual void	actOnAdditionalTasks () override

virtual void	addInitialConditions () override

virtual void	addFVBCs () override

virtual void	addMaterials () override

virtual void	addPostprocessors () override

void	addAxisymmetricViscousSource ()
	Adds the cylindrical source kernel for the radial momentum equation when requested and valid.

virtual void	addAxisymmetricViscousSourceKernel (const std::vector< SubdomainName > &, unsigned int)
	Derived classes must override this hook to add the actual object that implements the axisymmetric viscous source term for their formulation once the helper has determined the relevant blocks and radial component.

void	addPorousMediumSpeedMaterial ()
	Add material to define the local speed in porous medium flows.

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

void	addFluidPropertiesFunctorMaterial ()
	Function which adds the general functor fluid properties functor material to define fluid functor material property.

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

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

const WCNSFVTurbulencePhysicsBase *	getCoupledTurbulencePhysics () const
	Find the turbulence physics.

std::vector< SubdomainName >	getAxisymmetricRZBlocks () const
	Return the set of blocks restricted to an RZ coordinate system.

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

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

void	addFVAdvectedInterpolationMethod (const MooseEnum &interpolation_method)
	Add the FVInterpolationMethod object for an advected interpolation method if absent.

void	assertParamDefined (const std::string &param) const

bool	isTransient () const

Factory &	getFactory ()

Factory &	getFactory () const

virtual FEProblemBase &	getProblem ()

virtual const FEProblemBase &	getProblem () const

void	prepareCopyVariablesFromMesh () const

void	copyVariablesFromMesh (const std::vector< VariableName > &variables_to_copy, bool are_nonlinear=true)

std::string	prefix () const

void	addUserObject (const std::string &uo_type, const std::string &uo_name, InputParameters &params)

void	saveSolverVariableName (const VariableName &var_name)

void	saveAuxVariableName (const VariableName &var_name)

bool	variableExists (const VariableName &var_name, bool error_if_aux) const

bool	solverVariableExists (const VariableName &var_name) const

const SolverSystemName &	getSolverSystem (unsigned int variable_index) const

const SolverSystemName &	getSolverSystem (const VariableName &variable_name) const

void	addRequiredPhysicsTask (const std::string &task)

void	assignBlocks (InputParameters &params, const std::vector< SubdomainName > &blocks) const

bool	allMeshBlocks (const std::vector< SubdomainName > &blocks) const

bool	allMeshBlocks (const std::set< SubdomainName > &blocks) const

std::set< SubdomainID >	getSubdomainIDs (const std::set< SubdomainName > &blocks) const

std::vector< std::string >	getSubdomainNamesAndIDs (const std::set< SubdomainID > &blocks) const

void	addPetscPairsToPetscOptions (const std::vector< std::pair< MooseEnumItem, std::string > > &petsc_pair_options)

bool	isVariableFV (const VariableName &var_name) const

bool	isVariableScalar (const VariableName &var_name) const

bool	shouldCreateVariable (const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool error_if_aux)

bool	shouldCreateIC (const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool ic_is_default_ic, const bool error_if_already_defined) const

bool	shouldCreateTimeDerivative (const VariableName &var_name, const std::vector< SubdomainName > &blocks, const bool error_if_already_defined) const

void	reportPotentiallyMissedParameters (const std::vector< std::string > &param_names, const std::string &object_type, const std::string &object_name="") const

virtual void	checkIntegrity () const

void	associateWithParameter (const std::string &param_name, InputParameters &params) const

void	associateWithParameter (const InputParameters &from_params, const std::string &param_name, InputParameters &params) const

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const

std::string	timedSectionName (const std::string &section_name) const

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

void	checkParamsBothSetOrNotSet (const std::string &param1, const std::string &param2) const

void	checkSecondParamSetOnlyIfFirstOneTrue (const std::string &param1, const std::string &param2) const

void	checkSecondParamSetOnlyIfFirstOneSet (const std::string &param1, const std::string &param2) const

void	checkSecondParamNotSetIfFirstOneSet (const std::string &param1, const std::string &param2) const

void	checkVectorParamsSameLength (const std::string &param1, const std::string &param2) const

void	checkVectorParamAndMultiMooseEnumLength (const std::string &param1, const std::string &param2) const

void	checkTwoDVectorParamsSameLength (const std::string &param1, const std::string &param2) const

void	checkVectorParamsNoOverlap (const std::vector< std::string > &param_vecs) const

void	checkTwoDVectorParamsNoRespectiveOverlap (const std::vector< std::string > &param_vecs) const

void	checkTwoDVectorParamInnerSameLengthAsOneDVector (const std::string &param1, const std::string &param2) const

void	checkTwoDVectorParamMultiMooseEnumSameLength (const std::string &param1, const std::string &param2, const bool error_for_param2) const

void	checkVectorParamNotEmpty (const std::string &param1) const

void	checkVectorParamsSameLengthIfSet (const std::string &param1, const std::string &param2, const bool ignore_empty_default_param2=false) const

void	checkVectorParamLengthSameAsCombinedOthers (const std::string &param1, const std::string &param2, const std::string &param3) const

void	checkBlockwiseConsistency (const std::string &block_param_name, const std::vector< std::string > &parameter_names) const

bool	parameterConsistent (const InputParameters &other_param, const std::string &param_name) const

void	warnInconsistent (const InputParameters &parameters, const std::string &param_name) const

void	errorDependentParameter (const std::string &param1, const std::string &value_not_set, const std::vector< std::string > &dependent_params) const

void	errorInconsistentDependentParameter (const std::string &param1, const std::string &value_set, const std::vector< std::string > &dependent_params) const

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
const TagName	_pressure_tag = "p_tag"
	Name of the vector to hold pressure momentum equation contributions.

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

const bool	_add_rz_viscous_source
	Whether to automatically add the cylindrical viscous source term.

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

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

const bool	_porous_medium_treatment
	Whether to use the porous medium treatment.

const MooseFunctorName	_porosity_name
	Name of the porosity functor.

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

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

const NonlinearVariableName	_pressure_name
	Pressure name.

const NonlinearVariableName	_fluid_temperature_name
	Fluid temperature name.

const MooseFunctorName	_density_name
	Name of the density material property.

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

const MooseFunctorName	_dynamic_viscosity_name
	Name of the dynamic viscosity material property.

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

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

UserObjectName	_rc_uo_name
	name of the Rhie Chow user object

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

std::vector< SolverSystemName >	_system_names

std::vector< unsigned int >	_system_numbers

const bool	_verbose

const MooseEnum &	_preconditioning

std::vector< SubdomainName >	_blocks

std::string	_registered_identifier

std::string	_specific_task_name

std::set< std::string >	_all_tasks

ActionWarehouse &	_awh

const std::string &	_current_task

std::shared_ptr< MooseMesh > &	_mesh

std::shared_ptr< MooseMesh > &	_displaced_mesh

std::shared_ptr< FEProblemBase > &	_problem

PerfID	_act_timer

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

Private Member Functions
virtual void	addSolverVariables () override

virtual void	addFVInterpolationMethods () override

virtual void	addFVKernels () override

virtual void	addUserObjects () override

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

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

void	addMomentumFluxKernels ()

virtual void	addMomentumPressureKernels () override

virtual void	addMomentumGravityKernels () override

virtual void	addMomentumFrictionKernels () override

virtual void	addMomentumBoussinesqKernels () override

virtual void	addInletBC () override
	Functions adding boundary conditions for the flow simulation.

virtual void	addOutletBC () override

virtual void	addWallsBC () override

virtual void	addSeparatorBC () override

virtual bool	hasForchheimerFriction () const override
	Return whether a Forchheimer friction model is in use.

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

virtual void	addFunctorMaterials () override

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

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

virtual void	addRelationshipManagers (Moose::RelationshipManagerType input_rm_type) override

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

virtual void	addRelationshipManagers (Moose::RelationshipManagerType when_type)

void	initializePhysics ()

virtual void	checkIntegrityEarly () const

virtual void	addAuxiliaryVariables ()

virtual void	addFEKernels ()

virtual void	addNodalKernels ()

virtual void	addDiracKernels ()

virtual void	addDGKernels ()

virtual void	addScalarKernels ()

virtual void	addInterfaceKernels ()

virtual void	addFVInterfaceKernels ()

virtual void	addFEBCs ()

virtual void	addNodalBCs ()

virtual void	addPeriodicBCs ()

virtual void	addFunctions ()

virtual void	addAuxiliaryKernels ()

virtual void	addCorrectors ()

virtual void	addMultiApps ()

virtual void	addTransfers ()

virtual void	addVectorPostprocessors ()

virtual void	addReporters ()

virtual void	addOutputs ()

virtual void	addPreconditioning ()

virtual void	addExecutioner ()

virtual void	addExecutors ()

void	checkRequiredTasks () const

bool	addRelationshipManager (Moose::RelationshipManagerType input_rm_type, const InputParameters &moose_object_pars, std::string rm_name, Moose::RelationshipManagerType rm_type, Moose::RelationshipManagerInputParameterCallback rm_input_parameter_func, Moose::RMSystemType sys_type=Moose::RMSystemType::NONE)

virtual std::string	meshPropertyPrefix (const std::string &data_name) const

const RestartableDataValue &	getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const

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

const T &	forwardGetParam (const std::string &param_name) const

const InputParameters &	forwardParameters () const

bool	forwardIsParamSetByUser (const std::string &param_name) const

bool	forwardIsParamValid (const std::string &param_name) const

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

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

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

const std::string &	forwardType () const

virtual const std::string &	forwardName () const

const std::vector< SubdomainName > &	forwardBlocks () const

Static Private Member Functions
static const hit::Node *	getHitNode (const InputParameters &params)

static std::string	messagePrefix (const InputParameters &params, const bool hit_prefix)

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

MooseEnum	_is_transient

std::vector< VariableName >	_solver_var_names

std::vector< VariableName >	_aux_var_names

unsigned int	_dim

std::set< std::string >	_required_tasks

const ParallelParamObject &	_parent

MooseApp &	_meta_data_app

const MooseObject *const	_meta_data_object

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

const PhysicsBase *const	_customer_class

Detailed Description

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

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

  : WCNSFVFlowPhysicsBase(parameters),
    _non_orthogonal_correction(getParam<bool>("orthogonality_correction"))
{
  addRequiredPhysicsTask("add_interpolation_method_physics");
 
  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 289 of file WCNSFVFlowPhysicsBase.C.

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

◆ addAxisymmetricViscousSource()

void WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource ( )

protectedinherited

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

Definition at line 570 of file WCNSFVFlowPhysicsBase.C.

{
  if (!_has_flow_equations || !_add_rz_viscous_source)
    return;
 
  const auto rz_blocks = getAxisymmetricRZBlocks();
  if (rz_blocks.empty())
    return;
 
  const auto radial_index = getProblem().mesh().getAxisymmetricRadialCoord();
  addAxisymmetricViscousSourceKernel(rz_blocks, radial_index);
}

Referenced by WCNSFVFlowPhysics::addFVKernels().

◆ addAxisymmetricViscousSourceEnabled()

bool WCNSFVFlowPhysicsBase::addAxisymmetricViscousSourceEnabled ( ) const

inlineinherited

Whether the cylindrical viscous source helper is enabled.

Definition at line 39 of file WCNSFVFlowPhysicsBase.h.

39{ return _add_rz_viscous_source; }

Referenced by WCNSFVTurbulencePhysics::addAxisymmetricTurbulentViscousSource().

◆ addAxisymmetricViscousSourceKernel()

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

inlineprotectedvirtualinherited

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

Reimplemented in WCNSFVFlowPhysics.

Definition at line 127 of file WCNSFVFlowPhysicsBase.h.

129 {

130 }

Referenced by WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource().

◆ addFluidPropertiesFunctorMaterial()

void WCNSFVFlowPhysicsBase::addFluidPropertiesFunctorMaterial ( )

protectedinherited

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

Definition at line 354 of file WCNSFVFlowPhysicsBase.C.

{
  // Not very future-proof but it works
  const bool use_ad = !dynamic_cast<WCNSLinearFVFlowPhysics *>(this);
  const std::string class_name =
      use_ad ? "GeneralFunctorFluidProps" : "NonADGeneralFunctorFluidProps";
  InputParameters params = getFactory().getValidParams(class_name);
  assignBlocks(params, _blocks);
 
  params.set<MooseFunctorName>(NS::pressure) = _pressure_name;
  params.set<MooseFunctorName>(NS::T_fluid) = _fluid_temperature_name;
  params.set<MooseFunctorName>(NS::speed) = NS::speed;
  params.applySpecificParameters(parameters(), {NS::fluid, NS::density, "mu_rampdown"});
  if (!MooseUtils::parsesToReal(_density_name))
    params.set<bool>("force_define_density") = true;
  if (!_porous_medium_treatment)
  {
    params.set<MooseFunctorName>(NS::porosity) = "1";
    params.set<MooseFunctorName>("characteristic_length") = "1";
  }
  else
    // not implemented yet
    paramInfo(
        NS::fluid,
        "Specifying the fluid properties user object does not define the GeneralFunctorFluidProps "
        "when using the porous medium treatment. You have to define this object in the input");
 
  // Dynamic pressure
  params.set<bool>("solving_for_dynamic_pressure") = _solve_for_dynamic_pressure;
  if (_solve_for_dynamic_pressure)
  {
    params.set<Point>("reference_pressure_point") = getParam<Point>("reference_pressure_point");
    if (!isParamSetByUser("reference_pressure_point"))
      paramWarning("reference_pressure_point",
                   "Default value of (0,0,0) used. If this point is outside the flow domain, the "
                   "simulation will error");
    params.set<Real>("reference_pressure") = getParam<Real>("reference_pressure");
  }
  params.set<Point>("gravity") = getParam<RealVectorValue>("gravity");
 
  if (!_porous_medium_treatment)
    getProblem().addFunctorMaterial(class_name, prefix() + "functor_fluidprops", params);
}

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

◆ addFunctorMaterials()

void WCNSLinearFVFlowPhysics::addFunctorMaterials ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 708 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);
      }
  }
}

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

{
  const std::string method_name = interpolation_method;
  if (getProblem().hasFVInterpolationMethod(method_name))
    return;
 
  const auto method_type = NS::fvAdvectedInterpolationMethodType(interpolation_method);
 
  InputParameters params = getFactory().getValidParams(method_type);
  getProblem().addFVInterpolationMethod(method_type, method_name, params);
}

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

◆ addFVBCs()

void WCNSFVFlowPhysicsBase::addFVBCs ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 297 of file WCNSFVFlowPhysicsBase.C.

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

◆ addFVInterpolationMethods()

void WCNSLinearFVFlowPhysics::addFVInterpolationMethods ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 101 of file WCNSLinearFVFlowPhysics.C.

{
  if (!_has_flow_equations || isParamValid("momentum_advection_interpolation_method_name"))
    return;
 
  addFVAdvectedInterpolationMethod(getParam<MooseEnum>("momentum_advection_interpolation"));
}

◆ addFVKernels()

void WCNSLinearFVFlowPhysics::addFVKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 190 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 399 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("FVFunctionIC");
  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().addFVInitialCondition(
          "FVFunctionIC", 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().addFVInitialCondition("FVFunctionIC", 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 425 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 306 of file WCNSFVFlowPhysicsBase.C.

{
  if (hasForchheimerFriction() || _porous_medium_treatment)
    addPorousMediumSpeedMaterial();
  else
    addNonPorousMediumSpeedMaterial();
 
  if (isParamValid(NS::fluid))
    addFluidPropertiesFunctorMaterial();
}

◆ addMomentumBoussinesqKernels()

void WCNSLinearFVFlowPhysics::addMomentumBoussinesqKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 398 of file WCNSLinearFVFlowPhysics.C.

{
  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);
  }
}

Referenced by addFVKernels().

◆ addMomentumFluxKernels()

void WCNSLinearFVFlowPhysics::addMomentumFluxKernels ( )

private

Definition at line 268 of file WCNSLinearFVFlowPhysics.C.

{
  const auto momentum_advection_method_name =
      isParamValid("momentum_advection_interpolation_method_name")
          ? getParam<InterpolationMethodName>("momentum_advection_interpolation_method_name")
          : InterpolationMethodName(
                std::string(getParam<MooseEnum>("momentum_advection_interpolation")));
 
  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);
  if (!_turbulence_physics)
    params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
  else
    params.set<MooseFunctorName>(NS::mu) = NS::mu_eff;
 
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<InterpolationMethodName>("advected_interp_method_name") =
      momentum_advection_method_name;
  params.set<bool>("use_nonorthogonal_correction") = _non_orthogonal_correction;
  params.set<bool>("use_deviatoric_terms") = includeSymmetrizedViscousStress();
 
  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);
  }
}

Referenced by addFVKernels().

◆ addMomentumFrictionKernels()

void WCNSLinearFVFlowPhysics::addMomentumFrictionKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 324 of file WCNSLinearFVFlowPhysics.C.

{
  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("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);
    }
  }
}

Referenced by addFVKernels().

◆ addMomentumGravityKernels()

void WCNSLinearFVFlowPhysics::addMomentumGravityKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 374 of file WCNSLinearFVFlowPhysics.C.

{
  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);
      }
  }
}

Referenced by addFVKernels().

◆ addMomentumPressureKernels()

void WCNSLinearFVFlowPhysics::addMomentumPressureKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 306 of file WCNSLinearFVFlowPhysics.C.

{
  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);
  }
}

Referenced by addFVKernels().

◆ addMomentumTimeKernels()

void WCNSLinearFVFlowPhysics::addMomentumTimeKernels ( )

overrideprivatevirtual

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

Implements WCNSFVFlowPhysicsBase.

Definition at line 250 of file WCNSLinearFVFlowPhysics.C.

{
  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);
  }
}

Referenced by addFVKernels().

◆ addNonPorousMediumSpeedMaterial()

void WCNSFVFlowPhysicsBase::addNonPorousMediumSpeedMaterial ( )

protectedinherited

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

Definition at line 336 of file WCNSFVFlowPhysicsBase.C.

{
  // Not very future-proof but it works
  const bool use_ad = !dynamic_cast<WCNSLinearFVFlowPhysics *>(this);
  const std::string class_name =
      use_ad ? "ADVectorMagnitudeFunctorMaterial" : "VectorMagnitudeFunctorMaterial";
  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().addFunctorMaterial(class_name, prefix() + "ins_speed_material", params);
}

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

◆ addOutletBC()

void WCNSLinearFVFlowPhysics::addOutletBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 508 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 318 of file WCNSFVFlowPhysicsBase.C.

{
  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().addFunctorMaterial(
      "PINSFVSpeedFunctorMaterial", prefix() + "pins_speed_material", params);
}

Referenced by WCNSFVFlowPhysicsBase::addMaterials().

◆ addPostprocessors()

void WCNSFVFlowPhysicsBase::addPostprocessors ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 460 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 221 of file WCNSLinearFVFlowPhysics.C.

{
  {
    std::string kernel_type = "LinearFVPressureCorrectionDiffusion";
    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);
  }
}

Referenced by addFVKernels().

◆ addRhieChowUserObjects()

void WCNSLinearFVFlowPhysics::addRhieChowUserObjects ( )

overrideprivatevirtual

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

Implements WCNSFVFlowPhysicsBase.

Definition at line 650 of file WCNSLinearFVFlowPhysics.C.

{
  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;
  params.set<MooseEnum>("pressure_projection_method") =
      getParam<MooseEnum>("pressure_projection_method");
  params.set<MooseEnum>("pressure_diffusion_interpolation") =
      getParam<MooseEnum>("pressure_diffusion_interpolation");
 
  addUserObject(object_type, rhieChowUOName(), params);
}

Referenced by addUserObjects().

◆ addSeparatorBC()

virtual void WCNSLinearFVFlowPhysics::addSeparatorBC ( )

inlineoverrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 60 of file WCNSLinearFVFlowPhysics.h.

60{}

◆ addSolverVariables()

void WCNSLinearFVFlowPhysics::addSolverVariables ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 130 of file WCNSLinearFVFlowPhysics.C.

{
  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 642 of file WCNSLinearFVFlowPhysics.C.

{
  mooseAssert(!_porous_medium_treatment, "Not implemented");
  // Rhie Chow user object for interpolation velocities
  addRhieChowUserObjects();
}

◆ addWallsBC()

void WCNSLinearFVFlowPhysics::addWallsBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 559 of file WCNSLinearFVFlowPhysics.C.

{
  const std::string u_names[3] = {"u", "v", "w"};
  bool has_symmetry_bc = false;
 
  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 if (momentum_wall_type == "symmetry")
    {
      has_symmetry_bc = true;
      {
        const std::string bc_type = "LinearFVVelocitySymmetryBC";
        InputParameters params = getFactory().getValidParams(bc_type);
        params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
        for (unsigned int d = 0; d < dimension(); ++d)
          params.set<SolverVariableName>(u_names[d]) = _velocity_names[d];
 
        for (const auto d : make_range(dimension()))
        {
          params.set<LinearVariableName>("variable") = _velocity_names[d];
          params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
          getProblem().addLinearFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
        }
      }
      {
        const std::string bc_type = "LinearFVPressureSymmetryBC";
        InputParameters params = getFactory().getValidParams(bc_type);
        params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
        params.set<LinearVariableName>("variable") = _pressure_name;
        params.set<MooseFunctorName>("HbyA_flux") = "HbyA";
        getProblem().addLinearFVBC(bc_type, _pressure_name + "_" + boundary_name, params);
      }
    }
    else
      mooseError("Unsupported wall boundary condition type: " + std::string(momentum_wall_type));
  }
 
  if (getParam<bool>("pressure_two_term_bc_expansion"))
  {
    if (!has_symmetry_bc)
    {
      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);
    }
    else
      for (const auto & [boundary_name, momentum_wall_type] : _momentum_wall_types)
        if (momentum_wall_type != "symmetry")
        {
          const std::string bc_type = "LinearFVExtrapolatedPressureBC";
          InputParameters params = getFactory().getValidParams(bc_type);
          params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
          params.set<LinearVariableName>("variable") = _pressure_name;
          params.set<bool>("use_two_term_expansion") = true;
          getProblem().addLinearFVBC(
              bc_type, _pressure_name + "_extrapolation_" + boundary_name, params);
        }
  }
}

◆ compressibility()

const MooseEnum & WCNSFVFlowPhysicsBase::compressibility ( ) const

inlineinherited

Return the compressibility of the flow equations selected.

Definition at line 53 of file WCNSFVFlowPhysicsBase.h.

53{ return _compressibility; }

◆ densityName()

const MooseFunctorName & WCNSFVFlowPhysicsBase::densityName ( ) const

inlineinherited

Return the name of the density functor.

Definition at line 59 of file WCNSFVFlowPhysicsBase.h.

59{ return _density_name; }

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

◆ dynamicViscosityName()

const MooseFunctorName & WCNSFVFlowPhysicsBase::dynamicViscosityName ( ) const

inlineinherited

Return the name of the dynamic viscosity functor.

Definition at line 61 of file WCNSFVFlowPhysicsBase.h.

61{ return _dynamic_viscosity_name; }

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

◆ getAdditionalRMParams()

InputParameters NavierStokesPhysicsBase::getAdditionalRMParams ( ) const

overrideprotectedvirtualinherited

Parameters to change or add relationship managers.

Reimplemented from PhysicsBase.

Definition at line 56 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;
}

◆ getAxisymmetricRZBlocks()

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

protectedinherited

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

Definition at line 539 of file WCNSFVFlowPhysicsBase.C.

{
  std::vector<SubdomainName> rz_blocks;
  const auto & mesh = getProblem().mesh();
 
  const bool use_all_blocks =
      _blocks.empty() || allMeshBlocks(_blocks) ||
      std::find(_blocks.begin(), _blocks.end(), "ANY_BLOCK_ID") != _blocks.end();
 
  std::vector<SubdomainID> block_ids;
  if (use_all_blocks)
  {
    const auto & mesh_blocks = mesh.meshSubdomains();
    block_ids.insert(block_ids.end(), mesh_blocks.begin(), mesh_blocks.end());
  }
  else
    block_ids = mesh.getSubdomainIDs(_blocks);
 
  for (const auto subdomain_id : block_ids)
    if (mesh.getCoordSystem(subdomain_id) == Moose::COORD_RZ)
    {
      auto name = mesh.getSubdomainName(subdomain_id);
      if (name.empty())
        name = Moose::stringify(subdomain_id);
      rz_blocks.push_back(name);
    }
 
  return rz_blocks;
}

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

◆ getCoupledTurbulencePhysics()

const WCNSFVTurbulencePhysicsBase * WCNSFVFlowPhysicsBase::getCoupledTurbulencePhysics ( ) const

protectedinherited

Find the turbulence physics.

Definition at line 506 of file WCNSFVFlowPhysicsBase.C.

{
  // User passed it, just use that
  if (isParamValid("coupled_turbulence_physics"))
    return getCoupledPhysics<WCNSFVTurbulencePhysicsBase>(
        getParam<PhysicsName>("coupled_turbulence_physics"));
  // Look for any physics of the right type, and check the block restriction
  else
  {
    const auto all_turbulence_physics = getCoupledPhysics<const WCNSFVTurbulencePhysicsBase>(true);
    for (const auto physics : all_turbulence_physics)
    {
      if (checkBlockRestrictionIdentical(
              physics->name(), physics->blocks(), /*error_if_not_identical=*/false))
        return physics;
      else if (_verbose)
        mooseInfoRepeated("Detected Turbulence Physics '" + physics->name() +
                          "' with an incompatible block restriction. It will thus not be coupled "
                          "to this flow equations physics");
    }
  }
  // Did not find one
  return nullptr;
}

Referenced by WCNSFVFlowPhysicsBase::actOnAdditionalTasks().

◆ getFlowVariableName()

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

protectedinherited

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

Definition at line 480 of file WCNSFVFlowPhysicsBase.C.

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

48{ return _fluid_temperature_name; }

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

◆ getFluxInletDirections()

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

inlineinherited

Get the inlet direction if using a flux inlet.

Definition at line 89 of file WCNSFVFlowPhysicsBase.h.

89{ return _flux_inlet_directions; }

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

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

◆ getFluxInletPPs()

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

inlineinherited

Get the inlet flux postprocessor if using a flux inlet.

Definition at line 91 of file WCNSFVFlowPhysicsBase.h.

91{ return _flux_inlet_pps; }

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

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

◆ getHydraulicSeparators()

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

inlineinherited

Get the hydraulic separator boundaries.

Definition at line 81 of file WCNSFVFlowPhysicsBase.h.

81{ return _hydraulic_separators; }

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergySeparatorBC().

◆ getInletBoundaries()

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

inlineinherited

Get the inlet boundaries.

Definition at line 75 of file WCNSFVFlowPhysicsBase.h.

75{ return _inlet_boundaries; }

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

◆ getLinearFrictionCoefName()

virtual MooseFunctorName WCNSLinearFVFlowPhysics::getLinearFrictionCoefName ( ) const

inlineoverrideprivatevirtual

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

Implements WCNSFVFlowPhysicsBase.

Definition at line 67 of file WCNSLinearFVFlowPhysics.h.

  {
    mooseError("Not implemented");
  }

◆ getMomentumAdvectionFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getMomentumAdvectionFaceInterpolationMethod ( ) const

inlineinherited

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

Definition at line 65 of file WCNSFVFlowPhysicsBase.h.

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

  {
    return _momentum_face_interpolation;
  }

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

◆ getNumberAlgebraicGhostingLayersNeeded()

unsigned short WCNSLinearFVFlowPhysics::getNumberAlgebraicGhostingLayersNeeded ( ) const

overrideprivatevirtual

Return the number of algebraic ghosting layers needed.

Reimplemented from WCNSFVFlowPhysicsBase.

Definition at line 734 of file WCNSLinearFVFlowPhysics.C.

{
  return 1;
}

◆ getOutletBoundaries()

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

inlineinherited

Get the outlet boundaries.

Definition at line 77 of file WCNSFVFlowPhysicsBase.h.

77{ return _outlet_boundaries; }

WCNSFVFlowPhysicsBase::_outlet_boundaries

const std::vector< BoundaryName > _outlet_boundaries

Boundaries with a flow outlet specified on them.

Definition WCNSFVFlowPhysicsBase.h:230

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

◆ getPorosityFunctorName()

MooseFunctorName WCNSFVFlowPhysicsBase::getPorosityFunctorName ( const bool smoothed ) const

inherited

Definition at line 497 of file WCNSFVFlowPhysicsBase.C.

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

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

◆ getPressureName()

const NonlinearVariableName & WCNSFVFlowPhysicsBase::getPressureName ( ) const

inlineinherited

Definition at line 47 of file WCNSFVFlowPhysicsBase.h.

47{ return _pressure_name; }

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC(), WCNSFVFluidHeatTransferPhysics::addEnergyWallBC(), WCNSFVFluidHeatTransferPhysicsBase::addInitialConditions(), WCNSFVFluidHeatTransferPhysics::addMaterials(), WCNSLinearFVFluidHeatTransferPhysics::addMaterials(), and WCNSFVFlowPhysicsBase::getFlowVariableName().

◆ getSuppliedUserObjects()

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

overridevirtual

Reimplemented from PhysicsBase.

Definition at line 702 of file WCNSLinearFVFlowPhysics.C.

{
  return {rhieChowUOName()};
}

◆ getVelocityFaceInterpolationMethod()

const MooseEnum & WCNSFVFlowPhysicsBase::getVelocityFaceInterpolationMethod ( ) const

inlineinherited

Get the face interpolation method for velocity.

Definition at line 63 of file WCNSFVFlowPhysicsBase.h.

63{ return _velocity_interpolation; }

WCNSFVFlowPhysicsBase::_velocity_interpolation

const MooseEnum _velocity_interpolation

The velocity face interpolation method for advecting other quantities.

Definition WCNSFVFlowPhysicsBase.h:211

Referenced by WCNSFVTwoPhaseMixturePhysics::addAdvectionSlipTerm().

◆ getVelocityNames()

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

inlineinherited

To interface with other Physics.

Definition at line 46 of file WCNSFVFlowPhysicsBase.h.

46{ return _velocity_names; }

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

◆ getWallBoundaries()

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

inlineinherited

Get the wall boundaries.

Definition at line 79 of file WCNSFVFlowPhysicsBase.h.

79{ return _wall_boundaries; }

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

◆ gravityVector()

RealVectorValue WCNSFVFlowPhysicsBase::gravityVector ( ) const

inlineinherited

Return the gravity vector.

Definition at line 57 of file WCNSFVFlowPhysicsBase.h.

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

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

◆ hasFlowEquations()

bool WCNSFVFlowPhysicsBase::hasFlowEquations ( ) const

inlineinherited

Whether the physics is actually creating the flow equations.

Definition at line 37 of file WCNSFVFlowPhysicsBase.h.

37{ return _has_flow_equations; }

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

◆ hasForchheimerFriction()

virtual bool WCNSLinearFVFlowPhysics::hasForchheimerFriction ( ) const

inlineoverrideprivatevirtual

Return whether a Forchheimer friction model is in use.

Implements WCNSFVFlowPhysicsBase.

Definition at line 62 of file WCNSLinearFVFlowPhysics.h.

62{ return false; };

◆ hasTurbulencePhysics()

bool WCNSFVFlowPhysicsBase::hasTurbulencePhysics ( ) const

inlineprotectedinherited

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

Definition at line 156 of file WCNSFVFlowPhysicsBase.h.

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

Referenced by WCNSFVFlowPhysics::addWallsBC().

◆ includeIsotropicStress()

bool WCNSFVFlowPhysicsBase::includeIsotropicStress ( ) const

inlineinherited

Whether to include the isotropic viscous stress contribution.

Definition at line 43 of file WCNSFVFlowPhysicsBase.h.

43{ return _include_isotropic_viscous_stress; }

WCNSFVFlowPhysicsBase::_include_isotropic_viscous_stress

const bool _include_isotropic_viscous_stress

Whether to include the isotropic viscous stress contribution.

Definition WCNSFVFlowPhysicsBase.h:206

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

◆ includeSymmetrizedViscousStress()

bool WCNSFVFlowPhysicsBase::includeSymmetrizedViscousStress ( ) const

inlineinherited

Whether to include the symmetrized contribution in the viscous stress.

Definition at line 41 of file WCNSFVFlowPhysicsBase.h.

41{ return _include_symmetrized_viscous_stress; }

WCNSFVFlowPhysicsBase::_include_symmetrized_viscous_stress

const bool _include_symmetrized_viscous_stress

Whether to include the symmetrized viscous stress contribution.

Definition WCNSFVFlowPhysicsBase.h:204

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

◆ initializePhysicsAdditional()

void WCNSLinearFVFlowPhysics::initializePhysicsAdditional ( )

overrideprotectedvirtual

Reimplemented from WCNSFVFlowPhysicsBase.

Definition at line 110 of file WCNSLinearFVFlowPhysics.C.

{
  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"};
  }
}

◆ inletBoundaryType()

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

inlineinherited

Get the type of the inlet BC.

Definition at line 83 of file WCNSFVFlowPhysicsBase.h.

  {
    return NS::MomentumInletTypes(
        static_cast<int>(libmesh_map_find(_momentum_inlet_types, boundary_name)));
  }

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC().

◆ porousMediumTreatment()

bool WCNSFVFlowPhysicsBase::porousMediumTreatment ( ) const

inlineinherited

Return whether a porous medium treatment is applied.

Definition at line 55 of file WCNSFVFlowPhysicsBase.h.

55{ return _porous_medium_treatment; }

◆ rhieChowUOName()

const UserObjectName & WCNSFVFlowPhysicsBase::rhieChowUOName ( ) const

inherited

Return the name of the Rhie Chow user object.

Definition at line 532 of file WCNSFVFlowPhysicsBase.C.

{
  mooseAssert(!_rc_uo_name.empty(), "The Rhie-Chow user-object name should be set!");
  return _rc_uo_name;
}

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

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

◆ validParams()

InputParameters WCNSLinearFVFlowPhysics::validParams ( )

static

Definition at line 27 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.set<MooseEnum>("momentum_advection_interpolation") = NS::fvAdvectedInterpolationMethods();
  params.addParam<InterpolationMethodName>(
      "momentum_advection_interpolation_method_name",
      "Name of an externally defined FVInterpolationMethod to use for momentum advection. When "
      "provided, this overrides 'momentum_advection_interpolation'.");
 
  params.addParam<bool>(
      "orthogonality_correction", false, "Whether to use orthogonality correction");
  params.renameParam("orthogonality_correction", "use_nonorthogonal_correction", "");
  params.addParam<MooseEnum>(
      "pressure_diffusion_interpolation",
      NS::fvFaceInterpolationMethods(),
      "The face interpolation method for Ainv in the pressure correction diffusion term.");
  params.addParamNamesToGroup("pressure_diffusion_interpolation", "Numerical scheme");
  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");
 
  // Rhie-Chow
  params.transferParam<MooseEnum>(RhieChowMassFlux::validParams(), "pressure_projection_method");
  params.addParamNamesToGroup("momentum_advection_interpolation_method_name", "Numerical scheme");
 
  return params;
}

Member Data Documentation

◆ _add_rz_viscous_source

const bool WCNSFVFlowPhysicsBase::_add_rz_viscous_source

protectedinherited

Whether to automatically add the cylindrical viscous source term.

Definition at line 176 of file WCNSFVFlowPhysicsBase.h.

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

◆ _compressibility

const MooseEnum WCNSFVFlowPhysicsBase::_compressibility

protectedinherited

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

Definition at line 179 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysics::addFVKernels(), WCNSFVFlowPhysics::addMomentumBoussinesqKernels(), 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 48 of file NavierStokesPhysicsBase.h.

Referenced by WCNSFVTurbulencePhysics::addAuxiliaryVariables(), WCNSLinearFVFluidHeatTransferPhysics::addAuxiliaryVariables(), WCNSFVFlowPhysicsBase::addInitialConditions(), WCNSFVFluidHeatTransferPhysicsBase::addInitialConditions(), WCNSFVScalarTransportPhysicsBase::addInitialConditions(), WCNSFVTurbulencePhysicsBase::addInitialConditions(), WCNSFVFlowPhysics::addSolverVariables(), WCNSFVFluidHeatTransferPhysics::addSolverVariables(), WCNSFVTurbulencePhysics::addSolverVariables(), addSolverVariables(), WCNSLinearFVFluidHeatTransferPhysics::addSolverVariables(), WCNSLinearFVTurbulencePhysics::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 200 of file WCNSFVFlowPhysicsBase.h.

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

◆ _density_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_name

protectedinherited

Name of the density material property.

Definition at line 198 of file WCNSFVFlowPhysicsBase.h.

◆ _dynamic_viscosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_dynamic_viscosity_name

protectedinherited

Name of the dynamic viscosity material property.

Definition at line 202 of file WCNSFVFlowPhysicsBase.h.

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

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

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

◆ _flux_inlet_pps

std::vector<PostprocessorName> WCNSFVFlowPhysicsBase::_flux_inlet_pps

protectedinherited

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

Definition at line 244 of file WCNSFVFlowPhysicsBase.h.

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

◆ _friction_blocks

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

protectedinherited

Subdomains where we want to have volumetric friction.

Definition at line 221 of file WCNSFVFlowPhysicsBase.h.

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

◆ _friction_coeffs

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

protectedinherited

The coefficients used for each item if friction type.

Definition at line 225 of file WCNSFVFlowPhysicsBase.h.

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

◆ _friction_types

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

protectedinherited

The friction correlation types used for each block.

Definition at line 223 of file WCNSFVFlowPhysicsBase.h.

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

Referenced by WCNSFVFlowPhysicsBase::addAxisymmetricViscousSource(), WCNSFVFlowPhysics::addCorrectors(), addFVInterpolationMethods(), 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 234 of file WCNSFVFlowPhysicsBase.h.

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

◆ _include_isotropic_viscous_stress

const bool WCNSFVFlowPhysicsBase::_include_isotropic_viscous_stress

protectedinherited

Whether to include the isotropic viscous stress contribution.

Definition at line 206 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::includeIsotropicStress().

◆ _include_symmetrized_viscous_stress

const bool WCNSFVFlowPhysicsBase::_include_symmetrized_viscous_stress

protectedinherited

Whether to include the symmetrized viscous stress contribution.

Definition at line 204 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::includeSymmetrizedViscousStress().

◆ _inlet_boundaries

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

protectedinherited

Boundaries with a flow inlet specified on them.

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

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

◆ _momentum_face_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_momentum_face_interpolation

protectedinherited

The momentum face interpolation method for stress terms.

Definition at line 215 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::getMomentumFaceInterpolationMethod().

◆ _momentum_inlet_functors

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

protectedinherited

Functors describing the momentum inlet for each boundary.

Definition at line 249 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_inlet_types

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

protectedinherited

Momentum inlet boundary types.

Definition at line 237 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_outlet_types

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

protectedinherited

Momentum outlet boundary types.

Definition at line 239 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_wall_functors

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

protectedinherited

Functors describing the momentum for each wall boundary.

Definition at line 253 of file WCNSFVFlowPhysicsBase.h.

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

◆ _momentum_wall_types

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

protectedinherited

Momentum wall boundary types.

Definition at line 241 of file WCNSFVFlowPhysicsBase.h.

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

◆ _non_orthogonal_correction

const bool WCNSLinearFVFlowPhysics::_non_orthogonal_correction

private

Whether to use the correction term for non-orthogonality.

Definition at line 75 of file WCNSLinearFVFlowPhysics.h.

Referenced by addMomentumFluxKernels(), and addPressureCorrectionKernels().

◆ _outlet_boundaries

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

protectedinherited

Boundaries with a flow outlet specified on them.

Definition at line 230 of file WCNSFVFlowPhysicsBase.h.

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

◆ _porosity_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_porosity_name

protectedinherited

Name of the porosity functor.

Definition at line 186 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 184 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 251 of file WCNSFVFlowPhysicsBase.h.

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

◆ _pressure_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_pressure_name

protectedinherited

Pressure name.

Definition at line 193 of file WCNSFVFlowPhysicsBase.h.

◆ _pressure_tag

const TagName WCNSFVFlowPhysicsBase::_pressure_tag = "p_tag"

protectedinherited

Name of the vector to hold pressure momentum equation contributions.

Definition at line 171 of file WCNSFVFlowPhysicsBase.h.

◆ _rc_uo_name

UserObjectName WCNSFVFlowPhysicsBase::_rc_uo_name

protectedinherited

name of the Rhie Chow user object

Definition at line 209 of file WCNSFVFlowPhysicsBase.h.

Referenced by WCNSFVFlowPhysicsBase::rhieChowUOName().

◆ _solve_for_dynamic_pressure

const bool WCNSFVFlowPhysicsBase::_solve_for_dynamic_pressure

protectedinherited

Whether we are solving for the total or dynamic pressure.

Definition at line 181 of file WCNSFVFlowPhysicsBase.h.

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

◆ _turbulence_physics

const WCNSFVTurbulencePhysicsBase* WCNSFVFlowPhysicsBase::_turbulence_physics

protectedinherited

Can be set to a coupled turbulence physics.

Definition at line 218 of file WCNSFVFlowPhysicsBase.h.

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

◆ _velocity_interpolation

const MooseEnum WCNSFVFlowPhysicsBase::_velocity_interpolation

protectedinherited

The velocity face interpolation method for advecting other quantities.

Definition at line 211 of file WCNSFVFlowPhysicsBase.h.

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

◆ _velocity_names

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

protectedinherited

Velocity names.

Definition at line 191 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 232 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

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ WCNSLinearFVFlowPhysics()

Member Function Documentation

◆ actOnAdditionalTasks()

◆ addAxisymmetricViscousSource()

◆ addAxisymmetricViscousSourceEnabled()

◆ addAxisymmetricViscousSourceKernel()

◆ addFluidPropertiesFunctorMaterial()

◆ addFunctorMaterials()

◆ addFVAdvectedInterpolationMethod()

◆ addFVBCs()

◆ addFVInterpolationMethods()

◆ addFVKernels()

◆ addInitialConditions()

◆ addInletBC()

◆ addMaterials()

◆ addMomentumBoussinesqKernels()

◆ addMomentumFluxKernels()

◆ addMomentumFrictionKernels()

◆ addMomentumGravityKernels()

◆ addMomentumPressureKernels()

◆ addMomentumTimeKernels()

◆ addNonPorousMediumSpeedMaterial()

◆ addOutletBC()

◆ addPorousMediumSpeedMaterial()

◆ addPostprocessors()

◆ addPressureCorrectionKernels()

◆ addRhieChowUserObjects()

◆ addSeparatorBC()

◆ addSolverVariables()

◆ addUserObjects()

◆ addWallsBC()

◆ compressibility()

◆ densityName()

◆ dynamicViscosityName()

◆ getAdditionalRMParams()

◆ getAxisymmetricRZBlocks()

◆ getCoupledTurbulencePhysics()

◆ getFlowVariableName()

◆ getFluidTemperatureName()

◆ getFluxInletDirections()

◆ getFluxInletPPs()

◆ getHydraulicSeparators()

◆ getInletBoundaries()

◆ getLinearFrictionCoefName()

◆ getMomentumAdvectionFaceInterpolationMethod()

◆ getMomentumFaceInterpolationMethod()

◆ getNumberAlgebraicGhostingLayersNeeded()

◆ getOutletBoundaries()

◆ getPorosityFunctorName()

◆ getPressureName()

◆ getSuppliedUserObjects()

◆ getVelocityFaceInterpolationMethod()

◆ getVelocityNames()

◆ getWallBoundaries()

◆ gravityVector()

◆ hasFlowEquations()

◆ hasForchheimerFriction()

◆ hasTurbulencePhysics()

◆ includeIsotropicStress()

◆ includeSymmetrizedViscousStress()

◆ initializePhysicsAdditional()

◆ inletBoundaryType()

◆ porousMediumTreatment()

◆ rhieChowUOName()

◆ usingNavierStokesFVSyntax()

◆ validParams()

Member Data Documentation