Component1DConnection.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 "Component1DConnection.h"
#include "Component1D.h"

const std::map<std::string, Component1DConnection::EEndType>
    Component1DConnection::_end_type_to_enum{{"IN", IN}, {"OUT", OUT}};

template <>
Component1DConnection::EEndType
THM::stringToEnum(const std::string & s)
{
  return stringToEnum<Component1DConnection::EEndType>(s, Component1DConnection::_end_type_to_enum);
}

InputParameters
Component1DConnection::validParams()
{
  InputParameters params = Component::validParams();
  return params;
}

Component1DConnection::Component1DConnection(const InputParameters & params) : Component(params) {}

void
Component1DConnection::setupMesh()
{
  for (const auto & connection : _connections)
  {
    const std::string & comp_name = connection._component_name;

    if (hasComponentByName<Component1D>(comp_name))
    {
      _boundary_names.push_back(connection._boundary_name);

      const Component1D & comp = getComponentByName<Component1D>(comp_name);
      for (auto && conn : comp.getConnections(connection._end_type))
      {
        // get info from the connection
        _positions.push_back(conn._position);
        _elems.push_back(conn._elem);
        _sides.push_back(conn._side);
        _nodes.push_back(conn._node->id());
        _normals.push_back(conn._normal);
        _boundary_ids.push_back(conn._boundary_id);
        _directions.push_back(comp.getDirection());
      }
    }
    else
      logError("Trying to connect to a component '",
               comp_name,
               "', but there is no such component in the simulation. Please check your spelling.");
  }
}

void
Component1DConnection::init()
{
  Component::init();

  if (_connections.size() > 0)
  {
    std::vector<UserObjectName> fp_names;
    std::vector<THM::FlowModelID> flow_model_ids;
    for (const auto & connection : _connections)
    {
      const std::string comp_name = connection._component_name;
      if (hasComponentByName<Component1D>(comp_name))
      {
        const Component1D & comp = getTHMProblem().getComponentByName<Component1D>(comp_name);

        // add to list of subdomain names
        const std::vector<SubdomainName> & subdomain_names = comp.getSubdomainNames();
        _connected_subdomain_names.insert(
            _connected_subdomain_names.end(), subdomain_names.begin(), subdomain_names.end());
      }
    }
  }
  else
    logError("The component is not connected.");
}

void
Component1DConnection::check() const
{
  Component::check();

  for (const auto & comp_name : _connected_component_names)
    checkComponentOfTypeExistsByName<Component1D>(comp_name);
}

void
Component1DConnection::addConnection(const BoundaryName & boundary_name)
{
  const size_t oparenthesis_pos = boundary_name.find('(');
  if (oparenthesis_pos != std::string::npos)
  {
    logError("You are using the old connection format 'comp_name(end)'. Please update your input "
             "file to the new one 'comp_name:end'.");
  }
  else
  {
    const size_t colon_pos = boundary_name.rfind(':');
    // if it has a colon, assume 'component_name:end_type' format
    if (colon_pos != std::string::npos)
    {
      const std::string connected_component_name = boundary_name.substr(0, colon_pos);
      const std::string str_end =
          boundary_name.substr(colon_pos + 1, boundary_name.length() - colon_pos - 1);
      const EEndType end_type = THM::stringToEnum<EEndType>(str_end);

      _connections.push_back(Connection(boundary_name, connected_component_name, end_type));
      _connected_component_names.push_back(connected_component_name);

      // Add dependency because the connected component's setupMesh() must be called
      // before this component's setupMesh().
      addDependency(connected_component_name);
    }
    else
    {
      logError("Incorrect connection specified '",
               boundary_name,
               "'. Valid connection format is 'component_name:end_type'.");
    }
  }
}

void
Component1DConnection::checkNumberOfConnections(const unsigned int & n_connections) const
{
  if (_connections.size() != n_connections)
    logError("The number of connections (", _connections.size(), ") must equal ", n_connections);
}

const std::vector<dof_id_type> &
Component1DConnection::getNodeIDs() const
{
  checkSetupStatus(MESH_PREPARED);

  return _nodes;
}

const std::vector<BoundaryName> &
Component1DConnection::getBoundaryNames() const
{
  checkSetupStatus(MESH_PREPARED);

  return _boundary_names;
}