Creates all the objects needed to solve the Navier Stokes scalar transport equations using the nonlinear finite volume weakly-compressible discretization (WCNSFV) More...

#include <WCNSFVScalarTransportPhysics.h>

Inheritance diagram for WCNSFVScalarTransportPhysics:

[legend]

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	WCNSFVScalarTransportPhysics (const InputParameters &parameters)

const std::vector< NonlinearVariableName > &	getAdvectedScalarNames () const
	Get the names of the advected scalar quantity variables. More...

bool	hasScalarEquations () const
	Whether the physics is actually creating the scalar advection equations. More...

virtual void	act () override final

virtual void	actOnAdditionalTasks ()

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

virtual void	addFVBCs () override

virtual void	setSlipVelocityParams (InputParameters &) const

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
std::vector< NonlinearVariableName >	_passive_scalar_names
	Names of the passive scalar variables. More...

const bool	_has_scalar_equation
	A boolean to help compatibility with the old Modules/NavierStokesFV syntax or to deliberately skip adding the equations (for example for mixtures with a stationary phase) More...

MultiMooseEnum	_passive_scalar_inlet_types
	Passive scalar inlet boundary types. More...

std::vector< std::vector< MooseFunctorName > >	_passive_scalar_inlet_functors
	Functors describing the inlet boundary values. See passive_scalar_inlet_types for what the functors actually represent. More...

std::vector< MooseFunctorName >	_passive_scalar_sources
	Functors for the passive scalar sources. Indexing is scalar variable index. More...

std::vector< std::vector< MooseFunctorName > >	_passive_scalar_coupled_sources
	Functors for the passive scalar (coupled) sources. Inner indexing is scalar variable index. More...

std::vector< std::vector< Real > >	_passive_scalar_sources_coef
	Coefficients multiplying for the passive scalar sources. Inner indexing is scalar variable index. More...

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

std::vector< SolverSystemName >	_system_names

std::vector< unsigned int >	_system_numbers

const bool	_verbose

const MooseEnum &	_preconditioning

std::vector< SubdomainName >	_blocks

std::string	_registered_identifier

std::string	_specific_task_name

std::set< std::string >	_all_tasks

ActionWarehouse &	_awh

const std::string &	_current_task

std::shared_ptr< MooseMesh > &	_mesh

std::shared_ptr< MooseMesh > &	_displaced_mesh

std::shared_ptr< FEProblemBase > &	_problem

PerfID	_act_timer

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

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

virtual void	addScalarTimeKernels () override
	Functions adding kernels for the incompressible / weakly-compressible scalar transport equation If the material properties are not constant, some of these can be used for weakly-compressible simulations as well. More...

virtual void	addScalarDiffusionKernels () override

virtual void	addScalarAdvectionKernels () override

virtual void	addScalarSourceKernels () override
	Equivalent of NSFVAction addScalarCoupledSourceKernels. More...

virtual void	addScalarInletBC () override
	Functions adding boundary conditions for the incompressible simulation. More...

virtual void	addScalarWallBC () override

virtual void	addScalarOutletBC () override

Detailed Description

Creates all the objects needed to solve the Navier Stokes scalar transport equations using the nonlinear finite volume weakly-compressible discretization (WCNSFV)

Definition at line 24 of file WCNSFVScalarTransportPhysics.h.

Constructor & Destructor Documentation

◆ WCNSFVScalarTransportPhysics()

WCNSFVScalarTransportPhysics::WCNSFVScalarTransportPhysics ( const InputParameters & parameters )

Definition at line 29 of file WCNSFVScalarTransportPhysics.C.

   : WCNSFVScalarTransportPhysicsBase(parameters)
 {
 }

Member Function Documentation

◆ addFVBCs()

void WCNSFVScalarTransportPhysicsBase::addFVBCs ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Definition at line 130 of file WCNSFVScalarTransportPhysicsBase.C.

 {
   // For compatibility with Modules/NavierStokesFV syntax
   if (!_has_scalar_equation)
     return;
 
   addScalarInletBC();
   // There is typically no wall flux of passive scalars, similarly we rarely know
   // their concentrations at the outlet at the beginning of the simulation
   // TODO: we will know the outlet values in case of flow reversal. Implement scalar outlet
   addScalarWallBC();
   addScalarOutletBC();
 }

