Creates all the objects needed to solve the Navier Stokes mass and momentum equations. More...

#include <WCNSFVFlowPhysics.h>

Inheritance diagram for WCNSFVFlowPhysics:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	WCNSFVFlowPhysics (const InputParameters &parameters)

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

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	checkAtMostOneParamSetByUser (const std::vector< std::string > &params) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

const bool	_porous_medium_treatment
	Whether to use the porous medium treatment.

const MooseFunctorName	_porosity_name
	Name of the porosity functor.

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

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

const NonlinearVariableName	_pressure_name
	Pressure name.

const NonlinearVariableName	_fluid_temperature_name
	Fluid temperature name.

const MooseFunctorName	_density_name
	Name of the density material property.

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

const MooseFunctorName	_dynamic_viscosity_name
	Name of the dynamic viscosity material property.

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

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

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

virtual void	addUserObjects () override

virtual void	addCorrectors () override

void	addMassTimeKernels ()
	Function adding kernels for the time derivative term of the weakly compressible continuity equation.

void	addMassKernels ()
	Function adding kernels for the incompressible continuity equation.

void	addPressurePinKernel ()
	Function adding the pressure constraint.

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

void	addMomentumMixingLengthKernels ()

void	addMomentumAdvectionKernels ()

void	addMomentumPressureKernels () override

void	addMomentumGravityKernels () override

void	addMomentumBoussinesqKernels () override

void	addMomentumFrictionKernels () override

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

void	addOutletBC () override

void	addWallsBC () override

void	addSeparatorBC () override

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

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

void	checkRhieChowFunctorsDefined () const
	Checks that sufficient Rhie Chow coefficients have been defined for the given dimension, used for scalar or temperature advection by auxiliary variables.

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

virtual void	addRelationshipManagers (Moose::RelationshipManagerType input_rm_type) override

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

virtual void	addRelationshipManagers (Moose::RelationshipManagerType when_type)

void	initializePhysics ()

virtual void	checkIntegrityEarly () const

virtual void	addAuxiliaryVariables ()

virtual void	addFEKernels ()

virtual void	addFVInterpolationMethods ()

virtual void	addNodalKernels ()

virtual void	addDiracKernels ()

virtual void	addDGKernels ()

virtual void	addScalarKernels ()

virtual void	addInterfaceKernels ()

virtual void	addFVInterfaceKernels ()

virtual void	addFEBCs ()

virtual void	addNodalBCs ()

virtual void	addPeriodicBCs ()

virtual void	addFunctions ()

virtual void	addAuxiliaryKernels ()

virtual void	addFunctorMaterials ()

virtual void	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 unsigned	_porosity_smoothing_layers
	The number of smoothing layers if that treatment is used on porosity.

UserObjectName	_rc_uo_name
	Name of the user object in charge of computing the Rhie Chow coefficients.

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 mass and momentum equations.

Definition at line 18 of file WCNSFVFlowPhysics.h.

Constructor & Destructor Documentation

◆ WCNSFVFlowPhysics()

WCNSFVFlowPhysics::WCNSFVFlowPhysics ( const InputParameters & parameters )

Definition at line 89 of file WCNSFVFlowPhysics.C.

  : WCNSFVFlowPhysicsBase(parameters),
    _porosity_smoothing_layers(isParamValid("porosity_smoothing_layers")
                                   ? getParam<unsigned short>("porosity_smoothing_layers")
                                   : 0)
{
  _flow_porosity_functor_name = isParamValid("porosity_smoothing_layers") &&
                                        getParam<unsigned short>("porosity_smoothing_layers")
                                    ? NS::smoothed_porosity
                                    : _porosity_name;
 
  // Most likely to be a mistake
  if (getParam<bool>("pin_pressure") &&
      getParam<std::vector<MooseFunctorName>>("pressure_functors").size())
    paramError("pin_pressure", "Cannot pin the pressure if a pressure boundary exists");
 
  // Pressure pin checks
  checkSecondParamSetOnlyIfFirstOneTrue("pin_pressure", "pinned_pressure_type");
  checkSecondParamSetOnlyIfFirstOneTrue("pin_pressure", "pinned_pressure_value");
  if (getParam<bool>("pin_pressure"))
  {
    if ((std::string(getParam<MooseEnum>("pinned_pressure_type")).find("point") !=
         std::string::npos) &&
        !isParamSetByUser("pinned_pressure_point"))
      paramError("pinned_pressure_point",
                 "This parameter must be set to specify the pinned pressure point");
    else if ((std::string(getParam<MooseEnum>("pinned_pressure_type")).find("point") ==
              std::string::npos) &&
             isParamSetByUser("pinned_pressure_point"))
      paramError("pinned_pressure_point",
                 "This parameter should not be given by the user with the corresponding "
                 "pinned_pressure_type setting: " +
                     std::string(getParam<MooseEnum>("pinned_pressure_type")) + ".");
  }
 
  // Porosity correction checks
  checkSecondParamSetOnlyIfFirstOneTrue("porous_medium_treatment", "use_friction_correction");
  checkSecondParamSetOnlyIfFirstOneTrue("use_friction_correction", "consistent_scaling");
  checkSecondParamSetOnlyIfFirstOneTrue("porous_medium_treatment",
                                        "porosity_interface_pressure_treatment");
  if (getParam<MooseEnum>("porosity_interface_pressure_treatment") != "bernoulli")
    errorDependentParameter("porosity_interface_pressure_treatment",
                            "bernoulli",
                            {"pressure_allow_expansion_on_bernoulli_faces",
                             "pressure_drop_sidesets",
                             "pressure_drop_form_factors"});
 
  // Porous media parameters
  checkSecondParamSetOnlyIfFirstOneTrue("porous_medium_treatment", "porosity_smoothing_layers");
}

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

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

overrideprivatevirtual

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 from WCNSFVFlowPhysicsBase.

Definition at line 501 of file WCNSFVFlowPhysics.C.

{
  InputParameters params = getFactory().getValidParams("INSFVMomentumViscousSourceRZ");
  assignBlocks(params, rz_blocks);
  params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<MooseEnum>("momentum_component") = NS::directions[radial_index];
  params.set<bool>("complete_expansion") = includeSymmetrizedViscousStress();
  params.set<NonlinearVariableName>("variable") = _velocity_names[radial_index];
 
  getProblem().addFVKernel("INSFVMomentumViscousSourceRZ",
                           prefix() + "ins_momentum_viscous_source_rz_" +
                               NS::directions[radial_index],
                           params);
}

◆ addCorrectors()

void WCNSFVFlowPhysics::addCorrectors ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 1081 of file WCNSFVFlowPhysics.C.

