InletMassFlowRateTemperature1Phase.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "InletMassFlowRateTemperature1Phase.h"
#include "FlowModelSinglePhase.h"

registerMooseObject("ThermalHydraulicsApp", InletMassFlowRateTemperature1Phase);

InputParameters
InletMassFlowRateTemperature1Phase::validParams()
{
  InputParameters params = FlowBoundary1Phase::validParams();
  params.addRequiredParam<Real>("m_dot", "Prescribed mass flow rate [kg/s]");
  params.addRequiredParam<Real>("T", "Prescribed temperature [K]");
  params.addParam<std::vector<FunctionName>>(
      "passives", {}, "Prescribed passive transport functions [amount/m^3]");
  params.addParam<bool>("reversible", true, "True for reversible, false for pure inlet");
  params.declareControllable("m_dot T");
  params.addClassDescription("Boundary condition with prescribed mass flow rate and temperature "
                             "for 1-phase flow channels.");
  return params;
}

InletMassFlowRateTemperature1Phase::InletMassFlowRateTemperature1Phase(
    const InputParameters & params)
  : FlowBoundary1Phase(params)
{
}

void
InletMassFlowRateTemperature1Phase::check() const
{
  FlowBoundary1Phase::check();

  auto fm = dynamic_cast<const FlowModelSinglePhase *>(_flow_model.get());
  if (fm)
  {
    const auto n_passives = getParam<std::vector<FunctionName>>("passives").size();
    if (n_passives != fm->passiveTransportSolutionVariableNames().size())
      logError("The number of entries in 'passives' must match the number of passive transport "
               "variables provided in the connected flow channel.");
  }
  else
    logError("Incompatible flow model. Make sure you use this component with single phase flow "
             "channel.");
}

void
InletMassFlowRateTemperature1Phase::addMooseObjects()
{
  ExecFlagEnum userobject_execute_on(MooseUtils::getDefaultExecFlagEnum());
  userobject_execute_on = {EXEC_INITIAL, EXEC_LINEAR, EXEC_NONLINEAR};

  // boundary flux user object
  {
    const std::string class_name = "ADBoundaryFlux3EqnGhostMassFlowRateTemperature";
    InputParameters params = _factory.getValidParams(class_name);
    params.applyParameters(parameters());
    params.set<Real>("mass_flow_rate") = getParam<Real>("m_dot");
    params.set<Real>("normal") = _normal;
    params.set<UserObjectName>("numerical_flux") = _numerical_flux_name;
    params.set<UserObjectName>("fluid_properties") = _fp_name;
    params.set<ExecFlagEnum>("execute_on") = userobject_execute_on;
    getTHMProblem().addUserObject(class_name, _boundary_uo_name, params);
    connectObject(params, _boundary_uo_name, "m_dot", "mass_flow_rate");
    connectObject(params, _boundary_uo_name, "T");
  }

  // BCs
  addWeakBCs();
}