◆ addFVKernels()

void WCNSFVScalarTransportPhysicsBase::addFVKernels ( )

overrideprotectedvirtualinherited

Reimplemented from PhysicsBase.

Reimplemented in WCNSFVTwoPhaseMixturePhysics, and WCNSLinearFVTwoPhaseMixturePhysics.

Definition at line 114 of file WCNSFVScalarTransportPhysicsBase.C.

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

 {
   // For compatibility with Modules/NavierStokesFV syntax
   if (!_has_scalar_equation)
     return;
 
   if (isTransient())
     addScalarTimeKernels();
 
   addScalarAdvectionKernels();
   addScalarDiffusionKernels();
   if (_passive_scalar_sources.size() || _passive_scalar_coupled_sources.size())
     addScalarSourceKernels();
 }

◆ addScalarAdvectionKernels()

void WCNSFVScalarTransportPhysics::addScalarAdvectionKernels ( )

overrideprivatevirtual

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 83 of file WCNSFVScalarTransportPhysics.C.

 {
   const std::string kernel_type =
       _porous_medium_treatment ? "PINSFVScalarFieldAdvection" : "INSFVScalarFieldAdvection";
   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<MooseEnum>("advected_interp_method") =
       getParam<MooseEnum>("passive_scalar_advection_interpolation");
   setSlipVelocityParams(params);
   if (_porous_medium_treatment)
     params.set<MooseFunctorName>(NS::porosity) =
         _flow_equations_physics->getPorosityFunctorName(/*smoothed=*/true);
 
   for (const auto & vname : _passive_scalar_names)
   {
     params.set<NonlinearVariableName>("variable") = vname;
     getProblem().addFVKernel(kernel_type, prefix() + "ins_" + vname + "_advection", params);
   }
 }

◆ addScalarDiffusionKernels()

void WCNSFVScalarTransportPhysics::addScalarDiffusionKernels ( )

overrideprivatevirtual

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 107 of file WCNSFVScalarTransportPhysics.C.

 {
   // Direct specification of diffusion term
   const auto passive_scalar_diffusivities =
       getParam<std::vector<MooseFunctorName>>("passive_scalar_diffusivity");
 
   if (passive_scalar_diffusivities.size())
   {
     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") = passive_scalar_diffusivities[name_i];
       getProblem().addFVKernel(
           kernel_type, prefix() + "ins_" + _passive_scalar_names[name_i] + "_diffusion", params);
     }
   }
 }

◆ addScalarInletBC()

void WCNSFVScalarTransportPhysics::addScalarInletBC ( )

overrideprivatevirtual

Functions adding boundary conditions for the incompressible simulation.