{
  if (!_has_flow_equations)
    return;
 
  // Pressure pin
  if (getParam<bool>("pin_pressure"))
  {
    const auto pin_type = getParam<MooseEnum>("pinned_pressure_type");
    std::string object_type = "NSPressurePin";
 
    // No need for the user object
    if (pin_type == "point-value" || pin_type == "average")
      return;
 
    // Create the average value postprocessor if needed
    if (pin_type == "average-uo")
    {
      // Volume average by default, but we could do inlet or outlet for example
      InputParameters params = getFactory().getValidParams("ElementAverageValue");
      params.set<std::vector<VariableName>>("variable") = {_pressure_name};
      assignBlocks(params, _blocks);
      params.set<std::vector<OutputName>>("outputs") = {"none"};
      getProblem().addPostprocessor("ElementAverageValue", "ns_pressure_average", params);
    }
 
    InputParameters params = getFactory().getValidParams(object_type);
    if (pin_type == "point-value" || pin_type == "point-value-uo")
      params.set<MooseEnum>("pin_type") = "point-value";
    else
      params.set<MooseEnum>("pin_type") = "average";
 
    params.set<PostprocessorName>("phi0") = getParam<PostprocessorName>("pinned_pressure_value");
    params.set<NonlinearVariableName>("variable") = _pressure_name;
    if (pin_type == "point-value" || pin_type == "point-value-uo")
      params.set<Point>("point") = getParam<Point>("pinned_pressure_point");
    else if (pin_type == "average-uo")
      params.set<PostprocessorName>("pressure_average") = "ns_pressure_average";
 
    getProblem().addUserObject(object_type, prefix() + "ins_mass_pressure_pin", params);
  }
}

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

◆ 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 WCNSLinearFVFlowPhysics::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();
}

◆ addFVKernels()

void WCNSFVFlowPhysics::addFVKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 267 of file WCNSFVFlowPhysics.C.

{
  if (!_has_flow_equations)
    return;
 
  // Mass equation: time derivative
  if (_compressibility == "weakly-compressible" &&
      shouldCreateTimeDerivative(_pressure_name, _blocks, /*error if already defined*/ false))
    addMassTimeKernels();
 
  // Mass equation: divergence of momentum
  addMassKernels();
 
  // Pressure pin
  if (getParam<bool>("pin_pressure"))
    addPressurePinKernel();
 
  // Momentum equation: time derivative
  if (isTransient())
    addMomentumTimeKernels();
 
  // Momentum equation: momentum advection
  addMomentumAdvectionKernels();
 
  // Momentum equation: momentum viscous stress
  addMomentumViscousDissipationKernels();
  addAxisymmetricViscousSource();
 
  // Momentum equation: pressure term
  addMomentumPressureKernels();
 
  // Momentum equation: gravity source term
  addMomentumGravityKernels();
 
  // Momentum equation: friction kernels
  if (_friction_types.size())
    addMomentumFrictionKernels();
 
  // 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 WCNSFVFlowPhysics::addInletBC ( )

overrideprivatevirtual

Functions adding boundary conditions for the incompressible simulation.

These are used for weakly-compressible simulations as well.

Implements WCNSFVFlowPhysicsBase.

Definition at line 724 of file WCNSFVFlowPhysics.C.

{
  // Check the size of the BC parameters
  unsigned int num_velocity_functor_inlets = 0;
  for (const auto & [bdy, momentum_outlet_type] : _momentum_inlet_types)
    if (momentum_outlet_type == "fixed-velocity" || momentum_outlet_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 flux_bc_counter = 0;
  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 = "INSFVInletVelocityBC";
      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<NonlinearVariableName>("variable") = _velocity_names[d];
        params.set<MooseFunctorName>("functor") = momentum_functors[d];
 
        getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + inlet_bdy, params);
      }
      ++velocity_pressure_counter;
    }
    else if (momentum_inlet_type == "fixed-pressure")
    {
      const std::string bc_type = "INSFVOutletPressureBC";
      InputParameters params = getFactory().getValidParams(bc_type);
      params.set<NonlinearVariableName>("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<FunctionName>("function") =
          libmesh_map_find(_momentum_inlet_functors, inlet_bdy)[0];
      params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
 
      getProblem().addFVBC(bc_type, _pressure_name + "_" + inlet_bdy, params);
      ++velocity_pressure_counter;
    }
    else if (momentum_inlet_type == "flux-mass" || momentum_inlet_type == "flux-velocity")
    {
      {
        const std::string bc_type =
            _porous_medium_treatment ? "PWCNSFVMomentumFluxBC" : "WCNSFVMomentumFluxBC";
        InputParameters params = getFactory().getValidParams(bc_type);
 
        if (_flux_inlet_directions.size())
          params.set<Point>("direction") = _flux_inlet_directions[flux_bc_counter];
 
        params.set<MooseFunctorName>(NS::density) = _density_name;
        params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
        params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
        if (_porous_medium_treatment)
          params.set<MooseFunctorName>(NS::porosity) = _porosity_name;
        if (_flux_inlet_pps.size() < flux_bc_counter + 1)
          paramError("flux_inlet_pps",
                     "More inlet flux BCs than inlet flux pps (" +
                         std::to_string(_flux_inlet_pps.size()) + ")");
 
        if (momentum_inlet_type == "flux-mass")
        {
          params.set<PostprocessorName>("mdot_pp") = _flux_inlet_pps[flux_bc_counter];
          params.set<PostprocessorName>("area_pp") = "area_pp_" + inlet_bdy;
        }
        else
          params.set<PostprocessorName>("velocity_pp") = _flux_inlet_pps[flux_bc_counter];
 
        for (const auto d : make_range(dimension()))
          params.set<MooseFunctorName>(NS::velocity_vector[d]) = _velocity_names[d];
 
        for (const auto d : make_range(dimension()))
        {
          params.set<MooseEnum>("momentum_component") = NS::directions[d];
          params.set<NonlinearVariableName>("variable") = _velocity_names[d];
 
          getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + inlet_bdy, params);
        }
      }
      {
        const std::string bc_type = "WCNSFVMassFluxBC";
        InputParameters params = getFactory().getValidParams(bc_type);
        params.set<MooseFunctorName>(NS::density) = _density_name;
        params.set<NonlinearVariableName>("variable") = _pressure_name;
        params.set<std::vector<BoundaryName>>("boundary") = {inlet_bdy};
 
        if (_flux_inlet_directions.size())
          params.set<Point>("direction") = _flux_inlet_directions[flux_bc_counter];
 
        if (momentum_inlet_type == "flux-mass")
        {
          params.set<PostprocessorName>("mdot_pp") = _flux_inlet_pps[flux_bc_counter];
          params.set<PostprocessorName>("area_pp") = "area_pp_" + inlet_bdy;
        }
        else
          params.set<PostprocessorName>("velocity_pp") = _flux_inlet_pps[flux_bc_counter];
 
        for (const auto d : make_range(dimension()))
          params.set<MooseFunctorName>(NS::velocity_vector[d]) = _velocity_names[d];
 
        getProblem().addFVBC(bc_type, _pressure_name + "_" + inlet_bdy, params);
      }
 
      // need to increment flux_bc_counter
      ++flux_bc_counter;
    }
  }
}

◆ addMassKernels()

void WCNSFVFlowPhysics::addMassKernels ( )

private

Function adding kernels for the incompressible continuity equation.

Definition at line 332 of file WCNSFVFlowPhysics.C.

{
  std::string kernel_type = "INSFVMassAdvection";
  std::string kernel_name = prefix() + "ins_mass_advection";
 
  if (_porous_medium_treatment)
  {
    kernel_type = "PINSFVMassAdvection";
    kernel_name = prefix() + "pins_mass_advection";
  }
 
  InputParameters params = getFactory().getValidParams(kernel_type);
  assignBlocks(params, _blocks);
  params.set<NonlinearVariableName>("variable") = _pressure_name;
  params.set<MooseFunctorName>(NS::density) = _density_name;
  params.set<MooseEnum>("velocity_interp_method") = _velocity_interpolation;
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<MooseEnum>("advected_interp_method") =
      getParam<MooseEnum>("mass_advection_interpolation");
 
  getProblem().addFVKernel(kernel_type, kernel_name, params);
}

