Creates all the objects needed to add a turbulence model to an incompressible / weakly-compressible Navier Stokes finite volume flow simulation. More...

#include <WCNSFVTurbulencePhysics.h>

Inheritance diagram for WCNSFVTurbulencePhysics:

[legend]

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	WCNSFVTurbulencePhysics (const InputParameters &parameters)

bool	hasTurbulenceModel () const
	Whether a turbulence model is in use. More...

std::vector< BoundaryName >	turbulenceWalls () const
	The names of the boundaries with turbulence wall functions. More...

MooseEnum	turbulenceEpsilonWallTreatment () const
	The turbulence epsilon wall treatment (same for all turbulence walls currently) More...

MooseEnum	turbulenceTemperatureWallTreatment () const
	The turbulence temperature wall treatment (same for all turbulence walls currently) More...

MooseFunctorName	tkeName () const
	The name of the turbulent kinetic energy variable. More...

virtual void	act () override final

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

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

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

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

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

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

unsigned int	dimension () const

const ActionComponent &	getActionComponent (const ComponentName &comp_name) const

void	checkComponentType (const ActionComponent &component) const

virtual void	addComponent (const ActionComponent &component)

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

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

void	timedAct ()

MooseObjectName	uniqueActionName () const

const std::string &	specificTaskName () const

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

void	appendTask (const std::string &task)

MooseApp &	getMooseApp () const

const std::string &	type () const

virtual const std::string &	name () const

std::string	typeAndName () const

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

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

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

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

PerfGraph &	perfGraph ()

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

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

const WCNSFVFlowPhysicsBase *	getCoupledFlowPhysics () const

const WCNSFVTurbulencePhysics *	getCoupledTurbulencePhysics () const

MooseFunctorName	getPorosityFunctorName (bool smoothed) const
	Return the porosity functor name. More...

const MooseFunctorName &	densityName () const

const MooseFunctorName &	dynamicViscosityName () const

Static Public Member Functions
static InputParameters	validParams ()

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static constexpr auto	SYSTEM

static constexpr auto	NAME

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

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

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

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

bool	isTransient () const

Factory &	getFactory ()

Factory &	getFactory () const

virtual FEProblemBase &	getProblem ()

virtual const FEProblemBase &	getProblem () const

void	prepareCopyVariablesFromMesh () const

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

std::string	prefix () const

void	saveSolverVariableName (const VariableName &var_name)

void	saveAuxVariableName (const VariableName &var_name)

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

bool	solverVariableExists (const VariableName &var_name) const

const SolverSystemName &	getSolverSystem (unsigned int variable_index) const

const SolverSystemName &	getSolverSystem (const VariableName &variable_name) const

void	addRequiredPhysicsTask (const std::string &task)

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

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

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

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

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

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

bool	isVariableFV (const VariableName &var_name) const

bool	isVariableScalar (const VariableName &var_name) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

std::vector< SolverSystemName >	_system_names

std::vector< unsigned int >	_system_numbers

const bool	_verbose

const MooseEnum &	_preconditioning

std::vector< SubdomainName >	_blocks

std::string	_registered_identifier

std::string	_specific_task_name

std::set< std::string >	_all_tasks

ActionWarehouse &	_awh

const std::string &	_current_task

std::shared_ptr< MooseMesh > &	_mesh

std::shared_ptr< MooseMesh > &	_displaced_mesh

std::shared_ptr< FEProblemBase > &	_problem

PerfID	_act_timer

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

const NavierStokesPhysicsBase *	_advection_physics
	The Physics class using this helper. More...

const WCNSFVFlowPhysicsBase *	_flow_equations_physics
	Flow physics. More...

const WCNSFVTurbulencePhysics *	_turbulence_physics
	Turbulence. More...

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

const bool	_porous_medium_treatment
	Switch to show if porous medium treatment is requested or not. More...

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

const NonlinearVariableName	_pressure_name
	Pressure name. More...

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

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

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

Private Member Functions
virtual void	initializePhysicsAdditional () override

virtual void	actOnAdditionalTasks () override

void	retrieveCoupledPhysics ()
	Retrieve the other WCNSFVPhysics at play in the simulation to be able to add the relevant terms (turbulent diffusion notably) More...

virtual void	addSolverVariables () override

virtual void	addAuxiliaryVariables () override

virtual void	addFVKernels () override

virtual void	addFVBCs () override

virtual void	addInitialConditions () override

virtual void	addAuxiliaryKernels () override

virtual void	addMaterials () override

void	addFlowTurbulenceKernels ()
	Functions adding kernels for turbulence in the other equation(s) More...

void	addFluidEnergyTurbulenceKernels ()

void	addScalarAdvectionTurbulenceKernels ()

void	addKEpsilonTimeDerivatives ()
	Functions adding kernels for the k-epsilon to the k-epsilon equations. More...

void	addKEpsilonAdvection ()

void	addKEpsilonDiffusion ()

void	addKEpsilonSink ()

Private Attributes
const MooseEnum	_turbulence_model
	Turbulence model to create the equation(s) for. More...

bool	_has_flow_equations

bool	_has_energy_equation

bool	_has_scalar_equations

const WCNSFVFluidHeatTransferPhysics *	_fluid_energy_physics
	The heat advection physics to add turbulent mixing for. More...

const WCNSFVScalarTransportPhysics *	_scalar_transport_physics
	The scalar advection physics to add turbulent mixing for. More...

const VariableName	_mixing_length_name
	Name of the mixing length auxiliary variable. More...

std::vector< BoundaryName >	_turbulence_walls
	List of boundaries to act as walls for turbulence models. More...

MooseEnum	_wall_treatment_eps
	Turbulence wall treatment for epsilon (same for all walls currently) More...

MooseEnum	_wall_treatment_temp
	Turbulence wall treatment for temperature (same for all walls currently) More...

const VariableName	_tke_name
	Name of the turbulent kinetic energy. More...

const VariableName	_tked_name
	Name of the turbulent kinetic energy dissipation. More...

const VariableName	_turbulent_viscosity_name = NS::mu_t
	Name of the turbulence viscosity auxiliary variable (or property) More...

Detailed Description

Creates all the objects needed to add a turbulence model to an incompressible / weakly-compressible Navier Stokes finite volume flow simulation.

Definition at line 23 of file WCNSFVTurbulencePhysics.h.

Constructor & Destructor Documentation

◆ WCNSFVTurbulencePhysics()

WCNSFVTurbulencePhysics::WCNSFVTurbulencePhysics ( const InputParameters & parameters )

Definition at line 187 of file WCNSFVTurbulencePhysics.C.

   : NavierStokesPhysicsBase(parameters),
     WCNSFVCoupledAdvectionPhysicsHelper(this),
     _turbulence_model(getParam<MooseEnum>("turbulence_handling")),
     _mixing_length_name(getParam<AuxVariableName>("mixing_length_name")),
     _turbulence_walls(getParam<std::vector<BoundaryName>>("turbulence_walls")),
     _wall_treatment_eps(getParam<MooseEnum>("wall_treatment_eps")),
     _wall_treatment_temp(getParam<MooseEnum>("wall_treatment_T")),
     _tke_name(getParam<MooseFunctorName>("tke_name")),
     _tked_name(getParam<MooseFunctorName>("tked_name"))
 {
   if (_verbose && _turbulence_model != "none")
     _console << "Creating a " << std::string(_turbulence_model) << " turbulence model."
              << std::endl;
 
   // Keep track of the variable names, for loading variables from files notably
   if (_turbulence_model == "mixing-length")
     saveAuxVariableName(_mixing_length_name);
   else if (_turbulence_model == "k-epsilon")
   {
     saveSolverVariableName(_tke_name);
     saveSolverVariableName(_tked_name);
     if (getParam<bool>("mu_t_as_aux_variable"))
       saveAuxVariableName(_turbulent_viscosity_name);
     if (getParam<bool>("k_t_as_aux_variable"))
       saveAuxVariableName(NS::k_t);
   }
 
   // Parameter checks
   if (_turbulence_model == "none")
     errorInconsistentDependentParameter("turbulence_handling", "none", {"turbulence_walls"});
   if (_turbulence_model != "mixing-length")
     errorDependentParameter("turbulence_handling",
                             "mixing-length",
                             {"mixing_length_delta",
                              "mixing_length_aux_execute_on",
                              "von_karman_const",
                              "von_karman_const_0",
                              "mixing_length_two_term_bc_expansion"});
   if (_turbulence_model != "k-epsilon")
   {
     errorDependentParameter("turbulence_handling",
                             "k-epsilon",
                             {"C_mu",
                              "C1_eps",
                              "C2_eps",
                              "bulk_wall_treatment",
                              "tke_scaling",
                              "tke_face_interpolation",
                              "tke_two_term_bc_expansion",
                              "tked_scaling",
                              "tked_face_interpolation",
                              "tked_two_term_bc_expansion",
                              "turbulent_viscosity_two_term_bc_expansion"});
     checkSecondParamSetOnlyIfFirstOneTrue("mu_t_as_aux_variable", "initial_mu_t");
   }
 }

Member Function Documentation

