Adds objects to perform hydrodynamic coupling at injection and production points. More...

#include <CoupledInjectionProductionPhysics.h>

Inheritance diagram for CoupledInjectionProductionPhysics:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	CoupledInjectionProductionPhysics (const InputParameters &parameters)

virtual InputParameters	getAdditionalRMParams () const

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

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

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

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

unsigned int	dimension () const

const ActionComponent &	getActionComponent (const ComponentName &comp_name) const

void	checkComponentType (const ActionComponent &component) const

virtual void	addComponent (const ActionComponent &component)

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

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

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

void	timedAct ()

MooseObjectName	uniqueActionName () const

const std::string &	specificTaskName () const

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

void	appendTask (const std::string &task)

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

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

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

PerfGraph &	perfGraph ()

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

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

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

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::string	unique_action_name_param

static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Member Functions
void	addPPSourceDiracKernel (const Point &point, const NonlinearVariableName &var, const PostprocessorName &pp_name)
	Adds a PorousFlowPointSourceFromPostprocessor.

void	addReceiverPostprocessor (const PostprocessorName &pp_name)
	Adds a Receiver.

void	addPointValuePostprocessor (const VariableName &var, const Point &point, const PostprocessorName &pp_name)
	Adds a PointValue.

void	addPostprocessorTransfer (const PostprocessorName &pp_name, bool from_multi_app)
	Adds a MultiAppPostprocessorTransfer.

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

bool	isTransient () const

Factory &	getFactory ()

Factory &	getFactory () const

virtual FEProblemBase &	getProblem ()

virtual const FEProblemBase &	getProblem () const

void	prepareCopyVariablesFromMesh () const

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

std::string	prefix () const

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

void	saveSolverVariableName (const VariableName &var_name)

void	saveAuxVariableName (const VariableName &var_name)

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

bool	solverVariableExists (const VariableName &var_name) const

const SolverSystemName &	getSolverSystem (unsigned int variable_index) const

const SolverSystemName &	getSolverSystem (const VariableName &variable_name) const

void	addRequiredPhysicsTask (const std::string &task)

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

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

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

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

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

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

bool	isVariableFV (const VariableName &var_name) const

bool	isVariableScalar (const VariableName &var_name) const

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

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

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

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

virtual void	checkIntegrity () const

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

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

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

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

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

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

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

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

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

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

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

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

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Protected Attributes
const std::vector< Point > &	_injection_points
	Injection points.

const std::vector< Point > &	_production_points
	Production points.

std::vector< Point >	_points
	Injection and production points.

std::vector< std::string >	_labels
	Label for each point.

std::vector< SolverSystemName >	_system_names

std::vector< unsigned int >	_system_numbers

const bool	_verbose

const MooseEnum &	_preconditioning

std::vector< SubdomainName >	_blocks

std::string	_registered_identifier

std::string	_specific_task_name

std::set< std::string >	_all_tasks

ActionWarehouse &	_awh

const std::string &	_current_task

std::shared_ptr< MooseMesh > &	_mesh

std::shared_ptr< MooseMesh > &	_displaced_mesh

std::shared_ptr< FEProblemBase > &	_problem

PerfID	_act_timer

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

MooseApp &	_pg_moose_app

const std::string	_prefix

const Parallel::Communicator &	_communicator

Private Member Functions
void	addDiracKernels () override

void	addPostprocessors () override

void	addTransfers () override

virtual void	addRelationshipManagers (Moose::RelationshipManagerType input_rm_type) override

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

virtual void	addRelationshipManagers (Moose::RelationshipManagerType when_type)

void	initializePhysics ()

virtual void	initializePhysicsAdditional ()

virtual void	checkIntegrityEarly () const

virtual void	addSolverVariables ()

virtual void	addAuxiliaryVariables ()

virtual void	addInitialConditions ()

virtual void	addFEKernels ()

virtual void	addFVKernels ()