Referenced by addFVKernels().

◆ addMassTimeKernels()

void WCNSFVFlowPhysics::addMassTimeKernels ( )

private

Function adding kernels for the time derivative term of the weakly compressible continuity equation.

Definition at line 311 of file WCNSFVFlowPhysics.C.

{
  std::string mass_kernel_type = "WCNSFVMassTimeDerivative";
  std::string kernel_name = prefix() + "wcns_mass_time";
 
  if (_porous_medium_treatment)
  {
    mass_kernel_type = "PWCNSFVMassTimeDerivative";
    kernel_name = prefix() + "pwcns_mass_time";
  }
 
  InputParameters params = getFactory().getValidParams(mass_kernel_type);
  assignBlocks(params, _blocks);
  params.set<NonlinearVariableName>("variable") = _pressure_name;
  params.set<MooseFunctorName>(NS::time_deriv(NS::density)) = NS::time_deriv(_density_name);
  if (_porous_medium_treatment)
    params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
  getProblem().addFVKernel(mass_kernel_type, kernel_name, params);
}

Referenced by addFVKernels().

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

◆ addMomentumAdvectionKernels()

void WCNSFVFlowPhysics::addMomentumAdvectionKernels ( )

private

Definition at line 414 of file WCNSFVFlowPhysics.C.

{
  std::string kernel_type = "INSFVMomentumAdvection";
  std::string kernel_name = prefix() + "ins_momentum_advection_";
 
  if (_porous_medium_treatment)
  {
    kernel_type = "PINSFVMomentumAdvection";
    kernel_name = prefix() + "pins_momentum_advection_";
  }
 
  InputParameters params = getFactory().getValidParams(kernel_type);
  assignBlocks(params, _blocks);
  params.set<MooseFunctorName>(NS::density) = _density_name;
  params.set<MooseEnum>("velocity_interp_method") = _velocity_interpolation;
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<MooseEnum>("advected_interp_method") = _momentum_advection_interpolation;
  if (_porous_medium_treatment)
    params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
  params.applySpecificParameters(parameters(), INSFVMomentumAdvection::listOfCommonParams());
 
  for (const auto d : make_range(dimension()))
  {
    params.set<NonlinearVariableName>("variable") = _velocity_names[d];
    params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
    getProblem().addFVKernel(kernel_type, kernel_name + NS::directions[d], params);
  }
}

Referenced by addFVKernels().

◆ addMomentumBoussinesqKernels()

void WCNSFVFlowPhysics::addMomentumBoussinesqKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 583 of file WCNSFVFlowPhysics.C.

{
  if (_compressibility == "weakly-compressible")
    paramError("boussinesq_approximation",
               "We cannot use boussinesq approximation while running in weakly-compressible mode!");
 
  std::string kernel_type = "INSFVMomentumBoussinesq";
  std::string kernel_name = prefix() + "ins_momentum_boussinesq_";
 
  if (_porous_medium_treatment)
  {
    kernel_type = "PINSFVMomentumBoussinesq";
    kernel_name = prefix() + "pins_momentum_boussinesq_";
  }
 
  InputParameters params = getFactory().getValidParams(kernel_type);
  assignBlocks(params, _blocks);
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<MooseFunctorName>(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");
  if (_porous_medium_treatment)
    params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
  // User declared the flow to be incompressible, we have to trust them
  params.set<bool>("_override_constant_check") = true;
 
  for (const auto d : make_range(dimension()))
  {
    if (getParam<RealVectorValue>("gravity")(d) != 0)
    {
      params.set<MooseEnum>("momentum_component") = NS::directions[d];
      params.set<NonlinearVariableName>("variable") = _velocity_names[d];
 
      getProblem().addFVKernel(kernel_type, kernel_name + NS::directions[d], params);
    }
  }
}

Referenced by addFVKernels().

◆ addMomentumFrictionKernels()

void WCNSFVFlowPhysics::addMomentumFrictionKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 624 of file WCNSFVFlowPhysics.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 = "PINSFVMomentumFriction";
  InputParameters params = getFactory().getValidParams(kernel_type);
  params.set<MooseFunctorName>(NS::density) = _density_name;
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  if (hasForchheimerFriction())
    params.set<MooseFunctorName>(NS::speed) = NS::speed;
  params.set<bool>("standard_friction_formulation") =
      getParam<bool>("standard_friction_formulation");
  params.set<bool>("is_porous_medium") = _porous_medium_treatment;
 
  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<NonlinearVariableName>("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 if (upper_name == "FORCHHEIMER")
        {
          params.set<MooseFunctorName>(NS::speed) = NS::speed;
          params.set<MooseFunctorName>("Forchheimer_name") = _friction_coeffs[block_i][type_i];
        }
        else
          paramError("friction_types",
                     "Friction type '",
                     _friction_types[block_i][type_i],
                     "' is not implemented");
      }
 
      getProblem().addFVKernel(kernel_type,
                               prefix() + "momentum_friction_" + block_name + "_" +
                                   NS::directions[d],
                               params);
    }
 
    if (_porous_medium_treatment && getParam<bool>("use_friction_correction"))
    {
      const std::string correction_kernel_type = "PINSFVMomentumFrictionCorrection";
      InputParameters corr_params = getFactory().getValidParams(correction_kernel_type);
      if (num_friction_blocks)
        corr_params.set<std::vector<SubdomainName>>("block") = _friction_blocks[block_i];
      else
        assignBlocks(corr_params, _blocks);
      corr_params.set<MooseFunctorName>(NS::density) = _density_name;
      corr_params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
      corr_params.set<Real>("consistent_scaling") = getParam<Real>("consistent_scaling");
      for (const auto d : make_range(dimension()))
      {
        corr_params.set<NonlinearVariableName>("variable") = _velocity_names[d];
        corr_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")
          {
            corr_params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
            corr_params.set<MooseFunctorName>("Darcy_name") = _friction_coeffs[block_i][type_i];
          }
          else if (upper_name == "FORCHHEIMER")
          {
            corr_params.set<MooseFunctorName>(NS::speed) = NS::speed;
            corr_params.set<MooseFunctorName>("Forchheimer_name") =
                _friction_coeffs[block_i][type_i];
          }
        }
 
        getProblem().addFVKernel(correction_kernel_type,
                                 prefix() + "pins_momentum_friction_correction_" + block_name +
                                     "_" + NS::directions[d],
                                 corr_params);
      }
    }
  }
}

Referenced by addFVKernels().

◆ addMomentumGravityKernels()

void WCNSFVFlowPhysics::addMomentumGravityKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 548 of file WCNSFVFlowPhysics.C.

{
  if (parameters().isParamValid("gravity") && !_solve_for_dynamic_pressure)
  {
    std::string kernel_type = "INSFVMomentumGravity";
    std::string kernel_name = prefix() + "ins_momentum_gravity_";
 
    if (_porous_medium_treatment)
    {
      kernel_type = "PINSFVMomentumGravity";
      kernel_name = prefix() + "pins_momentum_gravity_";
    }
 
    InputParameters params = getFactory().getValidParams(kernel_type);
    assignBlocks(params, _blocks);
    params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
    params.set<MooseFunctorName>(NS::density) = _density_gravity_name;
    params.set<RealVectorValue>("gravity") = getParam<RealVectorValue>("gravity");
    if (_porous_medium_treatment)
      params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
 
    for (const auto d : make_range(dimension()))
    {
      if (getParam<RealVectorValue>("gravity")(d) != 0)
      {
        params.set<MooseEnum>("momentum_component") = NS::directions[d];
        params.set<NonlinearVariableName>("variable") = _velocity_names[d];
 
        getProblem().addFVKernel(kernel_type, kernel_name + NS::directions[d], params);
      }
    }
  }
}