◆ actOnAdditionalTasks()

void WCNSFVTurbulencePhysics::actOnAdditionalTasks ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 253 of file WCNSFVTurbulencePhysics.C.

 {
   // Other Physics may not exist or be initialized at construction time, so
   // we retrieve them now, on this task which occurs after 'init_physics'
   if (_current_task == "get_turbulence_physics")
     retrieveCoupledPhysics();
 }

◆ addAuxiliaryKernels()

void WCNSFVTurbulencePhysics::addAuxiliaryKernels ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 730 of file WCNSFVTurbulencePhysics.C.

 {
   // Note that if we are restarting this will overwrite the restarted mixing-length
   if (_turbulence_model == "mixing-length")
   {
     const std::string ml_kernel_type = "WallDistanceMixingLengthAux";
     InputParameters ml_params = getFactory().getValidParams(ml_kernel_type);
     assignBlocks(ml_params, _blocks);
     ml_params.set<AuxVariableName>("variable") = _mixing_length_name;
     ml_params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
     if (parameters().isParamValid("mixing_length_aux_execute_on"))
       ml_params.set<ExecFlagEnum>("execute_on") =
           getParam<ExecFlagEnum>("mixing_length_aux_execute_on");
     else
       ml_params.set<ExecFlagEnum>("execute_on") = {EXEC_INITIAL, EXEC_TIMESTEP_END};
     ml_params.set<MooseFunctorName>("von_karman_const") =
         getParam<MooseFunctorName>("von_karman_const");
     ml_params.set<MooseFunctorName>("von_karman_const_0") =
         getParam<MooseFunctorName>("von_karman_const_0");
     ml_params.set<MooseFunctorName>("delta") = getParam<MooseFunctorName>("mixing_length_delta");
 
     getProblem().addAuxKernel(ml_kernel_type, prefix() + "mixing_length_aux ", ml_params);
   }
 
   if (_turbulence_model == "k-epsilon" && getParam<bool>("mu_t_as_aux_variable"))
   {
     const std::string u_names[3] = {"u", "v", "w"};
     const std::string mut_kernel_type = "kEpsilonViscosityAux";
     InputParameters params = getFactory().getValidParams(mut_kernel_type);
     assignBlocks(params, _blocks);
     params.set<AuxVariableName>("variable") = _turbulent_viscosity_name;
     for (const auto d : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
     params.set<MooseFunctorName>(NS::TKE) = _tke_name;
     params.set<MooseFunctorName>(NS::TKED) = _tked_name;
     params.set<MooseFunctorName>(NS::density) = _flow_equations_physics->densityName();
     params.set<MooseFunctorName>(NS::mu) = _flow_equations_physics->dynamicViscosityName();
     params.set<Real>("C_mu") = getParam<Real>("C_mu");
     params.set<MooseEnum>("wall_treatment") = _wall_treatment_eps;
     params.set<bool>("bulk_wall_treatment") = getParam<bool>("bulk_wall_treatment");
     params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
     params.set<ExecFlagEnum>("execute_on") = {EXEC_NONLINEAR};
     params.set<bool>("newton_solve") = true;
     getProblem().addAuxKernel(mut_kernel_type, prefix() + "mixing_length_aux ", params);
   }
   if (_turbulence_model == "k-epsilon" && getParam<bool>("k_t_as_aux_variable") &&
       _has_energy_equation)
   {
     const std::string kt_kernel_type = "TurbulentConductivityAux";
     InputParameters params = getFactory().getValidParams(kt_kernel_type);
     assignBlocks(params, _blocks);
     params.set<AuxVariableName>("variable") = NS::k_t;
     params.set<MooseFunctorName>(NS::cp) = _fluid_energy_physics->getSpecificHeatName();
     params.set<MooseFunctorName>(NS::mu_t) = _turbulent_viscosity_name;
     params.set<MooseFunctorName>(NS::turbulent_Prandtl) = getParam<MooseFunctorName>("Pr_t");
     params.set<ExecFlagEnum>("execute_on") = {EXEC_INITIAL, EXEC_NONLINEAR};
     getProblem().addAuxKernel(kt_kernel_type, prefix() + "turbulent_conductivity_aux ", params);
   }
 }

◆ addAuxiliaryVariables()

void WCNSFVTurbulencePhysics::addAuxiliaryVariables ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 385 of file WCNSFVTurbulencePhysics.C.

 {
   if (_turbulence_model == "mixing-length" && _define_variables)
   {
     auto params = getFactory().getValidParams("MooseVariableFVReal");
     assignBlocks(params, _blocks);
     if (isParamValid("mixing_length_two_term_bc_expansion"))
       params.set<bool>("two_term_boundary_expansion") =
           getParam<bool>("mixing_length_two_term_bc_expansion");
     if (!shouldCreateVariable(_tke_name, _blocks, /*error if aux*/ false))
       reportPotentiallyMissedParameters({"mixing_length_two_term_bc_expansion"},
                                         "MooseVariableFVReal");
     else
       getProblem().addAuxVariable("MooseVariableFVReal", _mixing_length_name, params);
   }
   if (_turbulence_model == "k-epsilon" && getParam<bool>("mu_t_as_aux_variable"))
   {
     auto params = getFactory().getValidParams("MooseVariableFVReal");
     assignBlocks(params, _blocks);
     if (isParamValid("turbulent_viscosity_two_term_bc_expansion"))
       params.set<bool>("two_term_boundary_expansion") =
           getParam<bool>("turbulent_viscosity_two_term_bc_expansion");
     if (!shouldCreateVariable(_turbulent_viscosity_name, _blocks, /*error if aux*/ false))
       reportPotentiallyMissedParameters({"turbulent_viscosity_two_term_bc_expansion"},
                                         "MooseVariableFVReal");
     else
       getProblem().addAuxVariable("MooseVariableFVReal", _turbulent_viscosity_name, params);
   }
   if (_turbulence_model == "k-epsilon" && getParam<bool>("k_t_as_aux_variable"))
   {
     auto params = getFactory().getValidParams("MooseVariableFVReal");
     assignBlocks(params, _blocks);
     if (shouldCreateVariable(NS::k_t, _blocks, /*error if aux*/ false))
       getProblem().addAuxVariable("MooseVariableFVReal", NS::k_t, params);
   }
 }

◆ addFlowTurbulenceKernels()

void WCNSFVTurbulencePhysics::addFlowTurbulenceKernels ( )

private

Functions adding kernels for turbulence in the other equation(s)

Definition at line 448 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   if (_turbulence_model == "mixing-length")
   {
     const std::string u_names[3] = {"u", "v", "w"};
     const std::string kernel_type = "INSFVMixingLengthReynoldsStress";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<MooseFunctorName>(NS::density) = _density_name;
     params.set<MooseFunctorName>(NS::mixing_length) = _mixing_length_name;
 
     std::string kernel_name = prefix() + "ins_momentum_mixing_length_reynolds_stress_";
     if (_porous_medium_treatment)
       kernel_name = prefix() + "pins_momentum_mixing_length_reynolds_stress_";
 
     params.set<UserObjectName>("rhie_chow_user_object") = _flow_equations_physics->rhieChowUOName();
     for (const auto dim_i : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[dim_i]) = _velocity_names[dim_i];
 
     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);
     }
   }
   else if (_turbulence_model == "k-epsilon")
   {
     // We rely on using the turbulent viscosity in the flow equation
     // This check is rudimentary, we should think of a better way
     // We could also check for the use of 'mu_t' with the right parameters already
     if (_flow_equations_physics->dynamicViscosityName() != "mu" &&
         !MooseUtils::isFloat(_flow_equations_physics->dynamicViscosityName()))
       mooseError(
           "Regular fluid viscosity 'mu' should be used for the momentum diffusion term. You are "
           "currently using: " +
           _flow_equations_physics->dynamicViscosityName());
 
     const std::string u_names[3] = {"u", "v", "w"};
     const std::string kernel_type = "INSFVMomentumDiffusion";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<MooseFunctorName>("mu") = _turbulent_viscosity_name;
     params.set<MooseEnum>("mu_interp_method") =
         getParam<MooseEnum>("turbulent_viscosity_interp_method");
     params.set<MooseEnum>("variable_interp_method") =
         _flow_equations_physics->getMomentumFaceInterpolationMethod();
     params.set<bool>("complete_expansion") = true;
 
     std::string kernel_name = prefix() + "ins_momentum_k_epsilon_reynolds_stress_";
     if (_porous_medium_treatment)
       kernel_name = prefix() + "pins_momentum_k_epsilon_reynolds_stress_";
 
     params.set<UserObjectName>("rhie_chow_user_object") = _flow_equations_physics->rhieChowUOName();
     for (const auto dim_i : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[dim_i]) = _velocity_names[dim_i];
 
     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);
     }
   }
 }

◆ addFluidEnergyTurbulenceKernels()

void WCNSFVTurbulencePhysics::addFluidEnergyTurbulenceKernels ( )

private