These are used for weakly-compressible simulations as well.

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 162 of file WCNSFVScalarTransportPhysics.C.

 {
   const auto & inlet_boundaries = _flow_equations_physics->getInletBoundaries();
   if (inlet_boundaries.empty())
     return;
 
   // Boundary checks
   // TODO: once we have vectors of MooseEnum, we could use the same templated check for types and
   // functors
   if (inlet_boundaries.size() * _passive_scalar_names.size() != _passive_scalar_inlet_types.size())
     paramError(
         "passive_scalar_inlet_types",
         "The number of scalar inlet types (" + std::to_string(_passive_scalar_inlet_types.size()) +
             ") is not equal to the number of inlet boundaries (" +
             std::to_string(inlet_boundaries.size()) + ") times the number of passive scalars (" +
             std::to_string(_passive_scalar_names.size()) + ")");
   if (_passive_scalar_names.size() != _passive_scalar_inlet_functors.size())
     paramError("passive_scalar_inlet_functors",
                "The number of groups of inlet functors (" +
                    std::to_string(_passive_scalar_inlet_functors.size()) +
                    ") is not equal to the number of passive scalars (" +
                    std::to_string(_passive_scalar_names.size()) + ")");
 
   for (const auto name_i : index_range(_passive_scalar_names))
   {
     if (inlet_boundaries.size() != _passive_scalar_inlet_functors[name_i].size())
       paramError("passive_scalar_inlet_functors",
                  "The number of inlet boundary functors for scalar '" +
                      _passive_scalar_names[name_i] +
                      "' does not match the number of inlet boundaries (" +
                      std::to_string(_passive_scalar_inlet_functors[name_i].size()) + ")");
 
     unsigned int flux_bc_counter = 0;
     unsigned int num_inlets = inlet_boundaries.size();
     for (unsigned int bc_ind = 0; bc_ind < num_inlets; ++bc_ind)
     {
       if (_passive_scalar_inlet_types[name_i * num_inlets + bc_ind] == "fixed-value")
       {
         const std::string bc_type = "FVFunctionDirichletBC";
         InputParameters params = getFactory().getValidParams(bc_type);
         params.set<NonlinearVariableName>("variable") = _passive_scalar_names[name_i];
         params.set<FunctionName>("function") = _passive_scalar_inlet_functors[name_i][bc_ind];
         params.set<std::vector<BoundaryName>>("boundary") = {inlet_boundaries[bc_ind]};
 
         getProblem().addFVBC(
             bc_type, _passive_scalar_names[name_i] + "_" + inlet_boundaries[bc_ind], params);
       }
       else if (_passive_scalar_inlet_types[name_i * num_inlets + bc_ind] == "flux-mass" ||
                _passive_scalar_inlet_types[name_i * num_inlets + bc_ind] == "flux-velocity")
       {
         const auto flux_inlet_directions = _flow_equations_physics->getFluxInletDirections();
         const auto flux_inlet_pps = _flow_equations_physics->getFluxInletPPs();
 
         const std::string bc_type = "WCNSFVScalarFluxBC";
         InputParameters params = getFactory().getValidParams(bc_type);
         params.set<NonlinearVariableName>("variable") = _passive_scalar_names[name_i];
         params.set<MooseFunctorName>("passive_scalar") = _passive_scalar_names[name_i];
         if (flux_inlet_directions.size())
           params.set<Point>("direction") = flux_inlet_directions[flux_bc_counter];
         if (_passive_scalar_inlet_types[name_i * num_inlets + bc_ind] == "flux-mass")
         {
           params.set<PostprocessorName>("mdot_pp") = flux_inlet_pps[flux_bc_counter];
           params.set<PostprocessorName>("area_pp") = "area_pp_" + inlet_boundaries[bc_ind];
         }
         else
           params.set<PostprocessorName>("velocity_pp") = flux_inlet_pps[flux_bc_counter];
 
         params.set<MooseFunctorName>(NS::density) = _density_name;
         params.set<PostprocessorName>("scalar_value_pp") =
             _passive_scalar_inlet_functors[name_i][bc_ind];
         params.set<std::vector<BoundaryName>>("boundary") = {inlet_boundaries[bc_ind]};
 
         params.set<MooseFunctorName>(NS::velocity_x) = _velocity_names[0];
         if (dimension() > 1)
           params.set<MooseFunctorName>(NS::velocity_y) = _velocity_names[1];
         if (dimension() > 2)
           params.set<MooseFunctorName>(NS::velocity_z) = _velocity_names[2];
 
         getProblem().addFVBC(bc_type,
                              prefix() + _passive_scalar_names[name_i] + "_" +
                                  inlet_boundaries[bc_ind],
                              params);
         flux_bc_counter += 1;
       }
     }
   }
 }

◆ addScalarOutletBC()

void WCNSFVScalarTransportPhysics::addScalarOutletBC ( )

overrideprivatevirtual

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 251 of file WCNSFVScalarTransportPhysics.C.

 {
   // Advection outlet is naturally handled by the advection flux kernel
   return;
 }

◆ addScalarSourceKernels()

void WCNSFVScalarTransportPhysics::addScalarSourceKernels ( )

overrideprivatevirtual