Referenced by addFVKernels().

◆ addMomentumMixingLengthKernels()

void WCNSFVFlowPhysics::addMomentumMixingLengthKernels ( )

private

◆ addMomentumPressureKernels()

void WCNSFVFlowPhysics::addMomentumPressureKernels ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 519 of file WCNSFVFlowPhysics.C.

{
  std::string kernel_type = "INSFVMomentumPressure";
  std::string kernel_name = prefix() + "ins_momentum_pressure_";
 
  if (_porous_medium_treatment)
  {
    kernel_type = "PINSFVMomentumPressure";
    kernel_name = prefix() + "pins_momentum_pressure_";
  }
 
  InputParameters params = getFactory().getValidParams(kernel_type);
  assignBlocks(params, _blocks);
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<MooseFunctorName>("pressure") = _pressure_name;
  params.set<bool>("correct_skewness") =
      getParam<MooseEnum>("pressure_face_interpolation") == "skewness-corrected";
  if (_porous_medium_treatment)
    params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
 
  for (const auto d : make_range(dimension()))
  {
    params.set<MooseEnum>("momentum_component") = NS::directions[d];
    params.set<NonlinearVariableName>("variable") = _velocity_names[d];
    getProblem().addFVKernel(kernel_type, kernel_name + NS::directions[d], params);
  }
}

Referenced by addFVKernels().

◆ addMomentumTimeKernels()

void WCNSFVFlowPhysics::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 375 of file WCNSFVFlowPhysics.C.

{
  std::string kernel_type = (_compressibility == "weakly-compressible")
                                ? "WCNSFVMomentumTimeDerivative"
                                : "INSFVMomentumTimeDerivative";
  std::string kernel_name = prefix() +
                            ((_compressibility == "weakly-compressible") ? "wcns_" : "ins_") +
                            "momentum_time_";
 
  if (_porous_medium_treatment)
  {
    // Porosity does not appear in the term
    kernel_type = (_compressibility == "weakly-compressible") ? "WCNSFVMomentumTimeDerivative"
                                                              : "PINSFVMomentumTimeDerivative";
    kernel_name = prefix() + ((_compressibility == "weakly-compressible") ? "pwcns_" : "pins_") +
                  "momentum_time_";
  }
 
  InputParameters params = getFactory().getValidParams(kernel_type);
  assignBlocks(params, _blocks);
  params.set<MooseFunctorName>(NS::density) = _density_name;
  if (_compressibility == "weakly-compressible")
    params.set<MooseFunctorName>(NS::time_deriv(NS::density)) = NS::time_deriv(_density_name);
 
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<bool>("contribute_to_rc") =
      getParam<bool>("time_derivative_contributes_to_RC_coefficients");
 
  for (const auto d : make_range(dimension()))
  {
    params.set<NonlinearVariableName>("variable") = _velocity_names[d];
    params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
    if (shouldCreateTimeDerivative(_velocity_names[d], _blocks, /*error if already defined*/ false))
      getProblem().addFVKernel(kernel_type, kernel_name + _velocity_names[d], params);
  }
}

Referenced by addFVKernels().

◆ addMomentumViscousDissipationKernels()

void WCNSFVFlowPhysics::addMomentumViscousDissipationKernels ( )

private

Definition at line 445 of file WCNSFVFlowPhysics.C.

{
  std::string kernel_type = "INSFVMomentumDiffusion";
  std::string kernel_name = prefix() + "ins_momentum_diffusion_";
 
  if (_porous_medium_treatment)
  {
    kernel_type = "PINSFVMomentumDiffusion";
    kernel_name = prefix() + "pins_momentum_diffusion_";
  }
 
  InputParameters params = getFactory().getValidParams(kernel_type);
  assignBlocks(params, _blocks);
  params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
  params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
  const bool user_include_iso = includeIsotropicStress();
  if (user_include_iso && _porous_medium_treatment)
    paramWarning("include_isotropic_viscous_stress",
                 "Including the isotropic viscous stress is not supported with the porous medium "
                 "treatment. Ignoring the request.");
  const bool include_isotropic = (!_porous_medium_treatment) && user_include_iso;
  if (include_isotropic)
    params.set<bool>("include_isotropic_viscous_stress") = true;
  params.set<MooseEnum>("mu_interp_method") = getParam<MooseEnum>("mu_interp_method");
  params.set<MooseEnum>("variable_interp_method") =
      getParam<MooseEnum>("momentum_face_interpolation");
  bool include_symmetric = includeSymmetrizedViscousStress();
  if (include_symmetric && _porous_medium_treatment)
  {
    paramWarning("include_symmetrized_viscous_stress",
                 "Including the symmetrized viscous stress is not supported with the porous "
                 "medium treatment. Ignoring the request.");
    include_symmetric = false;
  }
  if (include_symmetric || include_isotropic)
  {
    params.set<bool>("complete_expansion") = true;
    const std::string u_names[3] = {"u", "v", "w"};
    for (unsigned int i = 0; i < dimension(); ++i)
      params.set<MooseFunctorName>(u_names[i]) = _velocity_names[i];
  }
 
  if (_porous_medium_treatment)
    params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
  // Currently only Newton method for WCNSFVFlowPhysics
  params.set<bool>("newton_solve") = true;
  for (const auto d : make_range(dimension()))
  {
    params.set<NonlinearVariableName>("variable") = _velocity_names[d];
    params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
    getProblem().addFVKernel(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 WCNSFVFlowPhysics::addOutletBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 858 of file WCNSFVFlowPhysics.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)
  {
    if (momentum_outlet_type == "zero-gradient" ||
        momentum_outlet_type == "fixed-pressure-zero-gradient")
    {
      {
        const std::string bc_type = _porous_medium_treatment ? "PINSFVMomentumAdvectionOutflowBC"
                                                             : "INSFVMomentumAdvectionOutflowBC";
        InputParameters params = getFactory().getValidParams(bc_type);
        params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
        if (_porous_medium_treatment)
          params.set<MooseFunctorName>(NS::porosity) = _flow_porosity_functor_name;
        params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
        params.set<MooseFunctorName>(NS::density) = _density_name;
 
        for (unsigned int i = 0; i < dimension(); ++i)
          params.set<MooseFunctorName>(u_names[i]) = _velocity_names[i];
 
        for (const auto d : make_range(dimension()))
        {
          params.set<NonlinearVariableName>("variable") = _velocity_names[d];
          params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
          getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + outlet_bdy, params);
        }
      }
    }
 
    if (momentum_outlet_type == "fixed-pressure" ||
        momentum_outlet_type == "fixed-pressure-zero-gradient")
    {
      const std::string bc_type = "INSFVOutletPressureBC";
      InputParameters params = getFactory().getValidParams(bc_type);
      params.set<NonlinearVariableName>("variable") = _pressure_name;
      params.set<MooseFunctorName>("functor") = libmesh_map_find(_pressure_functors, outlet_bdy);
      params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
 
      getProblem().addFVBC(bc_type, _pressure_name + "_" + outlet_bdy, params);
    }
    else if (momentum_outlet_type == "zero-gradient")
    {
      const std::string bc_type = "INSFVMassAdvectionOutflowBC";
      InputParameters params = getFactory().getValidParams(bc_type);
      params.set<NonlinearVariableName>("variable") = _pressure_name;
      params.set<MooseFunctorName>(NS::density) = _density_name;
      params.set<std::vector<BoundaryName>>("boundary") = {outlet_bdy};
 
      for (const auto d : make_range(dimension()))
        params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
 
      getProblem().addFVBC(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);
    }
}

◆ addPressurePinKernel()

void WCNSFVFlowPhysics::addPressurePinKernel ( )

private

Function adding the pressure constraint.

Definition at line 356 of file WCNSFVFlowPhysics.C.

{
  const auto pin_type = getParam<MooseEnum>("pinned_pressure_type");
  const auto object_type =
      (pin_type == "average") ? "FVIntegralValueConstraint" : "FVPointValueConstraint";
  InputParameters params = getFactory().getValidParams(object_type);
  if (pin_type != "point-value" && pin_type != "average")
    return;
 
  params.set<CoupledName>("lambda") = {"lambda"};
  params.set<PostprocessorName>("phi0") = getParam<PostprocessorName>("pinned_pressure_value");
  params.set<NonlinearVariableName>("variable") = _pressure_name;
  if (pin_type == "point-value")
    params.set<Point>("point") = getParam<Point>("pinned_pressure_point");
 
  getProblem().addFVKernel(object_type, prefix() + "ins_mass_pressure_pin", params);
}

Referenced by addFVKernels().

◆ addRhieChowUserObjects()

void WCNSFVFlowPhysics::addRhieChowUserObjects ( )

overrideprivatevirtual

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

Implements WCNSFVFlowPhysicsBase.

Definition at line 1149 of file WCNSFVFlowPhysics.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);
  bool have_matching_rc_uo = false;
  for (const auto & obj : objs)
    if (const auto * const rc_obj = dynamic_cast<INSFVRhieChowInterpolator *>(obj); rc_obj)
      // Latter check is for whether one of the RC user object is defined everywhere
      if (rc_obj->blocks() == _blocks || (rc_obj->blocks().size() == 0 || _blocks.size() == 0))
      {
        have_matching_rc_uo = true;
        _rc_uo_name = rc_obj->name();
        break;
      }
 
  if (have_matching_rc_uo)
    return;
 
  _rc_uo_name =
      _porous_medium_treatment ? +"pins_rhie_chow_interpolator" : "ins_rhie_chow_interpolator";
 
  const std::string u_names[3] = {"u", "v", "w"};
  const auto object_type =
      _porous_medium_treatment ? "PINSFVRhieChowInterpolator" : "INSFVRhieChowInterpolator";
 
  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;
 
  if (_porous_medium_treatment)
  {
    params.set<MooseFunctorName>(NS::porosity) = _porosity_name;
    unsigned short smoothing_layers = isParamValid("porosity_smoothing_layers")
                                          ? getParam<unsigned short>("porosity_smoothing_layers")
                                          : 0;
    params.set<unsigned short>("smoothing_layers") = smoothing_layers;
  }
 
  if (!_has_flow_equations)
  {
    checkRhieChowFunctorsDefined();
    params.set<MooseFunctorName>("a_u") = "ax";
    params.set<MooseFunctorName>("a_v") = "ay";
    params.set<MooseFunctorName>("a_w") = "az";
  }
 
  params.applySpecificParameters(parameters(), INSFVRhieChowInterpolator::listOfCommonParams());
  addUserObject(object_type, _rc_uo_name, params);
}