Definition at line 517 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   if (_turbulence_model == "mixing-length")
   {
     const std::string u_names[3] = {"u", "v", "w"};
     const std::string kernel_type = "WCNSFVMixingLengthEnergyDiffusion";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<MooseFunctorName>(NS::density) = _density_name;
     params.set<MooseFunctorName>(NS::cp) = _fluid_energy_physics->getSpecificHeatName();
     params.set<MooseFunctorName>(NS::mixing_length) = _mixing_length_name;
     params.set<Real>("schmidt_number") = getParam<Real>("turbulent_prandtl");
     params.set<NonlinearVariableName>("variable") =
         _flow_equations_physics->getFluidTemperatureName();
 
     for (const auto dim_i : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[dim_i]) = _velocity_names[dim_i];
 
     if (_porous_medium_treatment)
       getProblem().addFVKernel(
           kernel_type, prefix() + "pins_energy_mixing_length_diffusion", params);
     else
       getProblem().addFVKernel(
           kernel_type, prefix() + "ins_energy_mixing_length_diffusion", params);
   }
   else if (_turbulence_model == "k-epsilon")
   {
     const std::string kernel_type = "FVDiffusion";
     const auto T_fluid_name = _flow_equations_physics->getFluidTemperatureName();
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<NonlinearVariableName>("variable") = T_fluid_name;
     params.set<MooseFunctorName>("coeff") = NS::k_t;
     getProblem().addFVKernel(kernel_type, prefix() + T_fluid_name + "_turbulent_diffusion", params);
   }
 }

◆ addFVBCs()

void WCNSFVTurbulencePhysics::addFVBCs ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 791 of file WCNSFVTurbulencePhysics.C.

 {
   const std::string u_names[3] = {"u", "v", "w"};
 
   if (_turbulence_model == "k-epsilon" && getParam<bool>("mu_t_as_aux_variable"))
   {
     mooseAssert(_flow_equations_physics, "Should have a flow equation physics");
     const std::string bc_type = "INSFVTurbulentViscosityWallFunction";
     InputParameters params = getFactory().getValidParams(bc_type);
     params.set<std::vector<BoundaryName>>("boundary") = _turbulence_walls;
     params.set<NonlinearVariableName>("variable") = _turbulent_viscosity_name;
     params.set<MooseFunctorName>(NS::density) = _flow_equations_physics->densityName();
     params.set<MooseFunctorName>(NS::mu) = _flow_equations_physics->dynamicViscosityName();
     params.set<MooseFunctorName>(NS::mu_t) = _turbulent_viscosity_name;
     params.set<MooseFunctorName>(NS::TKE) = _tke_name;
     params.set<MooseEnum>("wall_treatment") = _wall_treatment_eps;
     for (const auto d : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
 
     getProblem().addFVBC(bc_type, prefix() + "turbulence_walls", params);
     // Energy wall function boundary conditions are added in the WCNSFVFluidEnergyPhysics
     // because it facilitates counting the number of walls, specifying energy wall functors
     // the same way as for boundary conditions
   }
 }

◆ addFVKernels()

void WCNSFVTurbulencePhysics::addFVKernels ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 423 of file WCNSFVTurbulencePhysics.C.

 {
   if (_turbulence_model == "none")
     return;
 
   // Turbulence terms in other equations
   if (_has_flow_equations)
     addFlowTurbulenceKernels();
   if (_has_energy_equation)
     addFluidEnergyTurbulenceKernels();
   if (_has_scalar_equations)
     addScalarAdvectionTurbulenceKernels();
 
   // Turbulence models with their own set of equations
   if (_turbulence_model == "k-epsilon")
   {
     if (isTransient())
       addKEpsilonTimeDerivatives();
     addKEpsilonAdvection();
     addKEpsilonDiffusion();
     addKEpsilonSink();
   }
 }

◆ addInitialConditions()

void WCNSFVTurbulencePhysics::addInitialConditions ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 818 of file WCNSFVTurbulencePhysics.C.

 {
   if (_turbulence_model == "mixing-length" || _turbulence_model == "none")
     return;
   const std::string ic_type = "FunctionIC";
   InputParameters params = getFactory().getValidParams(ic_type);
 
   // Parameter checking: error if initial conditions are provided but not going to be used
   if ((getParam<bool>("initialize_variables_from_mesh_file") || !_define_variables) &&
       ((getParam<bool>("mu_t_as_aux_variable") && isParamValid("initial_mu_t")) ||
        isParamSetByUser("initial_tke") || isParamSetByUser("initial_tked")))
     mooseError("inital_mu_t/tke/tked should not be provided if we are restarting from a mesh file "
                "or not defining variables in the Physics");
 
   // do not set initial conditions if we are not defining variables
   if (!_define_variables)
     return;
   // on regular restarts (from checkpoint), we obey the user specification of initial conditions
 
   if (getParam<bool>("mu_t_as_aux_variable"))
   {
     const auto rho_name = _flow_equations_physics->densityName();
     // If the user provided an initial value, we use that
     if (isParamValid("initial_mu_t"))
       params.set<FunctionName>("function") = getParam<FunctionName>("initial_mu_t");
     // If we can compute the initialization value from the user parameters, we do that
     else if (MooseUtils::isFloat(rho_name) &&
              MooseUtils::isFloat(getParam<FunctionName>("initial_tke")) &&
              MooseUtils::isFloat(getParam<FunctionName>("initial_tked")))
       params.set<FunctionName>("function") =
           std::to_string(std::atof(rho_name.c_str()) * getParam<Real>("C_mu") *
                          std::pow(std::atof(getParam<FunctionName>("initial_tke").c_str()), 2) /
                          std::atof(getParam<FunctionName>("initial_tked").c_str()));
     else
       paramError("initial_mu_t",
                  "Initial turbulent viscosity should be provided. A sensible value is "
                  "rho * C_mu TKE_initial^2 / TKED_initial");
 
     params.set<VariableName>("variable") = _turbulent_viscosity_name;
     // Always obey the user specification of an initial condition
     if (shouldCreateIC(_turbulent_viscosity_name,
                        _blocks,
                        /*whether IC is a default*/ !isParamSetByUser("initial_mu_t"),
                        /*error if already an IC*/ isParamSetByUser("initial_mu_t")))
       getProblem().addInitialCondition(ic_type, prefix() + "initial_mu_turb", params);
   }
   else if (isParamSetByUser("initial_mu_t"))
     paramError("initial_mu_t",
                "This parameter can only be specified if 'mu_t_as_aux_variable=true'");
 
   params.set<VariableName>("variable") = _tke_name;
   params.set<FunctionName>("function") = getParam<FunctionName>("initial_tke");
   if (shouldCreateIC(_tke_name,
                      _blocks,
                      /*whether IC is a default*/ !isParamSetByUser("initial_tke"),
                      /*error if already an IC*/ isParamSetByUser("initial_tke")))
     getProblem().addInitialCondition(ic_type, prefix() + "initial_tke", params);
   params.set<VariableName>("variable") = _tked_name;
   params.set<FunctionName>("function") = getParam<FunctionName>("initial_tked");
   if (shouldCreateIC(_tked_name,
                      _blocks,
                      /*whether IC is a default*/ !isParamSetByUser("initial_tked"),
                      /*error if already an IC*/ isParamSetByUser("initial_tked")))
     getProblem().addInitialCondition(ic_type, prefix() + "initial_tked", params);
 }

◆ addKEpsilonAdvection()

void WCNSFVTurbulencePhysics::addKEpsilonAdvection ( )

private

Definition at line 622 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   const std::string kernel_type = "INSFVTurbulentAdvection";
   InputParameters params = getFactory().getValidParams(kernel_type);
 
   assignBlocks(params, _blocks);
 
   params.set<MooseEnum>("velocity_interp_method") = _velocity_interpolation;
   params.set<UserObjectName>("rhie_chow_user_object") = _flow_equations_physics->rhieChowUOName();
   params.set<MooseFunctorName>(NS::density) = _flow_equations_physics->densityName();
   params.set<bool>("neglect_advection_derivatives") =
       getParam<bool>("neglect_advection_derivatives");
 
   params.set<MooseEnum>("advected_interp_method") =
       getParam<MooseEnum>("tke_advection_interpolation");
   params.set<NonlinearVariableName>("variable") = _tke_name;
   getProblem().addFVKernel(kernel_type, prefix() + "tke_advection", params);
   params.set<NonlinearVariableName>("variable") = _tked_name;
   params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
   params.set<MooseEnum>("advected_interp_method") =
       getParam<MooseEnum>("tked_advection_interpolation");
   getProblem().addFVKernel(kernel_type, prefix() + "tked_advection", params);
 }

◆ addKEpsilonDiffusion()

void WCNSFVTurbulencePhysics::addKEpsilonDiffusion ( )

private