virtual void	addFVInterpolationMethods ()

virtual void	addNodalKernels ()

virtual void	addDGKernels ()

virtual void	addScalarKernels ()

virtual void	addInterfaceKernels ()

virtual void	addFVInterfaceKernels ()

virtual void	addFEBCs ()

virtual void	addFVBCs ()

virtual void	addNodalBCs ()

virtual void	addPeriodicBCs ()

virtual void	addFunctions ()

virtual void	addAuxiliaryKernels ()

virtual void	addMaterials ()

virtual void	addFunctorMaterials ()

virtual void	addUserObjects ()

virtual void	addCorrectors ()

virtual void	addMultiApps ()

virtual void	addVectorPostprocessors ()

virtual void	addReporters ()

virtual void	addOutputs ()

virtual void	addPreconditioning ()

virtual void	addExecutioner ()

virtual void	addExecutors ()

void	checkRequiredTasks () const

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

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

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

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

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

const InputParameters &	forwardParameters () const

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

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

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

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

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

const std::string &	forwardType () const

virtual const std::string &	forwardName () const

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

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

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

Private Attributes
MooseEnum	_is_transient

std::vector< VariableName >	_solver_var_names

std::vector< VariableName >	_aux_var_names

unsigned int	_dim

std::set< std::string >	_required_tasks

const ParallelParamObject &	_parent

MooseApp &	_meta_data_app

const MooseObject *const	_meta_data_object

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

const PhysicsBase *const	_customer_class

Detailed Description

Adds objects to perform hydrodynamic coupling at injection and production points.

Definition at line 17 of file CoupledInjectionProductionPhysics.h.

Constructor & Destructor Documentation

◆ CoupledInjectionProductionPhysics()

CoupledInjectionProductionPhysics::CoupledInjectionProductionPhysics ( const InputParameters & parameters )

Definition at line 32 of file CoupledInjectionProductionPhysics.C.

  : PhysicsBase(parameters),
    _injection_points(getParam<std::vector<Point>>("injection_points")),
    _production_points(getParam<std::vector<Point>>("production_points"))
{
  _points = _injection_points;
  _points.insert(_points.end(), _production_points.begin(), _production_points.end());
 
  for (const auto i : index_range(_injection_points))
    _labels.push_back("inj" + std::to_string(i + 1));
  for (const auto i : index_range(_production_points))
    _labels.push_back("pro" + std::to_string(i + 1));
}

Member Function Documentation

◆ addDiracKernels()

void CoupledInjectionProductionPhysics::addDiracKernels ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 48 of file CoupledInjectionProductionPhysics.C.

{
  for (const auto i : index_range(_points))
  {
    addPPSourceDiracKernel(_points[i], "porepressure", "mass_rate_" + _labels[i]);
    addPPSourceDiracKernel(_points[i], "temperature", "energy_rate_" + _labels[i]);
  }
}

◆ addPointValuePostprocessor()

void CoupledInjectionProductionPhysics::addPointValuePostprocessor	(	const VariableName &	var,
		const Point &	point,
		const PostprocessorName &	pp_name
	)

protected

Adds a PointValue.

Definition at line 104 of file CoupledInjectionProductionPhysics.C.

{
  const std::string class_name = "PointValue";
  auto params = _factory.getValidParams(class_name);
  params.set<VariableName>("variable") = var;
  params.set<Point>("point") = point;
  params.set<ExecFlagEnum>("execute_on") = {EXEC_INITIAL, EXEC_TIMESTEP_END};
  getProblem().addPostprocessor(class_name, pp_name, params);
}

Referenced by addPostprocessors().

◆ addPostprocessors()

void CoupledInjectionProductionPhysics::addPostprocessors ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 58 of file CoupledInjectionProductionPhysics.C.