Equivalent of NSFVAction addScalarCoupledSourceKernels.

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 129 of file WCNSFVScalarTransportPhysics.C.

 {
   const std::string kernel_type = "FVCoupledForce";
   InputParameters params = getFactory().getValidParams(kernel_type);
   assignBlocks(params, _blocks);
 
   for (const auto scalar_i : index_range(_passive_scalar_names))
   {
     params.set<NonlinearVariableName>("variable") = _passive_scalar_names[scalar_i];
     if (_passive_scalar_sources.size())
     {
       // Added for backward compatibility with former Modules/NavierStokesFV syntax
       params.set<MooseFunctorName>("v") = _passive_scalar_sources[scalar_i];
       getProblem().addFVKernel(
           kernel_type, prefix() + "ins_" + _passive_scalar_names[scalar_i] + "_source", params);
     }
 
     // Sufficient for all intents and purposes
     if (_passive_scalar_coupled_sources.size())
       for (const auto i : index_range(_passive_scalar_coupled_sources[scalar_i]))
       {
         params.set<MooseFunctorName>("v") = _passive_scalar_coupled_sources[scalar_i][i];
         if (_passive_scalar_sources_coef.size())
           params.set<Real>("coef") = _passive_scalar_sources_coef[scalar_i][i];
         getProblem().addFVKernel(kernel_type,
                                  prefix() + "ins_" + _passive_scalar_names[scalar_i] +
                                      "_coupled_source_" + std::to_string(i),
                                  params);
       }
   }
 }

◆ addScalarTimeKernels()

void WCNSFVScalarTransportPhysics::addScalarTimeKernels ( )

overrideprivatevirtual

Functions adding kernels for the incompressible / weakly-compressible scalar transport equation If the material properties are not constant, some of these can be used for weakly-compressible simulations as well.

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 68 of file WCNSFVScalarTransportPhysics.C.

 {
   for (const auto & vname : _passive_scalar_names)
   {
     const std::string kernel_type = "FVFunctorTimeKernel";
     InputParameters params = getFactory().getValidParams(kernel_type);
     assignBlocks(params, _blocks);
     params.set<NonlinearVariableName>("variable") = vname;
 
     if (shouldCreateTimeDerivative(vname, _blocks, /* error if already defined */ false))
       getProblem().addFVKernel(kernel_type, prefix() + "ins_" + vname + "_time", params);
   }
 }

◆ addScalarWallBC()

virtual void WCNSFVScalarTransportPhysics::addScalarWallBC ( )

inlineoverrideprivatevirtual

Implements WCNSFVScalarTransportPhysicsBase.

Definition at line 49 of file WCNSFVScalarTransportPhysics.h.

49 {}

◆ addSolverVariables()

void WCNSFVScalarTransportPhysics::addSolverVariables ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 35 of file WCNSFVScalarTransportPhysics.C.

 {
   // For compatibility with Modules/NavierStokesFV syntax
   if (!_has_scalar_equation)
     return;
 
   auto params = getFactory().getValidParams("INSFVScalarFieldVariable");
   assignBlocks(params, _blocks);
   params.set<MooseEnum>("face_interp_method") =
       getParam<MooseEnum>("passive_scalar_face_interpolation");
   params.set<bool>("two_term_boundary_expansion") =
       getParam<bool>("passive_scalar_two_term_bc_expansion");
 
   for (const auto name_i : index_range(_passive_scalar_names))
   {
     // Dont add if the user already defined the variable
     if (!shouldCreateVariable(_passive_scalar_names[name_i], _blocks, /*error if aux*/ true))
     {
       reportPotentiallyMissedParameters({"system_names", "passive_scalar_scaling"},
                                         "INSFVScalarFieldVariable");
       continue;
     }
 
     params.set<SolverSystemName>("solver_sys") = getSolverSystem(name_i);
     if (isParamValid("passive_scalar_scaling"))
       params.set<std::vector<Real>>("scaling") = {
           getParam<std::vector<Real>>("passive_scalar_scaling")[name_i]};
 
     getProblem().addVariable("INSFVScalarFieldVariable", _passive_scalar_names[name_i], params);
   }
 }

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

◆ getAdvectedScalarNames()

const std::vector<NonlinearVariableName>& WCNSFVScalarTransportPhysicsBase::getAdvectedScalarNames ( ) const

inlineinherited

Get the names of the advected scalar quantity variables.

Definition at line 33 of file WCNSFVScalarTransportPhysicsBase.h.

Referenced by WCNSFVTurbulencePhysics::addScalarAdvectionTurbulenceKernels().

   {
     return _passive_scalar_names;
   }

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

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

◆ hasScalarEquations()