Definition at line 647 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   {
     const std::string kernel_type = "INSFVTurbulentDiffusion";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
 
     params.set<NonlinearVariableName>("variable") = _tke_name;
     params.set<MooseFunctorName>("coeff") = _flow_equations_physics->dynamicViscosityName();
     getProblem().addFVKernel(kernel_type, prefix() + "tke_diffusion_mu", params);
 
     params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
     params.set<NonlinearVariableName>("variable") = _tked_name;
     getProblem().addFVKernel(kernel_type, prefix() + "tked_diffusion_mu", params);
   }
 
   {
     const std::string kernel_type = "INSFVTurbulentDiffusion";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
 
     params.set<NonlinearVariableName>("variable") = _tke_name;
     params.set<MooseFunctorName>("coeff") = _turbulent_viscosity_name;
     params.set<MooseFunctorName>("scaling_coef") = getParam<MooseFunctorName>("sigma_k");
     params.set<MooseEnum>("coeff_interp_method") =
         getParam<MooseEnum>("turbulent_viscosity_interp_method");
     getProblem().addFVKernel(kernel_type, prefix() + "tke_diffusion_mu_turb", params);
 
     params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
     params.set<NonlinearVariableName>("variable") = _tked_name;
     params.set<MooseFunctorName>("scaling_coef") = getParam<MooseFunctorName>("sigma_eps");
     getProblem().addFVKernel(kernel_type, prefix() + "tked_diffusion_mu_turb", params);
   }
 }

◆ addKEpsilonSink()

void WCNSFVTurbulencePhysics::addKEpsilonSink ( )

private

Definition at line 683 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   const std::string u_names[3] = {"u", "v", "w"};
   {
     const std::string kernel_type = "INSFVTKESourceSink";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<NonlinearVariableName>("variable") = _tke_name;
     params.set<MooseFunctorName>(NS::TKED) = _tked_name;
     params.set<MooseFunctorName>(NS::density) = _flow_equations_physics->densityName();
     params.set<MooseFunctorName>(NS::mu) = _flow_equations_physics->dynamicViscosityName();
     params.set<MooseFunctorName>(NS::mu_t) = _turbulent_viscosity_name;
     params.set<Real>("C_pl") = getParam<Real>("C_pl");
     params.set<bool>("linearized_model") = getParam<bool>("linearize_sink_sources");
     params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
     params.set<MooseEnum>("wall_treatment") = _wall_treatment_eps;
     // Currently only Newton method for WCNSFVTurbulencePhysics
     params.set<bool>("newton_solve") = true;
     for (const auto d : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
     getProblem().addFVKernel(kernel_type, prefix() + "tke_source_sink", params);
   }
 
   {
     const std::string kernel_type = "INSFVTKEDSourceSink";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<NonlinearVariableName>("variable") = _tked_name;
     params.set<MooseFunctorName>(NS::TKE) = _tke_name;
     params.set<MooseFunctorName>(NS::density) = _flow_equations_physics->densityName();
     params.set<MooseFunctorName>(NS::mu) = _flow_equations_physics->dynamicViscosityName();
     params.set<MooseFunctorName>(NS::mu_t) = _turbulent_viscosity_name;
     params.set<Real>("C_pl") = getParam<Real>("C_pl");
     params.set<bool>("linearized_model") = getParam<bool>("linearize_sink_sources");
     params.set<std::vector<BoundaryName>>("walls") = _turbulence_walls;
     params.set<MooseEnum>("wall_treatment") = _wall_treatment_eps;
     params.set<Real>("C1_eps") = getParam<Real>("C1_eps");
     params.set<Real>("C2_eps") = getParam<Real>("C2_eps");
     // Currently only Newton method for WCNSFVTurbulencePhysics
     params.set<bool>("newton_solve") = true;
     for (const auto d : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
     getProblem().addFVKernel(kernel_type, prefix() + "tked_source_sink", params);
   }
 }

◆ addKEpsilonTimeDerivatives()

void WCNSFVTurbulencePhysics::addKEpsilonTimeDerivatives ( )

private

Functions adding kernels for the k-epsilon to the k-epsilon equations.

Definition at line 607 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   const std::string kernel_type = "FVFunctorTimeKernel";
   InputParameters params = getFactory().getValidParams(kernel_type);
   assignBlocks(params, _blocks);
 
   params.set<NonlinearVariableName>("variable") = _tke_name;
   if (shouldCreateTimeDerivative(_tke_name, _blocks, /*error if already defined*/ false))
     getProblem().addFVKernel(kernel_type, prefix() + "tke_time", params);
   params.set<NonlinearVariableName>("variable") = _tked_name;
   if (shouldCreateTimeDerivative(_tked_name, _blocks, /*error if already defined*/ false))
     getProblem().addFVKernel(kernel_type, prefix() + "tked_time", params);
 }

◆ addMaterials()

void WCNSFVTurbulencePhysics::addMaterials ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 885 of file WCNSFVTurbulencePhysics.C.

 {
   if (_turbulence_model == "mixing-length")
   {
     const std::string u_names[3] = {"u", "v", "w"};
     InputParameters params =
         getFactory().getValidParams("MixingLengthTurbulentViscosityFunctorMaterial");
     assignBlocks(params, _blocks);
 
     for (const auto d : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[d]) = _velocity_names[d];
 
     params.set<MooseFunctorName>(NS::mixing_length) = _mixing_length_name;
     params.set<MooseFunctorName>(NS::density) = _density_name;
     params.set<MooseFunctorName>(NS::mu) = _dynamic_viscosity_name;
 
     getProblem().addMaterial("MixingLengthTurbulentViscosityFunctorMaterial",
                              prefix() + "mixing_length_material",
                              params);
   }
   else if (_turbulence_model == "k-epsilon")
   {
     if (!getProblem().hasFunctor(NS::mu_eff, /*thread_id=*/0))
     {
       InputParameters params = getFactory().getValidParams("ADParsedFunctorMaterial");
       assignBlocks(params, _blocks);
       const auto mu_name = _flow_equations_physics->dynamicViscosityName();
 
       // Avoid defining floats as functors in the parsed expression
       if (!MooseUtils::isFloat(_turbulent_viscosity_name) && !MooseUtils::isFloat(mu_name))
         params.set<std::vector<std::string>>("functor_names") = {_turbulent_viscosity_name,
                                                                  mu_name};
       else if (MooseUtils::isFloat(_turbulent_viscosity_name) && !MooseUtils::isFloat(mu_name))
         params.set<std::vector<std::string>>("functor_names") = {mu_name};
       else if (!MooseUtils::isFloat(_turbulent_viscosity_name) && MooseUtils::isFloat(mu_name))
         params.set<std::vector<std::string>>("functor_names") = {_turbulent_viscosity_name};
 
       params.set<std::string>("expression") =
           _flow_equations_physics->dynamicViscosityName() + " + " + _turbulent_viscosity_name;
       params.set<std::string>("property_name") = NS::mu_eff;
       getProblem().addMaterial("ADParsedFunctorMaterial", prefix() + "effective_viscosity", params);
     }
     if (!getParam<bool>("mu_t_as_aux_variable"))
     {
       InputParameters params = getFactory().getValidParams("INSFVkEpsilonViscosityFunctorMaterial");
       params.set<MooseFunctorName>(NS::TKE) = _tke_name;
       params.set<MooseFunctorName>(NS::TKED) = _tked_name;
       params.set<MooseFunctorName>(NS::density) = _density_name;
       params.set<ExecFlagEnum>("execute_on") = {EXEC_NONLINEAR};
       if (getParam<bool>("output_mu_t"))
         params.set<std::vector<OutputName>>("outputs") = {"all"};
       getProblem().addMaterial(
           "INSFVkEpsilonViscosityFunctorMaterial", prefix() + "compute_mu_t", params);
     }
 
     if (_has_energy_equation && !getProblem().hasFunctor(NS::k_t, /*thread_id=*/0))
     {
       mooseAssert(!getParam<bool>("k_t_as_aux_variable"), "k_t should exist");
       InputParameters params = getFactory().getValidParams("ADParsedFunctorMaterial");
       assignBlocks(params, _blocks);
       const auto mu_t_name = NS::mu_t;
       const auto cp_name = _fluid_energy_physics->getSpecificHeatName();
       const auto Pr_t_name = getParam<MooseFunctorName>("Pr_t");
 
       // Avoid defining floats as functors in the parsed expression
       if (!MooseUtils::isFloat(cp_name) && !MooseUtils::isFloat(Pr_t_name))
         params.set<std::vector<std::string>>("functor_names") = {cp_name, Pr_t_name, mu_t_name};
       else if (MooseUtils::isFloat(cp_name) && !MooseUtils::isFloat(Pr_t_name))
         params.set<std::vector<std::string>>("functor_names") = {Pr_t_name, mu_t_name};
       else if (!MooseUtils::isFloat(cp_name) && MooseUtils::isFloat(Pr_t_name))
         params.set<std::vector<std::string>>("functor_names") = {cp_name, mu_t_name};
       else
         params.set<std::vector<std::string>>("functor_names") = {mu_t_name};
 
       params.set<std::string>("expression") = mu_t_name + "*" + cp_name + "/" + Pr_t_name;
       params.set<std::string>("property_name") = NS::k_t;
       params.set<ExecFlagEnum>("execute_on") = {EXEC_NONLINEAR};
       params.set<std::vector<OutputName>>("outputs") = {"all"};
       getProblem().addMaterial(
           "ADParsedFunctorMaterial", prefix() + "turbulent_heat_eff_conductivity", params);
     }
 
     if (_has_scalar_equations && !getProblem().hasFunctor(NS::mu_t_passive_scalar, /*thread_id=*/0))
     {
       InputParameters params = getFactory().getValidParams("ADParsedFunctorMaterial");
       assignBlocks(params, _blocks);
       const auto & rho_name = _flow_equations_physics->densityName();
 
       // Avoid defining floats as functors in the parsed expression
       if (!MooseUtils::isFloat(_turbulent_viscosity_name) && !MooseUtils::isFloat(rho_name))
         params.set<std::vector<std::string>>("functor_names") = {_turbulent_viscosity_name,
                                                                  rho_name};
       else if (MooseUtils::isFloat(_turbulent_viscosity_name) && !MooseUtils::isFloat(rho_name))
         params.set<std::vector<std::string>>("functor_names") = {rho_name};
       else if (!MooseUtils::isFloat(_turbulent_viscosity_name) && MooseUtils::isFloat(rho_name))
         params.set<std::vector<std::string>>("functor_names") = {_turbulent_viscosity_name};
 
       const auto turbulent_schmidt_number = getParam<std::vector<Real>>("Sc_t");
       if (turbulent_schmidt_number.size() != 1)
         paramError(
             "passive_scalar_schmidt_number",
             "Only one passive scalar turbulent Schmidt number can be specified with k-epsilon");
       params.set<std::string>("expression") = _turbulent_viscosity_name + "/" + rho_name + " / " +
                                               std::to_string(turbulent_schmidt_number[0]);
       params.set<std::string>("property_name") = NS::mu_t_passive_scalar;
       getProblem().addMaterial(
           "ADParsedFunctorMaterial", prefix() + "scalar_turbulent_diffusivity", params);
     }
   }
 }

