SimpleTurbine1Phase.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 "SimpleTurbine1Phase.h"
#include "FlowModelSinglePhase.h"

registerMooseObject("ThermalHydraulicsApp", SimpleTurbine1Phase);

InputParameters
SimpleTurbine1Phase::validParams()
{
  InputParameters params = JunctionParallelChannels1Phase::validParams();

  params.addRequiredParam<Real>("power", "Turbine power [W]");
  params.addRequiredParam<bool>("on", "Flag determining if turbine is operating or not [-]");

  params.declareControllable("power on");

  params.addClassDescription(
      "Simple turbine model that extracts prescribed power from the working fluid");

  return params;
}

SimpleTurbine1Phase::SimpleTurbine1Phase(const InputParameters & params)
  : JunctionParallelChannels1Phase(params),
    _on(getParam<bool>("on")),
    _power(getParam<Real>("power")),
    _W_dot_var_name(junctionVariableName("W_dot"))
{
}

void
SimpleTurbine1Phase::addVariables()
{
  JunctionParallelChannels1Phase::addVariables();

  addJunctionVariable(false, _W_dot_var_name);
}

void
SimpleTurbine1Phase::buildVolumeJunctionUserObject()
{
  ExecFlagEnum execute_on(MooseUtils::getDefaultExecFlagEnum());
  execute_on = {EXEC_INITIAL, EXEC_LINEAR, EXEC_NONLINEAR};

  {
    const std::string class_name = "ADSimpleTurbine1PhaseUserObject";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<bool>("use_scalar_variables") = false;
    params.set<subdomain_id_type>("junction_subdomain_id") = _junction_subdomain_id;
    params.set<std::vector<BoundaryName>>("boundary") = _boundary_names;
    params.set<std::vector<Real>>("normals") = _normals;
    params.set<std::vector<processor_id_type>>("processor_ids") = getConnectedProcessorIDs();
    params.set<std::vector<UserObjectName>>("numerical_flux_names") = _numerical_flux_names;
    params.set<Real>("volume") = _volume;
    params.set<std::string>("component_name") = name();
    params.set<std::vector<VariableName>>("A") = {FlowModel::AREA};
    params.set<std::vector<VariableName>>("rhoA") = {FlowModelSinglePhase::RHOA};
    params.set<std::vector<VariableName>>("rhouA") = {FlowModelSinglePhase::RHOUA};
    params.set<std::vector<VariableName>>("rhoEA") = {FlowModelSinglePhase::RHOEA};
    params.set<std::vector<VariableName>>("rhoV") = {_rhoV_var_name};
    params.set<std::vector<VariableName>>("rhouV") = {_rhouV_var_name};
    params.set<std::vector<VariableName>>("rhovV") = {_rhovV_var_name};
    params.set<std::vector<VariableName>>("rhowV") = {_rhowV_var_name};
    params.set<std::vector<VariableName>>("rhoEV") = {_rhoEV_var_name};
    params.set<RealVectorValue>("dir_c0") = _directions[0];
    params.set<Real>("K") = _K;
    params.set<Real>("A_ref") = _A_ref;
    params.set<bool>("on") = _on;
    params.set<Real>("W_dot") = _power;
    params.set<UserObjectName>("fp") = _fp_name;
    params.set<ExecFlagEnum>("execute_on") = execute_on;
    getTHMProblem().addUserObject(class_name, _junction_uo_name, params);
    connectObject(params, _junction_uo_name, "power", "W_dot");
    connectObject(params, _junction_uo_name, "on");
    connectObject(params, _junction_uo_name, "K");
  }
}

void
SimpleTurbine1Phase::addMooseObjects()
{
  JunctionParallelChannels1Phase::addMooseObjects();

  {
    const std::string nm = genName(name(), "W_dot_aux");
    const std::string class_name = "SimpleTurbinePowerFieldAux";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<AuxVariableName>("variable") = _W_dot_var_name;
    params.set<Real>("value") = _power;
    params.set<bool>("on") = _on;
    params.set<std::vector<SubdomainName>>("block") = getSubdomainNames();
    getTHMProblem().addAuxKernel(class_name, nm, params);
    connectObject(params, nm, "power", "value");
    connectObject(params, nm, "on");
  }
}