HeatTransferBase.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 "HeatTransferBase.h"
#include "FlowModelSinglePhase.h"
#include "FlowChannelBase.h"
#include "ClosuresBase.h"
#include "MooseUtils.h"

InputParameters
HeatTransferBase::validParams()
{
  InputParameters params = ConnectorBase::validParams();
  params.addRequiredParam<std::string>("flow_channel",
                                       "Name of flow channel component to connect to");
  params.addParam<bool>(
      "P_hf_transferred", false, "Is heat flux perimeter transferred from an external source?");
  params.addParam<FunctionName>("P_hf", "Heat flux perimeter [m]");
  return params;
}

HeatTransferBase::HeatTransferBase(const InputParameters & parameters)
  : ConnectorBase(parameters),
    _flow_channel_name(getParam<std::string>("flow_channel")),
    _P_hf_transferred(getParam<bool>("P_hf_transferred")),
    _P_hf_provided(isParamValid("P_hf"))
{
  addDependency(_flow_channel_name);
}

void
HeatTransferBase::init()
{
  ConnectorBase::init();

  checkComponentOfTypeExistsByName<FlowChannelBase>(_flow_channel_name);

  if (hasComponentByName<FlowChannelBase>(_flow_channel_name))
  {
    const FlowChannelBase & flow_channel = getComponentByName<FlowChannelBase>(_flow_channel_name);

    // add the name of this heat transfer component to list for flow channel
    flow_channel.addHeatTransferName(name());

    // get various data from flow channel
    _flow_channel_subdomains = flow_channel.getSubdomainNames();
    _model_type = flow_channel.getFlowModelID();
    _fp_name = flow_channel.getFluidPropertiesName();
    _A_fn_name = flow_channel.getAreaFunctionName();
    _closures_objects = flow_channel.getClosuresObjects();
  }
}

void
HeatTransferBase::initSecondary()
{
  ConnectorBase::initSecondary();

  // determine names of heat transfer variables
  if (hasComponentByName<FlowChannelBase>(_flow_channel_name))
  {
    const FlowChannelBase & flow_channel = getComponentByName<FlowChannelBase>(_flow_channel_name);

    const std::string suffix = flow_channel.getHeatTransferNamesSuffix(name());

    _P_hf_name = FlowModel::HEAT_FLUX_PERIMETER + suffix;
    _T_wall_name = FlowModel::TEMPERATURE_WALL + suffix;
    _T_wall_mat_name = FlowModel::TEMPERATURE_WALL + suffix;
    _q_wall_name = FlowModel::HEAT_FLUX_WALL + suffix;
  }
}

void
HeatTransferBase::check() const
{
  ConnectorBase::check();
}

void
HeatTransferBase::addVariables()
{
  // heat flux perimeter variable
  if (!_P_hf_transferred)
    addHeatedPerimeter();
}

void
HeatTransferBase::addMooseObjects()
{
  // create heat flux perimeter aux if not transferred from external app
  if (!_P_hf_transferred)
  {
    const std::string class_name = "FunctionAux";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<AuxVariableName>("variable") = {_P_hf_name};
    params.set<std::vector<SubdomainName>>("block") = _flow_channel_subdomains;
    params.set<FunctionName>("function") = _P_hf_fn_name;

    ExecFlagEnum execute_on(MooseUtils::getDefaultExecFlagEnum());
    execute_on = {EXEC_TIMESTEP_BEGIN, EXEC_INITIAL};
    params.set<ExecFlagEnum>("execute_on") = execute_on;

    getTHMProblem().addAuxKernel(class_name, genName(name(), "P_hf_auxkernel"), params);
  }
}

void
HeatTransferBase::addHeatedPerimeter()
{
  getTHMProblem().addSimVariable(
      false, _P_hf_name, getTHMProblem().getFlowFEType(), _flow_channel_subdomains);

  // create heat flux perimeter variable if not transferred from external app
  if (!_P_hf_transferred)
  {
    if (_P_hf_provided)
    {
      _P_hf_fn_name = getParam<FunctionName>("P_hf");
    }
    // create heat flux perimeter function if not provided; assume circular flow channel
    else
    {
      _P_hf_fn_name = genName(name(), "P_hf_fn");

      const std::string class_name = "GeneralizedCircumference";
      InputParameters params = _factory.getValidParams(class_name);
      params.set<FunctionName>("area_function") = _A_fn_name;
      getTHMProblem().addFunction(class_name, _P_hf_fn_name, params);
    }

    if (!_app.isRestarting())
      getTHMProblem().addFunctionIC(_P_hf_name, _P_hf_fn_name, _flow_channel_subdomains);
  }
}

const VariableName &
HeatTransferBase::getHeatedPerimeterName() const
{
  checkSetupStatus(INITIALIZED_SECONDARY);

  return _P_hf_name;
}

const VariableName &
HeatTransferBase::getWallTemperatureName() const
{
  checkSetupStatus(INITIALIZED_SECONDARY);

  return _T_wall_name;
}

const MaterialPropertyName &
HeatTransferBase::getWallTemperatureMatName() const
{
  checkSetupStatus(INITIALIZED_SECONDARY);

  return _T_wall_mat_name;
}

const MaterialPropertyName &
HeatTransferBase::getWallHeatFluxName() const
{
  return _q_wall_name;
}

const UserObjectName &
HeatTransferBase::getFluidPropertiesName() const
{
  checkSetupStatus(INITIALIZED_PRIMARY);

  return _fp_name;
}