◆ addScalarAdvectionTurbulenceKernels()

void WCNSFVTurbulencePhysics::addScalarAdvectionTurbulenceKernels ( )

private

Definition at line 555 of file WCNSFVTurbulencePhysics.C.

Referenced by addFVKernels().

 {
   const auto & passive_scalar_names = _scalar_transport_physics->getAdvectedScalarNames();
   const auto & passive_scalar_schmidt_number = getParam<std::vector<Real>>("Sc_t");
   if (passive_scalar_schmidt_number.size() != passive_scalar_names.size() &&
       passive_scalar_schmidt_number.size() != 1)
     paramError(
         "Sc_t",
         "The number of turbulent Schmidt numbers defined is not equal to the number of passive "
         "scalar fields!");
 
   if (_turbulence_model == "mixing-length")
   {
     const std::string u_names[3] = {"u", "v", "w"};
     const std::string kernel_type = "INSFVMixingLengthScalarDiffusion";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<MooseFunctorName>(NS::mixing_length) = _mixing_length_name;
     for (const auto dim_i : make_range(dimension()))
       params.set<MooseFunctorName>(u_names[dim_i]) = _velocity_names[dim_i];
 
     for (const auto & name_i : index_range(passive_scalar_names))
     {
       params.set<NonlinearVariableName>("variable") = passive_scalar_names[name_i];
       if (passive_scalar_schmidt_number.size() > 1)
         params.set<Real>("schmidt_number") = passive_scalar_schmidt_number[name_i];
       else if (passive_scalar_schmidt_number.size() == 1)
         params.set<Real>("schmidt_number") = passive_scalar_schmidt_number[0];
       else
         params.set<Real>("schmidt_number") = 1.0;
 
       getProblem().addFVKernel(
           kernel_type, prefix() + passive_scalar_names[name_i] + "_mixing_length", params);
     }
   }
   else if (_turbulence_model == "k-epsilon")
   {
     const std::string kernel_type = "FVDiffusion";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
 
     for (const auto & name_i : index_range(passive_scalar_names))
     {
       params.set<NonlinearVariableName>("variable") = passive_scalar_names[name_i];
       params.set<MooseFunctorName>("coeff") = NS::mu_t_passive_scalar;
       getProblem().addFVKernel(
           kernel_type, prefix() + passive_scalar_names[name_i] + "_turbulent_diffusion", params);
     }
   }
 }

◆ addSolverVariables()

void WCNSFVTurbulencePhysics::addSolverVariables ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 333 of file WCNSFVTurbulencePhysics.C.

 {
   if (_turbulence_model == "mixing-length" || _turbulence_model == "none")
     return;
   else if (_turbulence_model == "k-epsilon")
   {
     // Dont add if the user already defined the variable
     // Add turbulent kinetic energy variable
     if (!shouldCreateVariable(_tke_name, _blocks, /*error if aux*/ true))
       reportPotentiallyMissedParameters(
           {"system_names", "tke_scaling", "tke_face_interpolation", "tke_two_term_bc_expansion"},
           "INSFVEnergyVariable");
     else if (_define_variables)
     {
       auto params = getFactory().getValidParams("INSFVEnergyVariable");
       assignBlocks(params, _blocks);
       params.set<std::vector<Real>>("scaling") = {getParam<Real>("tke_scaling")};
       params.set<MooseEnum>("face_interp_method") = getParam<MooseEnum>("tke_face_interpolation");
       params.set<bool>("two_term_boundary_expansion") = getParam<bool>("tke_two_term_bc_expansion");
       params.set<SolverSystemName>("solver_sys") = getSolverSystem(_tke_name);
 
       getProblem().addVariable("INSFVEnergyVariable", _tke_name, params);
     }
     else
       paramError("turbulence_kinetic_energy_variable",
                  "Variable (" + _tke_name +
                      ") supplied to the WCNSFVTurbulencePhysics does not exist!");
 
     // Add turbulent kinetic energy dissipation variable
     if (!shouldCreateVariable(_tked_name, _blocks, /*error if aux*/ true))
       reportPotentiallyMissedParameters(
           {"system_names", "tked_scaling", "tked_face_interpolation", "tked_two_term_bc_expansion"},
           "INSFVEnergyVariable");
     else if (_define_variables)
     {
       auto params = getFactory().getValidParams("INSFVEnergyVariable");
       assignBlocks(params, _blocks);
       params.set<std::vector<Real>>("scaling") = {getParam<Real>("tked_scaling")};
       params.set<MooseEnum>("face_interp_method") = getParam<MooseEnum>("tked_face_interpolation");
       params.set<bool>("two_term_boundary_expansion") =
           getParam<bool>("tked_two_term_bc_expansion");
       params.set<SolverSystemName>("solver_sys") = getSolverSystem(_tked_name);
       getProblem().addVariable("INSFVEnergyVariable", _tked_name, params);
     }
     else
       paramError("turbulence_kinetic_energy_dissipation_variable",
                  "Variable (" + _tked_name +
                      ") supplied to the WCNSFVTurbulencePhysics does not exist!");
   }
 }

◆ densityName()

const MooseFunctorName& WCNSFVCoupledAdvectionPhysicsHelper::densityName ( ) const

inlineinherited

Definition at line 36 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

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

36 { return _density_name; }

WCNSFVCoupledAdvectionPhysicsHelper::_density_name

const MooseFunctorName _density_name

Name of the density material property.

Definition: WCNSFVCoupledAdvectionPhysicsHelper.h:59

◆ dynamicViscosityName()

const MooseFunctorName& WCNSFVCoupledAdvectionPhysicsHelper::dynamicViscosityName ( ) const

inlineinherited

Definition at line 37 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

37 { return _dynamic_viscosity_name; }

WCNSFVCoupledAdvectionPhysicsHelper::_dynamic_viscosity_name

const MooseFunctorName _dynamic_viscosity_name

Name of the dynamic viscosity material property.

Definition: WCNSFVCoupledAdvectionPhysicsHelper.h:61

◆ getAdditionalRMParams()

InputParameters NavierStokesPhysicsBase::getAdditionalRMParams ( ) const

overrideprotectedvirtualinherited

Parameters to change or add relationship managers.

Reimplemented from PhysicsBase.

Definition at line 42 of file NavierStokesPhysicsBase.C.

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

◆ getCoupledFlowPhysics()

const WCNSFVFlowPhysicsBase * WCNSFVCoupledAdvectionPhysicsHelper::getCoupledFlowPhysics ( ) const

inherited

Definition at line 53 of file WCNSFVCoupledAdvectionPhysicsHelper.C.

 {
   // User passed it, just use that
   if (_advection_physics->isParamValid("coupled_flow_physics"))
     return _advection_physics->getCoupledPhysics<WCNSFVFlowPhysicsBase>(
         _advection_physics->getParam<PhysicsName>("coupled_flow_physics"));
   // Look for any physics of the right type, and check the block restriction
   else
   {
     const auto all_flow_physics =
         _advection_physics->getCoupledPhysics<const WCNSFVFlowPhysicsBase>();
     for (const auto physics : all_flow_physics)
       if (_advection_physics->checkBlockRestrictionIdentical(
               physics->name(), physics->blocks(), /*error_if_not_identical=*/false))
       {
         return physics;
       }
   }
   mooseError("No coupled flow Physics found of type derived from 'WCNSFVFlowPhysicsBase'. Use the "
              "'coupled_flow_physics' parameter to give the name of the desired "
              "WCNSFVFlowPhysicsBase-derived Physics to couple with");
 }