Referenced by addUserObjects().

◆ addSeparatorBC()

void WCNSFVFlowPhysics::addSeparatorBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 1049 of file WCNSFVFlowPhysics.C.

{
  if (_hydraulic_separators.size())
  {
    std::string bc_type = "INSFVVelocityHydraulicSeparatorBC";
    InputParameters params = getFactory().getValidParams(bc_type);
    params.set<std::vector<BoundaryName>>("boundary") = _hydraulic_separators;
    params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
 
    for (const auto d : make_range(dimension()))
    {
      params.set<NonlinearVariableName>("variable") = _velocity_names[d];
      params.set<MooseEnum>("momentum_component") = NS::directions[d];
      getProblem().addFVBC(bc_type, prefix() + _velocity_names[d] + "_separators", params);
    }
 
    bc_type = "INSFVScalarFieldSeparatorBC";
    params = getFactory().getValidParams(bc_type);
    params.set<std::vector<BoundaryName>>("boundary") = _hydraulic_separators;
    params.set<NonlinearVariableName>("variable") = _pressure_name;
    getProblem().addFVBC(bc_type, prefix() + _pressure_name + "_separators", params);
  }
}

◆ addSolverVariables()

void WCNSFVFlowPhysics::addSolverVariables ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 141 of file WCNSFVFlowPhysics.C.

{
  if (!_has_flow_equations)
    return;
 
  for (const auto d : make_range(dimension()))
    saveSolverVariableName(_velocity_names[d]);
  saveSolverVariableName(_pressure_name);
 
  // 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",
                                         "momentum_scaling",
                                         "momentum_face_interpolation",
                                         "momentum_two_term_bc_expansion"},
                                        "INSFVVelocityVariable");
    else if (_define_variables)
    {
      std::string variable_type = "INSFVVelocityVariable";
      if (_porous_medium_treatment)
        variable_type = "PINSFVSuperficialVelocityVariable";
 
      auto params = getFactory().getValidParams(variable_type);
      assignBlocks(params, _blocks); // TODO: check wrt components
      params.set<std::vector<Real>>("scaling") = {getParam<Real>("momentum_scaling")};
      params.set<MooseEnum>("face_interp_method") =
          getParam<MooseEnum>("momentum_face_interpolation");
      params.set<bool>("two_term_boundary_expansion") =
          getParam<bool>("momentum_two_term_bc_expansion");
 
      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 WCNSFVFlowPhysics does not exist!");
  }
 
  // Pressure
  const bool using_bernouilli_pressure_var =
      _porous_medium_treatment &&
      getParam<MooseEnum>("porosity_interface_pressure_treatment") != "automatic";
  const auto pressure_type =
      using_bernouilli_pressure_var ? "BernoulliPressureVariable" : "INSFVPressureVariable";
  if (!shouldCreateVariable(_pressure_name, _blocks, /*error if aux*/ true))
  {
    std::vector<std::string> potentially_missed = {"system_names",
                                                   "mass_scaling",
                                                   "pressure_face_interpolation",
                                                   "pressure_two_term_bc_expansion"};
    if (using_bernouilli_pressure_var)
    {
      std::vector<std::string> other_missed = {"pressure_allow_expansion_on_bernoulli_faces",
                                               "pressure_drop_sidesets",
                                               "pressure_drop_form_factors"};
      potentially_missed.insert(potentially_missed.end(), other_missed.begin(), other_missed.end());
    }
    reportPotentiallyMissedParameters(potentially_missed, pressure_type);
  }
  else if (_define_variables)
  {
    auto params = getFactory().getValidParams(pressure_type);
    assignBlocks(params, _blocks);
    params.set<std::vector<Real>>("scaling") = {getParam<Real>("mass_scaling")};
    params.set<MooseEnum>("face_interp_method") =
        getParam<MooseEnum>("pressure_face_interpolation");
    params.set<bool>("two_term_boundary_expansion") =
        getParam<bool>("pressure_two_term_bc_expansion");
 
    if (using_bernouilli_pressure_var)
    {
      params.set<MooseFunctorName>("u") = _velocity_names[0];
      if (dimension() >= 2)
        params.set<MooseFunctorName>("v") = _velocity_names[1];
      if (dimension() == 3)
        params.set<MooseFunctorName>("w") = _velocity_names[2];
      params.set<MooseFunctorName>(NS::porosity) = _porosity_name;
      params.set<MooseFunctorName>(NS::density) = _density_name;
      params.set<bool>("allow_two_term_expansion_on_bernoulli_faces") =
          getParam<bool>("pressure_allow_expansion_on_bernoulli_faces");
      params.set<std::vector<BoundaryName>>("pressure_drop_sidesets") =
          getParam<std::vector<BoundaryName>>("pressure_drop_sidesets");
      params.set<std::vector<Real>>("pressure_drop_form_factors") =
          getParam<std::vector<Real>>("pressure_drop_form_factors");
    }
    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 WCNSFVFlowPhysics does not exist!");
 
  // Add lagrange multiplier for pinning pressure, if needed
  if (getParam<bool>("pin_pressure"))
  {
    auto type = getParam<MooseEnum>("pinned_pressure_type");
    auto lm_params = getFactory().getValidParams("MooseVariableScalar");
    lm_params.set<MooseEnum>("family") = "scalar";
    lm_params.set<MooseEnum>("order") = "first";
 
    if ((type == "point-value" || type == "average"))
    {
      if (!_problem->hasScalarVariable("lambda"))
      {
        lm_params.set<SolverSystemName>("solver_sys") = getSolverSystem("lambda");
        getProblem().addVariable("MooseVariableScalar", "lambda", lm_params);
      }
      else
        reportPotentiallyMissedParameters({"system_names"}, "MooseVariableScalar");
    }
  }
}