{
  for (const auto i : index_range(_points))
  {
    addPointValuePostprocessor("porepressure", _points[i], "p_" + _labels[i]);
    addPointValuePostprocessor("temperature", _points[i], "T_" + _labels[i]);
    addReceiverPostprocessor("mass_rate_" + _labels[i]);
    addReceiverPostprocessor("energy_rate_" + _labels[i]);
  }
}

◆ addPostprocessorTransfer()

void CoupledInjectionProductionPhysics::addPostprocessorTransfer	(	const PostprocessorName &	pp_name,
		bool	from_multi_app
	)

protected

Adds a MultiAppPostprocessorTransfer.

Definition at line 117 of file CoupledInjectionProductionPhysics.C.

{
  const std::string class_name = "MultiAppPostprocessorTransfer";
  auto params = _factory.getValidParams(class_name);
  if (from_multi_app)
    params.set<MultiAppName>("from_multi_app") = getParam<MultiAppName>("multi_app");
  else
    params.set<MultiAppName>("to_multi_app") = getParam<MultiAppName>("multi_app");
  params.set<PostprocessorName>("from_postprocessor") = pp_name;
  params.set<PostprocessorName>("to_postprocessor") = pp_name;
  if (from_multi_app)
    params.set<MooseEnum>("reduction_type") = "average";
  params.set<ExecFlagEnum>("execute_on") = {EXEC_INITIAL, EXEC_TIMESTEP_END};
  getProblem().addTransfer(class_name, pp_name + "_transfer", params);
}

Referenced by addTransfers().

◆ addPPSourceDiracKernel()

void CoupledInjectionProductionPhysics::addPPSourceDiracKernel	(	const Point &	point,
		const NonlinearVariableName &	var,
		const PostprocessorName &	pp_name
	)

protected

Adds a PorousFlowPointSourceFromPostprocessor.

Definition at line 83 of file CoupledInjectionProductionPhysics.C.

{
  const std::string class_name = "PorousFlowPointSourceFromPostprocessor";
  auto params = _factory.getValidParams(class_name);
  params.set<NonlinearVariableName>("variable") = var;
  params.set<Point>("point") = point;
  params.set<PostprocessorName>("mass_flux") = pp_name;
  getProblem().addDiracKernel(class_name, "dirac_" + pp_name, params);
}

Referenced by addDiracKernels().

◆ addReceiverPostprocessor()

void CoupledInjectionProductionPhysics::addReceiverPostprocessor ( const PostprocessorName & pp_name )

protected

Adds a Receiver.

Definition at line 96 of file CoupledInjectionProductionPhysics.C.

{
  const std::string class_name = "Receiver";
  auto params = _factory.getValidParams(class_name);
  getProblem().addPostprocessor(class_name, pp_name, params);
}

Referenced by addPostprocessors().

◆ addTransfers()

void CoupledInjectionProductionPhysics::addTransfers ( )

overrideprivatevirtual

Reimplemented from PhysicsBase.

Definition at line 70 of file CoupledInjectionProductionPhysics.C.

{
  if (isParamValid("multi_app"))
    for (const auto i : index_range(_points))
    {
      addPostprocessorTransfer("p_" + _labels[i], false);
      addPostprocessorTransfer("T_" + _labels[i], false);
      addPostprocessorTransfer("mass_rate_" + _labels[i], true);
      addPostprocessorTransfer("energy_rate_" + _labels[i], true);
    }
}

◆ validParams()

InputParameters CoupledInjectionProductionPhysics::validParams ( )

static

Definition at line 18 of file CoupledInjectionProductionPhysics.C.

{
  InputParameters params = PhysicsBase::validParams();
 
  params.addRequiredParam<std::vector<Point>>("injection_points", "List of injection points [m]");
  params.addRequiredParam<std::vector<Point>>("production_points", "List of production points [m]");
  params.addParam<MultiAppName>("multi_app", "MultiApp to transfer to and from");
 
  params.addClassDescription(
      "Adds objects to perform hydrodynamic coupling at injection and production points.");
 
  return params;
}