◆ getCoupledTurbulencePhysics()

const WCNSFVTurbulencePhysics * WCNSFVCoupledAdvectionPhysicsHelper::getCoupledTurbulencePhysics ( ) const

inherited

Definition at line 77 of file WCNSFVCoupledAdvectionPhysicsHelper.C.

Referenced by WCNSFVFluidHeatTransferPhysicsBase::actOnAdditionalTasks().

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

◆ getNumberAlgebraicGhostingLayersNeeded()

unsigned short WCNSFVTurbulencePhysics::getNumberAlgebraicGhostingLayersNeeded ( ) const

overrideprotectedvirtual

Return the number of ghosting layers needed.

Implements NavierStokesPhysicsBase.

Definition at line 997 of file WCNSFVTurbulencePhysics.C.

 {
   unsigned short ghost_layers = _flow_equations_physics->getNumberAlgebraicGhostingLayersNeeded();
   // due to the computation of the eddy-diffusivity from the strain tensor
   if (_turbulence_model == "mixing-length")
     ghost_layers = std::max(ghost_layers, (unsigned short)3);
   return ghost_layers;
 }

◆ getPorosityFunctorName()

MooseFunctorName WCNSFVCoupledAdvectionPhysicsHelper::getPorosityFunctorName ( bool smoothed ) const

inherited

Return the porosity functor name.

It is important to forward to the Physics so we do not get the smoothing status wrong

Definition at line 47 of file WCNSFVCoupledAdvectionPhysicsHelper.C.

 {
   return _flow_equations_physics->getPorosityFunctorName(smoothed);
 }

◆ hasTurbulenceModel()

bool WCNSFVTurbulencePhysics::hasTurbulenceModel ( ) const

inline

Whether a turbulence model is in use.

Definition at line 32 of file WCNSFVTurbulencePhysics.h.

Referenced by WCNSFVFlowPhysicsBase::hasTurbulencePhysics().

32 { return _turbulence_model != "none"; }

WCNSFVTurbulencePhysics::_turbulence_model

const MooseEnum _turbulence_model

Turbulence model to create the equation(s) for.

Definition: WCNSFVTurbulencePhysics.h:76

◆ initializePhysicsAdditional()

void WCNSFVTurbulencePhysics::initializePhysicsAdditional ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 246 of file WCNSFVTurbulencePhysics.C.

 {
   if (_turbulence_model == "k-epsilon")
     getProblem().needSolutionState(2, Moose::SolutionIterationType::Nonlinear);
 }

◆ retrieveCoupledPhysics()

void WCNSFVTurbulencePhysics::retrieveCoupledPhysics ( )

private

Retrieve the other WCNSFVPhysics at play in the simulation to be able to add the relevant terms (turbulent diffusion notably)

Definition at line 262 of file WCNSFVTurbulencePhysics.C.

Referenced by actOnAdditionalTasks().

 {
   // _flow_equations_physics is initialized by 'WCNSFVCoupledAdvectionPhysicsHelper'
   if (_flow_equations_physics && _flow_equations_physics->hasFlowEquations())
     _has_flow_equations = true;
   else
     _has_flow_equations = false;
 
   // Sanity check for interaction for fluid heat transfer physics
   if (isParamValid("fluid_heat_transfer_physics") && _turbulence_model != "none")
   {
     _fluid_energy_physics = getCoupledPhysics<WCNSFVFluidHeatTransferPhysics>(
         getParam<PhysicsName>("fluid_heat_transfer_physics"), true);
     // Check for a missing parameter / do not support isolated physics for now
     if (!_fluid_energy_physics &&
         !getCoupledPhysics<const WCNSFVFluidHeatTransferPhysics>(true).empty())
       paramError("fluid_heat_transfer_physics",
                  "We currently do not support creating both turbulence physics and fluid heat "
                  "transfer physics that are not coupled together. Use "
                  "'fluid_heat_transfer_physics' to explicitly specify the coupling");
     if (_fluid_energy_physics && _fluid_energy_physics->hasEnergyEquation())
       _has_energy_equation = true;
     else
       _has_energy_equation = false;
   }
   else
   {
     _has_energy_equation = false;
     _fluid_energy_physics = nullptr;
   }
 
   // Sanity check for interaction with scalar transport physics
   if (isParamValid("scalar_transport_physics") && _turbulence_model != "none")
   {
     _scalar_transport_physics = getCoupledPhysics<WCNSFVScalarTransportPhysics>(
         getParam<PhysicsName>("scalar_transport_physics"), true);
     if (!_scalar_transport_physics &&
         !getCoupledPhysics<const WCNSFVScalarTransportPhysics>(true).empty())
       paramError(
           "scalar_transport_physics",
           "We currently do not support creating both turbulence physics and scalar transport "
           "physics that are not coupled together");
     if (_scalar_transport_physics && _scalar_transport_physics->hasScalarEquations())
       _has_scalar_equations = true;
     else
       _has_scalar_equations = false;
   }
   else
   {
     _has_scalar_equations = false;
     _scalar_transport_physics = nullptr;
   }
 
   // To help remediate the danger of the parameter setup
   if (_verbose)
   {
     if (_has_energy_equation)
       mooseInfoRepeated("Coupling turbulence physics with fluid heat transfer physics " +
                         _fluid_energy_physics->name());
     else
       mooseInfoRepeated("No fluid heat transfer equation considered by this turbulence "
                         "physics.");
     if (_has_scalar_equations)
       mooseInfoRepeated("Coupling turbulence physics with scalar transport physics " +
                         _scalar_transport_physics->name());
     else
       mooseInfoRepeated("No scalar transport equations considered by this turbulence physics.");
   }
 }

◆ tkeName()

MooseFunctorName WCNSFVTurbulencePhysics::tkeName ( ) const

inline

The name of the turbulent kinetic energy variable.

Definition at line 40 of file WCNSFVTurbulencePhysics.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyWallBC().

40 { return _tke_name; }

WCNSFVTurbulencePhysics::_tke_name

const VariableName _tke_name

Name of the turbulent kinetic energy.

Definition: WCNSFVTurbulencePhysics.h:97

◆ turbulenceEpsilonWallTreatment()

MooseEnum WCNSFVTurbulencePhysics::turbulenceEpsilonWallTreatment ( ) const

inline

The turbulence epsilon wall treatment (same for all turbulence walls currently)

Definition at line 36 of file WCNSFVTurbulencePhysics.h.

36 { return _wall_treatment_eps; }

WCNSFVTurbulencePhysics::_wall_treatment_eps

MooseEnum _wall_treatment_eps

Turbulence wall treatment for epsilon (same for all walls currently)

Definition: WCNSFVTurbulencePhysics.h:93

◆ turbulenceTemperatureWallTreatment()

MooseEnum WCNSFVTurbulencePhysics::turbulenceTemperatureWallTreatment ( ) const

inline

The turbulence temperature wall treatment (same for all turbulence walls currently)

Definition at line 38 of file WCNSFVTurbulencePhysics.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyWallBC().

38 { return _wall_treatment_temp; }

WCNSFVTurbulencePhysics::_wall_treatment_temp

MooseEnum _wall_treatment_temp

Turbulence wall treatment for temperature (same for all walls currently)

Definition: WCNSFVTurbulencePhysics.h:95

◆ turbulenceWalls()

std::vector<BoundaryName> WCNSFVTurbulencePhysics::turbulenceWalls ( ) const

inline

The names of the boundaries with turbulence wall functions.

Definition at line 34 of file WCNSFVTurbulencePhysics.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyWallBC().

34 { return _turbulence_walls; }

WCNSFVTurbulencePhysics::_turbulence_walls

std::vector< BoundaryName > _turbulence_walls

List of boundaries to act as walls for turbulence models.

Definition: WCNSFVTurbulencePhysics.h:91

◆ usingNavierStokesFVSyntax()

bool NavierStokesPhysicsBase::usingNavierStokesFVSyntax ( ) const

inlineprotectedinherited

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

Definition at line 32 of file NavierStokesPhysicsBase.h.

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

◆ validParams()

InputParameters WCNSFVTurbulencePhysics::validParams ( )

static