◆ addUserObjects()

void WCNSFVFlowPhysics::addUserObjects ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 1074 of file WCNSFVFlowPhysics.C.

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

◆ addWallsBC()

void WCNSFVFlowPhysics::addWallsBC ( )

overrideprivatevirtual

Implements WCNSFVFlowPhysicsBase.

Definition at line 931 of file WCNSFVFlowPhysics.C.

{
  const std::string u_names[3] = {"u", "v", "w"};
 
  // Count the number of fixed velocity wall boundaries (moving walls)
  unsigned int num_functor_walls = 0;
  for (const auto & [boundary_name, momentum_wall_type] : _momentum_wall_types)
    if (momentum_wall_type == "noslip")
      num_functor_walls++;
  if (_momentum_wall_functors.size() && num_functor_walls != _momentum_wall_functors.size())
    paramError("momentum_wall_functors",
               "If any wall functors are specified, the number of boundaries requiring a momentum "
               "functor (" +
                   std::to_string(num_functor_walls) + ") and the number of functors specified (" +
                   std::to_string(_momentum_wall_functors.size()) + ") must match");
  for (const auto & wall_functors : _momentum_wall_functors)
    if (wall_functors.second.size() != dimension())
      paramError("momentum_wall_functors",
                 "Number of wall functors (" + std::to_string(wall_functors.second.size()) +
                     ") must match dimension (" + std::to_string(dimension()) +
                     ").\nFunctors currently specified:" + Moose::stringify(wall_functors.second));
 
  for (const auto & [boundary_name, wall_type] : _momentum_wall_types)
  {
    if (wall_type == "noslip")
    {
      const std::string bc_type = "INSFVNoSlipWallBC";
      InputParameters params = getFactory().getValidParams(bc_type);
      params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
 
      for (const auto d : make_range(dimension()))
      {
        params.set<NonlinearVariableName>("variable") = _velocity_names[d];
        if (_momentum_wall_functors.count(boundary_name) == 0)
          params.set<FunctionName>("function") = "0";
        else
          params.set<FunctionName>("function") = _momentum_wall_functors[boundary_name][d];
 
        getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
      }
    }
    else if (wall_type == "wallfunction")
    {
      const std::string bc_type = "INSFVWallFunctionBC";
      InputParameters params = getFactory().getValidParams(bc_type);
      params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
      params.set<MooseFunctorName>(NS::density) = _density_name;
      params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
      params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
 
      for (const auto d : make_range(dimension()))
        params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
 
      for (const auto d : make_range(dimension()))
      {
        params.set<NonlinearVariableName>("variable") = _velocity_names[d];
        params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
        getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
      }
    }
    else if (wall_type == "slip")
    {
      const std::string bc_type = "INSFVNaturalFreeSlipBC";
      InputParameters params = getFactory().getValidParams(bc_type);
      params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
      params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
 
      for (const auto d : make_range(dimension()))
      {
        params.set<NonlinearVariableName>("variable") = _velocity_names[d];
        params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
        getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
      }
    }
    else if (wall_type == "symmetry")
    {
      {
        std::string bc_type;
        if (_porous_medium_treatment)
          bc_type = "PINSFVSymmetryVelocityBC";
        else
          bc_type = "INSFVSymmetryVelocityBC";
 
        InputParameters params = getFactory().getValidParams(bc_type);
        params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
 
        MooseFunctorName viscosity_name = _dynamic_viscosity_name;
        if (hasTurbulencePhysics())
          viscosity_name = NS::total_viscosity;
        params.set<MooseFunctorName>(NS::mu) = viscosity_name;
        params.set<UserObjectName>("rhie_chow_user_object") = rhieChowUOName();
 
        for (const auto d : make_range(dimension()))
          params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
 
        for (const auto d : make_range(dimension()))
        {
          params.set<NonlinearVariableName>("variable") = _velocity_names[d];
          params.set<MooseEnum>("momentum_component") = NS::directions[d];
 
          getProblem().addFVBC(bc_type, _velocity_names[d] + "_" + boundary_name, params);
        }
      }
      {
        const std::string bc_type = "INSFVSymmetryPressureBC";
        InputParameters params = getFactory().getValidParams(bc_type);
        params.set<NonlinearVariableName>("variable") = _pressure_name;
        params.set<std::vector<BoundaryName>>("boundary") = {boundary_name};
 
        getProblem().addFVBC(bc_type, _pressure_name + "_" + boundary_name, params);
      }
    }
  }
}

◆ checkRhieChowFunctorsDefined()

void WCNSFVFlowPhysics::checkRhieChowFunctorsDefined ( ) const

private

Checks that sufficient Rhie Chow coefficients have been defined for the given dimension, used for scalar or temperature advection by auxiliary variables.

Definition at line 1222 of file WCNSFVFlowPhysics.C.

{
  if (!getProblem().hasFunctor("ax", /*thread_id=*/0))
    mooseError("Rhie Chow coefficient ax must be provided for advection by auxiliary velocities");
  if (dimension() >= 2 && !getProblem().hasFunctor("ay", /*thread_id=*/0))
    mooseError("Rhie Chow coefficient ay must be provided for advection by auxiliary velocities");
  if (dimension() == 3 && !getProblem().hasFunctor("az", /*thread_id=*/0))
    mooseError("Rhie Chow coefficient az must be provided for advection by auxiliary velocities");
}

Referenced by addRhieChowUserObjects().

◆ 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; }

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

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

MooseFunctorName WCNSFVFlowPhysics::getLinearFrictionCoefName ( ) const

overridevirtual

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

Implements WCNSFVFlowPhysicsBase.

Definition at line 1233 of file WCNSFVFlowPhysics.C.