bool WCNSFVScalarTransportPhysicsBase::hasScalarEquations ( ) const

inlineinherited

Whether the physics is actually creating the scalar advection equations.

Definition at line 39 of file WCNSFVScalarTransportPhysicsBase.h.

Referenced by WCNSFVTurbulencePhysics::retrieveCoupledPhysics().

39 { return _has_scalar_equation; }

WCNSFVScalarTransportPhysicsBase::_has_scalar_equation

const bool _has_scalar_equation

A boolean to help compatibility with the old Modules/NavierStokesFV syntax or to deliberately skip ad...

Definition: WCNSFVScalarTransportPhysicsBase.h:50

◆ setSlipVelocityParams()

virtual void WCNSFVScalarTransportPhysicsBase::setSlipVelocityParams ( InputParameters & ) const

inlineprotectedvirtualinherited

Reimplemented in WCNSFVTwoPhaseMixturePhysics, and WCNSLinearFVTwoPhaseMixturePhysics.

Definition at line 44 of file WCNSFVScalarTransportPhysicsBase.h.

Referenced by WCNSLinearFVScalarTransportPhysics::addScalarAdvectionKernels(), and addScalarAdvectionKernels().

44 {}

◆ 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 WCNSFVScalarTransportPhysics::validParams ( )

static

Definition at line 19 of file WCNSFVScalarTransportPhysics.C.

Referenced by WCNSFVTwoPhaseMixturePhysics::validParams().

 {
   InputParameters params = WCNSFVScalarTransportPhysicsBase::validParams();
   params.addClassDescription("Define the Navier Stokes weakly-compressible scalar field transport "
                              "equation(s) using the nonlinear finite volume discretization");
   params.transferParam<MooseEnum>(NSFVBase::validParams(), "passive_scalar_face_interpolation");
   params.addParamNamesToGroup("passive_scalar_face_interpolation", "Numerical scheme");
   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(), WCNSFVTurbulencePhysics::addAuxiliaryVariables(), WCNSFVTurbulencePhysics::addInitialConditions(), WCNSFVFluidHeatTransferPhysicsBase::addInitialConditions(), WCNSFVScalarTransportPhysicsBase::addInitialConditions(), WCNSFVFlowPhysicsBase::addInitialConditions(), WCNSFVFluidHeatTransferPhysics::addSolverVariables(), WCNSLinearFVFluidHeatTransferPhysics::addSolverVariables(), WCNSFVFlowPhysics::addSolverVariables(), WCNSLinearFVFlowPhysics::addSolverVariables(), WCNSFVTurbulencePhysics::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(), WCNSFVTurbulencePhysics::addFlowTurbulenceKernels(), WCNSFVTurbulencePhysics::addFluidEnergyTurbulenceKernels(), WCNSFVTurbulencePhysics::addMaterials(), WCNSFVFluidHeatTransferPhysicsBase::addMaterials(), 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(), WCNSFVTurbulencePhysics::addMaterials(), and WCNSFVCoupledAdvectionPhysicsHelper::dynamicViscosityName().

◆ _flow_equations_physics

const WCNSFVFlowPhysicsBase* WCNSFVCoupledAdvectionPhysicsHelper::_flow_equations_physics

protectedinherited

Flow physics.

Definition at line 43 of file WCNSFVCoupledAdvectionPhysicsHelper.h.

◆ _has_scalar_equation

const bool WCNSFVScalarTransportPhysicsBase::_has_scalar_equation

protectedinherited

A boolean to help compatibility with the old Modules/NavierStokesFV syntax or to deliberately skip adding the equations (for example for mixtures with a stationary phase)

Definition at line 50 of file WCNSFVScalarTransportPhysicsBase.h.

Referenced by WCNSFVScalarTransportPhysicsBase::addFVBCs(), WCNSFVScalarTransportPhysicsBase::addFVKernels(), WCNSFVScalarTransportPhysicsBase::addInitialConditions(), WCNSLinearFVScalarTransportPhysics::addSolverVariables(), addSolverVariables(), WCNSFVScalarTransportPhysicsBase::hasScalarEquations(), and WCNSFVScalarTransportPhysicsBase::WCNSFVScalarTransportPhysicsBase().