Definition at line 30 of file WCNSFVTurbulencePhysics.C.

 {
   InputParameters params = NavierStokesPhysicsBase::validParams();
   params += WCNSFVCoupledAdvectionPhysicsHelper::validParams();
   params.addClassDescription(
       "Define a turbulence model for a incompressible or weakly-compressible Navier Stokes "
       "flow with a finite volume discretization");
 
   MooseEnum turbulence_type("mixing-length k-epsilon none", "none");
   params.addParam<MooseEnum>(
       "turbulence_handling",
       turbulence_type,
       "The way turbulent diffusivities are determined in the turbulent regime.");
   params += NSFVBase::commonTurbulenceParams();
   params.transferParam<bool>(NSFVBase::validParams(), "mixing_length_two_term_bc_expansion");
 
   // TODO Added to facilitate transition, remove default once NavierStokesFV action is removed
   params.addParam<AuxVariableName>(
       "mixing_length_name", "mixing_length", "Name of the mixing length auxiliary variable");
   params.deprecateParam("mixing_length_walls", "turbulence_walls", "");
 
   // Not implemented, re-enable with k-epsilon
   params.suppressParameter<MooseEnum>("preconditioning");
 
   // K-Epsilon parameters
   params.addParam<MooseFunctorName>(
       "tke_name", NS::TKE, "Name of the turbulent kinetic energy variable");
   params.addParam<MooseFunctorName>(
       "tked_name", NS::TKED, "Name of the turbulent kinetic energy dissipation variable");
   params.addParam<FunctionName>(
       "initial_tke", "0", "Initial value for the turbulence kinetic energy");
   params.addParam<FunctionName>(
       "initial_tked", "0", "Initial value for the turbulence kinetic energy dissipation");
   params.addParam<FunctionName>("initial_mu_t", "Initial value for the turbulence viscosity");
 
   params.addParam<Real>("C1_eps",
                         "C1 coefficient for the turbulent kinetic energy dissipation equation");
   params.addParam<Real>("C2_eps",
                         "C2 coefficient for the turbulent kinetic energy dissipation equation");
   params.addParam<MooseFunctorName>(
       "sigma_k", "Scaling coefficient for the turbulent kinetic energy diffusion term");
   params.addParam<MooseFunctorName>(
       "sigma_eps",
       "Scaling coefficient for the turbulent kinetic energy dissipation diffusion term");
   params.addParam<MooseFunctorName>(
       NS::turbulent_Prandtl, NS::turbulent_Prandtl, "Turbulent Prandtl number");
   params.transferParam<Real>(INSFVTKESourceSink::validParams(), "C_pl");
 
   // Boundary parameters
   params.addParam<bool>("bulk_wall_treatment", true, "Whether to treat the wall cell as bulk");
   MooseEnum wall_treatment("eq_newton eq_incremental eq_linearized neq", "neq");
   params.addParam<MooseEnum>("wall_treatment_eps",
                              wall_treatment,
                              "The method used for computing the epsilon wall functions");
   params.addParam<MooseEnum>("wall_treatment_T",
                              wall_treatment,
                              "The method used for computing the temperature wall functions");
   params.transferParam<Real>(INSFVTurbulentViscosityWallFunction::validParams(), "C_mu");
 
   // K-Epsilon numerical scheme parameters
   MooseEnum face_interpol_types("average skewness-corrected", "average");
   MooseEnum adv_interpol_types("average upwind", "upwind");
   params.addRangeCheckedParam<Real>(
       "tke_scaling",
       1.0,
       "tke_scaling > 0.0",
       "The scaling factor for the turbulent kinetic energy equation.");
   params.addParam<MooseEnum>("tke_face_interpolation",
                              face_interpol_types,
                              "The numerical scheme to interpolate the TKE to the "
                              "face (separate from the advected quantity interpolation).");
   params.addParam<MooseEnum>("tke_advection_interpolation",
                              adv_interpol_types,
                              "The numerical scheme to interpolate the TKE to the "
                              "face when in the advection kernel.");
   // Better Jacobian if not linearizing sink and sources
   params.addParam<bool>("linearize_sink_sources", false, "Whether to linearize the source term");
   // Better convergence on some cases when neglecting advection derivatives
   params.addParam<bool>(
       "neglect_advection_derivatives",
       false,
       "Whether to remove the off-diagonal velocity term in the TKE and TKED advection term");
   params.addParam<bool>(
       "tke_two_term_bc_expansion",
       true,
       "If a two-term Taylor expansion is needed for the determination of the boundary values"
       "of the turbulent kinetic energy.");
   params.addRangeCheckedParam<Real>(
       "tked_scaling",
       1.0,
       "tked_scaling > 0.0",
       "The scaling factor for the turbulent kinetic energy dissipation equation.");
   params.addParam<MooseEnum>("tked_face_interpolation",
                              face_interpol_types,
                              "The numerical scheme to interpolate the TKED to the "
                              "face (separate from the advected quantity interpolation).");
   params.addParam<MooseEnum>("tked_advection_interpolation",
                              adv_interpol_types,
                              "The numerical scheme to interpolate the TKED to the "
                              "face when in the advection kernel.");
   params.addParam<bool>(
       "tked_two_term_bc_expansion",
       true,
       "If a two-term Taylor expansion is needed for the determination of the boundary values"
       "of the turbulent kinetic energy dissipation.");
   params.addParam<bool>(
       "turbulent_viscosity_two_term_bc_expansion",
       true,
       "If a two-term Taylor expansion is needed for the determination of the boundary values"
       "of the turbulent viscosity.");
   MooseEnum coeff_interp_method("average harmonic", "harmonic");
   params.addParam<MooseEnum>("turbulent_viscosity_interp_method",
                              coeff_interp_method,
                              "Face interpolation method for the turbulent viscosity");
   params.addParam<bool>("mu_t_as_aux_variable",
                         false,
                         "Whether to use an auxiliary variable instead of a functor material "
                         "property for the turbulent viscosity");
   params.addParam<bool>("output_mu_t", true, "Whether to add mu_t to the field outputs");
   params.addParam<bool>("k_t_as_aux_variable",
                         false,
                         "Whether to use an auxiliary variable for the turbulent conductivity");
 
   // Add the coupled physics
   // TODO Remove the defaults once NavierStokesFV action is removed
   // It is a little risky right now because the user could forget to pass the parameter and
   // be missing the influence of turbulence on either of these physics. There is a check in the
   // constructor to present this from happening
   params.addParam<PhysicsName>(
       "fluid_heat_transfer_physics",
       "NavierStokesFV",
       "WCNSFVFluidHeatTransferPhysics generating the heat advection equations");
   params.addParam<PhysicsName>(
       "scalar_transport_physics",
       "NavierStokesFV",
       "WCNSFVScalarTransportPhysics generating the scalar advection equations");
 
   // Parameter groups
   params.addParamNamesToGroup("mixing_length_name mixing_length_two_term_bc_expansion",
                               "Mixing length model");
   params.addParamNamesToGroup("fluid_heat_transfer_physics turbulent_prandtl "
                               "scalar_transport_physics Sc_t",
                               "Coupled Physics");
   params.addParamNamesToGroup("initial_tke initial_tked C1_eps C2_eps sigma_k sigma_eps",
                               "K-Epsilon model");
   params.addParamNamesToGroup("C_mu bulk_wall_treatment wall_treatment_eps wall_treatment_T",
                               "K-Epsilon wall function");
   params.addParamNamesToGroup(
       "tke_scaling tke_face_interpolation tke_two_term_bc_expansion tked_scaling "
       "tked_face_interpolation tked_two_term_bc_expansion "
       "turbulent_viscosity_two_term_bc_expansion turbulent_viscosity_interp_method "
       "mu_t_as_aux_variable k_t_as_aux_variable linearize_sink_sources",
       "K-Epsilon model numerical");
 
   return params;
 }

Member Data Documentation

◆ _advection_physics

const NavierStokesPhysicsBase* WCNSFVCoupledAdvectionPhysicsHelper::_advection_physics

protectedinherited

The Physics class using this helper.

Definition at line 41 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVCoupledAdvectionPhysicsHelper::getCoupledFlowPhysics(), and WCNSFVCoupledAdvectionPhysicsHelper::getCoupledTurbulencePhysics().

◆ _compressibility

const MooseEnum WCNSFVCoupledAdvectionPhysicsHelper::_compressibility

protectedinherited

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

Definition at line 48 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyTimeKernels().

◆ _define_variables

bool NavierStokesPhysicsBase::_define_variables

protectedinherited

Whether to define variables if they do not exist.

Definition at line 44 of file NavierStokesPhysicsBase.h.

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

◆ _density_name

const MooseFunctorName WCNSFVCoupledAdvectionPhysicsHelper::_density_name

protectedinherited

Name of the density material property.

Definition at line 59 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyInletBC(), WCNSFVFluidHeatTransferPhysics::addEnergyTimeKernels(), WCNSLinearFVFluidHeatTransferPhysics::addEnergyTimeKernels(), WCNSFVFluidHeatTransferPhysics::addEnergyWallBC(), addFlowTurbulenceKernels(), addFluidEnergyTurbulenceKernels(), addMaterials(), WCNSFVFluidHeatTransferPhysicsBase::addMaterials(), WCNSFVScalarTransportPhysics::addScalarInletBC(), and WCNSFVCoupledAdvectionPhysicsHelper::densityName().

◆ _dynamic_viscosity_name

const MooseFunctorName WCNSFVCoupledAdvectionPhysicsHelper::_dynamic_viscosity_name

protectedinherited

Name of the dynamic viscosity material property.

Definition at line 61 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyWallBC(), addMaterials(), and WCNSFVCoupledAdvectionPhysicsHelper::dynamicViscosityName().

◆ _flow_equations_physics

const WCNSFVFlowPhysicsBase* WCNSFVCoupledAdvectionPhysicsHelper::_flow_equations_physics

protectedinherited

Flow physics.

Definition at line 43 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