{
  // Check all blocks. If more than one block, they would need to be consolidated #include in
  // a single functor material. We won't implement this for now
  if (_friction_types.empty())
    return "";
  else if (_friction_types.size() == 1)
  {
    for (const auto & type_i : index_range(_friction_types[0]))
    {
      const auto upper_name = MooseUtils::toUpper(_friction_types[0][type_i]);
      if (upper_name == "DARCY")
        return _friction_coeffs[0][type_i];
    }
    // No linear type found
    return "";
  }
  else if (_friction_types.size() > 1)
  {
    bool linear_friction_factor_found = false;
    MooseFunctorName linear_friction_factor;
    for (const auto block_i : index_range(_friction_types))
      for (const auto type_i : index_range(_friction_types[block_i]))
      {
        const auto upper_name = MooseUtils::toUpper(_friction_types[block_i][type_i]);
        if (upper_name == "DARCY" && !linear_friction_factor_found)
        {
          linear_friction_factor_found = true;
          linear_friction_factor = _friction_types[block_i][type_i];
        }
        else if (upper_name == "DARCY" && !linear_friction_factor_found)
          if (linear_friction_factor != _friction_types[block_i][type_i])
            mooseError("Multiple linear friction factor with different names have been specified. "
                       "This is not currently supported as a single name should be retrievable. "
                       "Use a PiecewiseByBlockFunctorMaterial to consolidate them.");
      }
    if (linear_friction_factor_found)
      return linear_friction_factor;
    else
      return "";
  }
  mooseError("Should not get here");
}

◆ 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 WCNSFVFlowPhysics::getNumberAlgebraicGhostingLayersNeeded ( ) const

overridevirtual

Return the number of algebraic ghosting layers needed.

Implements NavierStokesPhysicsBase.

Definition at line 1135 of file WCNSFVFlowPhysics.C.

{
  auto ghost_layers = WCNSFVFlowPhysicsBase::getNumberAlgebraicGhostingLayersNeeded();
  if (_porous_medium_treatment && isParamValid("porosity_smoothing_layers"))
    ghost_layers = std::max(getParam<unsigned short>("porosity_smoothing_layers"), ghost_layers);
  if ((_porous_medium_treatment &&
       getParam<MooseEnum>("porosity_interface_pressure_treatment") != "automatic") ||
      getParam<MooseEnum>("momentum_face_interpolation") == "skewness-corrected" ||
      getParam<MooseEnum>("pressure_face_interpolation") == "skewness-corrected")
    ghost_layers = std::max(ghost_layers, (unsigned short)3);
  return ghost_layers;
}

◆ 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 > WCNSFVFlowPhysics::getSuppliedUserObjects ( ) const

overridevirtual

Reimplemented from PhysicsBase.

Definition at line 1216 of file WCNSFVFlowPhysics.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; }

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

WCNSFVFlowPhysicsBase::_wall_boundaries

const std::vector< BoundaryName > _wall_boundaries

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

Definition WCNSFVFlowPhysicsBase.h:232

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

bool WCNSFVFlowPhysics::hasForchheimerFriction ( ) const

overrideprivatevirtual

Return whether a Forchheimer friction model is in use.

Implements WCNSFVFlowPhysicsBase.

Definition at line 1125 of file WCNSFVFlowPhysics.C.

{
  for (const auto block_i : index_range(_friction_types))
    for (const auto type_i : index_range(_friction_types[block_i]))
      if (MooseUtils::toUpper(_friction_types[block_i][type_i]) == "FORCHHEIMER")
        return true;
  return false;
}

Referenced by addMomentumFrictionKernels().

◆ 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 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 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(), addAxisymmetricViscousSourceKernel(), WCNSLinearFVFlowPhysics::addMomentumFluxKernels(), and addMomentumViscousDissipationKernels().

◆ initializePhysicsAdditional()

void WCNSFVFlowPhysicsBase::initializePhysicsAdditional ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Reimplemented in WCNSLinearFVFlowPhysics.

Definition at line 283 of file WCNSFVFlowPhysicsBase.C.

{
  getProblem().needFV();
}

Referenced by WCNSLinearFVFlowPhysics::initializePhysicsAdditional().

◆ 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;
}

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

◆ usingNavierStokesFVSyntax()

bool NavierStokesPhysicsBase::usingNavierStokesFVSyntax ( ) const

inlineprotectedinherited

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

Definition at line 33 of file NavierStokesPhysicsBase.h.

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

◆ validParams()

InputParameters WCNSFVFlowPhysics::validParams ( )

static

Definition at line 24 of file WCNSFVFlowPhysics.C.

{
  InputParameters params = WCNSFVFlowPhysicsBase::validParams();
  params.addClassDescription(
      "Define the Navier Stokes weakly-compressible mass and momentum equations");
 
  // Rhie Chow interpolation parameters
  params.transferParam<Real>(INSFVMomentumAdvection::validParams(), "characteristic_speed");
  params.addParam<bool>(
      "time_derivative_contributes_to_RC_coefficients",
      true,
      "Whether the time derivative term should contribute to the Rhie Chow coefficients. This adds "
      "stabilization, but makes the solution dependent on the time step size");
  params.addParamNamesToGroup("time_derivative_contributes_to_RC_coefficients characteristic_speed",
                              "Numerical scheme");
 
  // Used for flow mixtures, where one phase is solid / not moving under the action of gravity
  params.addParam<MooseFunctorName>(
      "density_for_gravity_terms",
      "If specified, replaces the 'density' for the Boussinesq and gravity momentum kernels");
 
  // Additional porous media parameters
  params.transferParam<unsigned short>(NSFVBase::validParams(), "porosity_smoothing_layers");
 
  // Techniques to limit or remove oscillations at porosity jump interfaces
  params.transferParam<MooseEnum>(NSFVBase::validParams(), "porosity_interface_pressure_treatment");
  params.transferParam<std::vector<BoundaryName>>(NSFVBase::validParams(),
                                                  "pressure_drop_sidesets");
  params.transferParam<std::vector<Real>>(NSFVBase::validParams(), "pressure_drop_form_factors");
 
  // Friction correction, a technique to limit oscillations at friction interfaces
  params.transferParam<bool>(NSFVBase::validParams(), "use_friction_correction");
  params.transferParam<Real>(NSFVBase::validParams(), "consistent_scaling");
 
  // Couple to turbulence physics
  params.addParam<PhysicsName>("coupled_turbulence_physics",
                               "Turbulence Physics coupled with the flow");
 
  // Spatial discretization scheme
  // Specify the numerical schemes for interpolations of velocity and pressure
  params.transferParam<MooseEnum>(NSFVBase::validParams(), "pressure_face_interpolation");
  params.transferParam<MooseEnum>(NSFVBase::validParams(), "mass_advection_interpolation");
  params.transferParam<bool>(NSFVBase::validParams(),
                             "pressure_allow_expansion_on_bernoulli_faces");
 
  // Nonlinear solver parameters
  params.transferParam<Real>(NSFVBase::validParams(), "mass_scaling");
  params.transferParam<Real>(NSFVBase::validParams(), "momentum_scaling");
 
  // Parameter groups
  params.addParamNamesToGroup("coupled_turbulence_physics", "Coupled Physics");
  params.addParamNamesToGroup(
      "porosity_interface_pressure_treatment pressure_allow_expansion_on_bernoulli_faces "
      "porosity_smoothing_layers use_friction_correction consistent_scaling "
      "pressure_drop_sidesets pressure_drop_form_factors",
      "Flow medium discontinuity treatment");
  params.addParamNamesToGroup(
      "pressure_face_interpolation mass_advection_interpolation mass_scaling momentum_scaling",
      "Numerical scheme");
 
  // TODO Add default preconditioning and move scaling parameters to a preconditioning group
 
  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 addFVKernels(), addMomentumBoussinesqKernels(), WCNSLinearFVFlowPhysics::addMomentumBoussinesqKernels(), 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(), addSolverVariables(), WCNSFVFluidHeatTransferPhysics::addSolverVariables(), WCNSFVTurbulencePhysics::addSolverVariables(), WCNSLinearFVFlowPhysics::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 WCNSLinearFVFlowPhysics::addFunctorMaterials(), addMomentumBoussinesqKernels(), WCNSLinearFVFlowPhysics::addMomentumBoussinesqKernels(), and addMomentumGravityKernels().

◆ _density_name

const MooseFunctorName WCNSFVFlowPhysicsBase::_density_name

protectedinherited

Name of the density material property.

Definition at line 198 of file WCNSFVFlowPhysicsBase.h.

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

◆ _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 addAxisymmetricViscousSourceKernel(), WCNSLinearFVFlowPhysics::addMomentumFluxKernels(), addMomentumFrictionKernels(), WCNSLinearFVFlowPhysics::addMomentumFrictionKernels(), addMomentumViscousDissipationKernels(), 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 addMassTimeKernels(), addMomentumAdvectionKernels(), addMomentumBoussinesqKernels(), addMomentumGravityKernels(), addMomentumPressureKernels(), addMomentumViscousDissipationKernels(), addOutletBC(), WCNSFVFlowPhysicsBase::addPorousMediumSpeedMaterial(), WCNSFVFlowPhysicsBase::getPorosityFunctorName(), and 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(), addMomentumBoussinesqKernels(), WCNSLinearFVFlowPhysics::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 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 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 addMomentumFrictionKernels(), WCNSLinearFVFlowPhysics::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 addMomentumFrictionKernels(), WCNSLinearFVFlowPhysics::addMomentumFrictionKernels(), and 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 addFVKernels(), WCNSLinearFVFlowPhysics::addFVKernels(), addMomentumFrictionKernels(), WCNSLinearFVFlowPhysics::addMomentumFrictionKernels(), getLinearFrictionCoefName(), and 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(), addCorrectors(), WCNSLinearFVFlowPhysics::addFVInterpolationMethods(), addFVKernels(), WCNSLinearFVFlowPhysics::addFVKernels(), addRhieChowUserObjects(), addSolverVariables(), WCNSLinearFVFlowPhysics::addSolverVariables(), WCNSFVFlowPhysicsBase::hasFlowEquations(), and WCNSLinearFVFlowPhysics::WCNSLinearFVFlowPhysics().

◆ _hydraulic_separators

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

protectedinherited

Hydraulic separator boundaries.

Definition at line 234 of file WCNSFVFlowPhysicsBase.h.

Referenced by addSeparatorBC(), WCNSFVFlowPhysicsBase::getHydraulicSeparators(), and WCNSLinearFVFlowPhysics::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 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 addInletBC(), WCNSLinearFVFlowPhysics::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 addInletBC(), WCNSLinearFVFlowPhysics::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 addOutletBC(), WCNSLinearFVFlowPhysics::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 addWallsBC(), WCNSLinearFVFlowPhysics::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 addWallsBC(), WCNSLinearFVFlowPhysics::addWallsBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _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 addInletBC(), addRhieChowUserObjects(), addSolverVariables(), WCNSFVFlowPhysicsBase::getPorosityFunctorName(), and WCNSFVFlowPhysics().

◆ _porosity_smoothing_layers

const unsigned WCNSFVFlowPhysics::_porosity_smoothing_layers

private

The number of smoothing layers if that treatment is used on porosity.

Definition at line 76 of file WCNSFVFlowPhysics.h.

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

Referenced by WCNSFVFlowPhysicsBase::addFluidPropertiesFunctorMaterial(), addInletBC(), addMassKernels(), addMassTimeKernels(), WCNSFVFlowPhysicsBase::addMaterials(), addMomentumAdvectionKernels(), addMomentumBoussinesqKernels(), addMomentumFrictionKernels(), addMomentumGravityKernels(), addMomentumPressureKernels(), addMomentumTimeKernels(), addMomentumViscousDissipationKernels(), addOutletBC(), WCNSFVFlowPhysicsBase::addPorousMediumSpeedMaterial(), addRhieChowUserObjects(), addSolverVariables(), WCNSLinearFVFlowPhysics::addUserObjects(), addWallsBC(), getNumberAlgebraicGhostingLayersNeeded(), WCNSFVFlowPhysicsBase::porousMediumTreatment(), WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase(), and WCNSLinearFVFlowPhysics::WCNSLinearFVFlowPhysics().

◆ _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 addOutletBC(), WCNSLinearFVFlowPhysics::addOutletBC(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

◆ _pressure_name

const NonlinearVariableName WCNSFVFlowPhysicsBase::_pressure_name

protectedinherited

◆ _pressure_tag

const TagName WCNSFVFlowPhysicsBase::_pressure_tag = "p_tag"

protectedinherited

Name of the vector to hold pressure momentum equation contributions.

Definition at line 171 of file WCNSFVFlowPhysicsBase.h.

◆ _rc_uo_name

UserObjectName WCNSFVFlowPhysics::_rc_uo_name

private

Name of the user object in charge of computing the Rhie Chow coefficients.

Definition at line 79 of file WCNSFVFlowPhysics.h.

Referenced by addRhieChowUserObjects().

◆ _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(), addMomentumGravityKernels(), WCNSLinearFVFlowPhysics::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(), WCNSLinearFVFlowPhysics::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 addMassKernels(), 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 WCNSLinearFVFlowPhysics::addWallsBC(), WCNSFVFlowPhysicsBase::getWallBoundaries(), and WCNSFVFlowPhysicsBase::WCNSFVFlowPhysicsBase().

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

navier_stokes/include/physics/WCNSFVFlowPhysics.h
navier_stokes/src/physics/WCNSFVFlowPhysics.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

◆ WCNSFVFlowPhysics()

Member Function Documentation

◆ actOnAdditionalTasks()

◆ addAxisymmetricViscousSource()

◆ addAxisymmetricViscousSourceEnabled()

◆ addAxisymmetricViscousSourceKernel()

◆ addCorrectors()

◆ addFluidPropertiesFunctorMaterial()

◆ addFVAdvectedInterpolationMethod()

◆ addFVBCs()

◆ addFVKernels()

◆ addInitialConditions()

◆ addInletBC()

◆ addMassKernels()

◆ addMassTimeKernels()

◆ addMaterials()

◆ addMomentumAdvectionKernels()

◆ addMomentumBoussinesqKernels()

◆ addMomentumFrictionKernels()

◆ addMomentumGravityKernels()

◆ addMomentumMixingLengthKernels()

◆ addMomentumPressureKernels()

◆ addMomentumTimeKernels()

◆ addMomentumViscousDissipationKernels()

◆ addNonPorousMediumSpeedMaterial()

◆ addOutletBC()

◆ addPorousMediumSpeedMaterial()

◆ addPostprocessors()

◆ addPressurePinKernel()

◆ addRhieChowUserObjects()

◆ addSeparatorBC()

◆ addSolverVariables()

◆ addUserObjects()

◆ addWallsBC()

◆ checkRhieChowFunctorsDefined()

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