◆ _fluid_energy_physics

const WCNSFVFluidHeatTransferPhysics* WCNSFVTurbulencePhysics::_fluid_energy_physics

private

The heat advection physics to add turbulent mixing for.

Definition at line 83 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addFluidEnergyTurbulenceKernels(), addMaterials(), and retrieveCoupledPhysics().

◆ _has_energy_equation

bool WCNSFVTurbulencePhysics::_has_energy_equation

private

Definition at line 79 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addFVKernels(), addMaterials(), and retrieveCoupledPhysics().

◆ _has_flow_equations

bool WCNSFVTurbulencePhysics::_has_flow_equations

private

Definition at line 78 of file WCNSFVTurbulencePhysics.h.

Referenced by addFVKernels(), and retrieveCoupledPhysics().

◆ _has_scalar_equations

bool WCNSFVTurbulencePhysics::_has_scalar_equations

private

Definition at line 80 of file WCNSFVTurbulencePhysics.h.

Referenced by addFVKernels(), addMaterials(), and retrieveCoupledPhysics().

◆ _mixing_length_name

const VariableName WCNSFVTurbulencePhysics::_mixing_length_name

private

Name of the mixing length auxiliary variable.

Definition at line 89 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addAuxiliaryVariables(), addFlowTurbulenceKernels(), addFluidEnergyTurbulenceKernels(), addMaterials(), addScalarAdvectionTurbulenceKernels(), and WCNSFVTurbulencePhysics().

◆ _porous_medium_treatment

const bool WCNSFVCoupledAdvectionPhysicsHelper::_porous_medium_treatment

protectedinherited

Switch to show if porous medium treatment is requested or not.

Definition at line 51 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyAdvectionKernels(), WCNSFVFluidHeatTransferPhysics::addEnergyHeatConductionKernels(), WCNSFVFluidHeatTransferPhysics::addEnergyTimeKernels(), addFlowTurbulenceKernels(), addFluidEnergyTurbulenceKernels(), WCNSFVScalarTransportPhysics::addScalarAdvectionKernels(), WCNSFVFluidHeatTransferPhysicsBase::processThermalConductivity(), WCNSLinearFVFluidHeatTransferPhysics::WCNSLinearFVFluidHeatTransferPhysics(), and WCNSLinearFVScalarTransportPhysics::WCNSLinearFVScalarTransportPhysics().

◆ _pressure_name

const NonlinearVariableName WCNSFVCoupledAdvectionPhysicsHelper::_pressure_name

protectedinherited

Pressure name.

Definition at line 56 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

◆ _scalar_transport_physics

const WCNSFVScalarTransportPhysics* WCNSFVTurbulencePhysics::_scalar_transport_physics

private

The scalar advection physics to add turbulent mixing for.

Definition at line 85 of file WCNSFVTurbulencePhysics.h.

Referenced by addScalarAdvectionTurbulenceKernels(), and retrieveCoupledPhysics().

◆ _tke_name

const VariableName WCNSFVTurbulencePhysics::_tke_name

private

Name of the turbulent kinetic energy.

Definition at line 97 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addAuxiliaryVariables(), addFVBCs(), addInitialConditions(), addKEpsilonAdvection(), addKEpsilonDiffusion(), addKEpsilonSink(), addKEpsilonTimeDerivatives(), addMaterials(), addSolverVariables(), tkeName(), and WCNSFVTurbulencePhysics().

◆ _tked_name

const VariableName WCNSFVTurbulencePhysics::_tked_name

private

Name of the turbulent kinetic energy dissipation.

Definition at line 99 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addInitialConditions(), addKEpsilonAdvection(), addKEpsilonDiffusion(), addKEpsilonSink(), addKEpsilonTimeDerivatives(), addMaterials(), addSolverVariables(), and WCNSFVTurbulencePhysics().

◆ _turbulence_model

const MooseEnum WCNSFVTurbulencePhysics::_turbulence_model

private

Turbulence model to create the equation(s) for.

Definition at line 76 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addAuxiliaryVariables(), addFlowTurbulenceKernels(), addFluidEnergyTurbulenceKernels(), addFVBCs(), addFVKernels(), addInitialConditions(), addMaterials(), addScalarAdvectionTurbulenceKernels(), addSolverVariables(), getNumberAlgebraicGhostingLayersNeeded(), hasTurbulenceModel(), initializePhysicsAdditional(), retrieveCoupledPhysics(), and WCNSFVTurbulencePhysics().

◆ _turbulence_physics

const WCNSFVTurbulencePhysics* WCNSFVCoupledAdvectionPhysicsHelper::_turbulence_physics

protectedinherited

Turbulence.

Definition at line 45 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVFluidHeatTransferPhysicsBase::actOnAdditionalTasks(), and WCNSFVFluidHeatTransferPhysics::addEnergyWallBC().

◆ _turbulence_walls

std::vector<BoundaryName> WCNSFVTurbulencePhysics::_turbulence_walls

private

List of boundaries to act as walls for turbulence models.

Definition at line 91 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addFVBCs(), addKEpsilonAdvection(), addKEpsilonDiffusion(), addKEpsilonSink(), and turbulenceWalls().

◆ _turbulent_viscosity_name

const VariableName WCNSFVTurbulencePhysics::_turbulent_viscosity_name = NS::mu_t

private

Name of the turbulence viscosity auxiliary variable (or property)

Definition at line 101 of file WCNSFVTurbulencePhysics.h.

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

◆ _velocity_interpolation

const MooseEnum WCNSFVCoupledAdvectionPhysicsHelper::_velocity_interpolation

protectedinherited

The velocity / momentum face interpolation method for advecting other quantities.

Definition at line 64 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by WCNSFVFluidHeatTransferPhysics::addEnergyAdvectionKernels(), addKEpsilonAdvection(), and WCNSFVScalarTransportPhysics::addScalarAdvectionKernels().

◆ _velocity_names

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

protectedinherited

Velocity names.

Definition at line 54 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

Referenced by addAuxiliaryKernels(), WCNSFVFluidHeatTransferPhysics::addEnergyInletBC(), WCNSFVFluidHeatTransferPhysics::addEnergyWallBC(), addFlowTurbulenceKernels(), addFluidEnergyTurbulenceKernels(), addFVBCs(), addKEpsilonSink(), addMaterials(), addScalarAdvectionTurbulenceKernels(), and WCNSFVScalarTransportPhysics::addScalarInletBC().

◆ _wall_treatment_eps

MooseEnum WCNSFVTurbulencePhysics::_wall_treatment_eps

private

Turbulence wall treatment for epsilon (same for all walls currently)

Definition at line 93 of file WCNSFVTurbulencePhysics.h.

Referenced by addAuxiliaryKernels(), addFVBCs(), addKEpsilonSink(), and turbulenceEpsilonWallTreatment().

◆ _wall_treatment_temp

MooseEnum WCNSFVTurbulencePhysics::_wall_treatment_temp

private

Turbulence wall treatment for temperature (same for all walls currently)

Definition at line 95 of file WCNSFVTurbulencePhysics.h.

Referenced by turbulenceTemperatureWallTreatment().

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

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

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ WCNSFVTurbulencePhysics()

Member Function Documentation

◆ actOnAdditionalTasks()

◆ addAuxiliaryKernels()

◆ addAuxiliaryVariables()

◆ addFlowTurbulenceKernels()

◆ addFluidEnergyTurbulenceKernels()

◆ addFVBCs()

◆ addFVKernels()

◆ addInitialConditions()

◆ addKEpsilonAdvection()

◆ addKEpsilonDiffusion()

◆ addKEpsilonSink()

◆ addKEpsilonTimeDerivatives()

◆ addMaterials()

◆ addScalarAdvectionTurbulenceKernels()

◆ addSolverVariables()

◆ densityName()

◆ dynamicViscosityName()

◆ getAdditionalRMParams()

◆ getCoupledFlowPhysics()

◆ getCoupledTurbulencePhysics()

◆ getNumberAlgebraicGhostingLayersNeeded()

◆ getPorosityFunctorName()

◆ hasTurbulenceModel()

◆ initializePhysicsAdditional()

◆ retrieveCoupledPhysics()

◆ tkeName()

◆ turbulenceEpsilonWallTreatment()

◆ turbulenceTemperatureWallTreatment()

◆ turbulenceWalls()

◆ usingNavierStokesFVSyntax()

◆ validParams()

Member Data Documentation

◆ _advection_physics

◆ _compressibility

◆ _define_variables

◆ _density_name

◆ _dynamic_viscosity_name

◆ _flow_equations_physics

◆ _fluid_energy_physics

◆ _has_energy_equation

◆ _has_flow_equations

◆ _has_scalar_equations

◆ _mixing_length_name

◆ _porous_medium_treatment

◆ _pressure_name

◆ _scalar_transport_physics

◆ _tke_name

◆ _tked_name

◆ _turbulence_model

◆ _turbulence_physics

◆ _turbulence_walls

◆ _turbulent_viscosity_name

◆ _velocity_interpolation

◆ _velocity_names

◆ _wall_treatment_eps

◆ _wall_treatment_temp