FlowChannel1PhaseBase.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 "FlowChannel1PhaseBase.h"
#include "HeatTransfer1PhaseBase.h"
#include "SlopeReconstruction1DInterface.h"
#include "THMNames.h"

InputParameters
FlowChannel1PhaseBase::validParams()
{
  InputParameters params = FlowChannelBase::validParams();

  params.addParam<FunctionName>("initial_p", "Initial pressure in the flow channel [Pa]");
  params.addParam<FunctionName>("initial_vel", "Initial velocity in the flow channel [m/s]");
  params.addParam<FunctionName>("initial_T", "Initial temperature in the flow channel [K]");
  params.addParam<FunctionName>("D_h", "Hydraulic diameter [m]");
  params.addParam<bool>("enable_heat_conduction", false, "Enable heat conduction?");
  params.addParam<MooseEnum>(
      "rdg_slope_reconstruction",
      SlopeReconstruction1DInterface<true>::getSlopeReconstructionMooseEnum("None"),
      "Slope reconstruction type for rDG spatial discretization");

  params.addParam<Real>("p_ref", 101.325e3, "Reference pressure [Pa]");
  params.addParam<Real>("T_ref", 273.15, "Reference temperature [K]");
  params.addParam<Real>("vel_ref", 1.0, "Reference velocity [m/s]");

  params.declareControllable("initial_p initial_T initial_vel");
  params.addParamNamesToGroup("initial_p initial_T initial_vel", "Variable initialization");
  params.addParamNamesToGroup("rdg_slope_reconstruction", "Numerical scheme");

  return params;
}

FlowChannel1PhaseBase::FlowChannel1PhaseBase(const InputParameters & params)
  : FlowChannelBase(params),

    _numerical_flux_name(genName(name(), "numerical_flux")),
    _rdg_slope_reconstruction(getParam<MooseEnum>("rdg_slope_reconstruction"))
{
}

void
FlowChannel1PhaseBase::init()
{
  FlowChannelBase::init();
}

std::shared_ptr<FlowModel>
FlowChannel1PhaseBase::buildFlowModel()
{
  const std::string class_name = flowModelClassName();
  InputParameters pars = _factory.getValidParams(class_name);
  pars.set<THMProblem *>("_thm_problem") = &getTHMProblem();
  pars.set<FlowChannelBase *>("_flow_channel") = this;
  pars.set<UserObjectName>("numerical_flux") = _numerical_flux_name;
  pars.set<bool>("output_vector_velocity") = getTHMProblem().getVectorValuedVelocity();
  pars.applyParameters(parameters());
  return _factory.create<FlowModel>(class_name, name(), pars, 0);
}

void
FlowChannel1PhaseBase::check() const
{
  FlowChannelBase::check();

  checkFluidProperties();

  // only 1-phase flow compatible heat transfers are allowed
  for (unsigned int i = 0; i < _heat_transfer_names.size(); i++)
  {
    if (!hasComponentByName<HeatTransfer1PhaseBase>(_heat_transfer_names[i]))
      logError("Coupled heat source '",
               _heat_transfer_names[i],
               "' is not compatible with single phase flow channel. Use single phase heat transfer "
               "component instead.");
  }

  bool ics_set = true;
  for (const auto & ic_param : ICParameters())
    ics_set = ics_set && isParamValid(ic_param);
  ics_set = ics_set || getTHMProblem().hasInitialConditionsFromFile();

  if (!ics_set && !_app.isRestarting())
  {
    // create a list of the missing IC parameters
    std::ostringstream oss;
    for (const auto & ic_param : ICParameters())
      if (!isParamValid(ic_param))
        oss << " " << ic_param;

    logError("The following initial condition parameters are missing:", oss.str());
  }
}

void
FlowChannel1PhaseBase::addMooseObjects()
{
  FlowChannelBase::addMooseObjects();

  if (!_pipe_pars_transferred)
    addHydraulicDiameterMaterial();
}

void
FlowChannel1PhaseBase::addHydraulicDiameterMaterial()
{
  const std::string mat_name = genName(name(), "D_h_material");

  if (isParamValid("D_h"))
  {
    const FunctionName & D_h_fn_name = getParam<FunctionName>("D_h");

    const std::string class_name = "ADGenericFunctionMaterial";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<std::vector<SubdomainName>>("block") = getSubdomainNames();
    params.set<std::vector<std::string>>("prop_names") = {THM::HYDRAULIC_DIAMETER};
    params.set<std::vector<FunctionName>>("prop_values") = {D_h_fn_name};
    getTHMProblem().addMaterial(class_name, mat_name, params);
  }
  else
  {
    const std::string class_name = "ADHydraulicDiameterCircularMaterial";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<std::vector<SubdomainName>>("block") = getSubdomainNames();
    params.set<MaterialPropertyName>("D_h_name") = THM::HYDRAULIC_DIAMETER;
    params.set<std::vector<VariableName>>("A") = {THM::AREA};
    getTHMProblem().addMaterial(class_name, mat_name, params);
  }
}

void
FlowChannel1PhaseBase::getHeatTransferVariableNames()
{
  FlowChannelBase::getHeatTransferVariableNames();

  for (unsigned int i = 0; i < _n_heat_transfer_connections; i++)
  {
    const HeatTransfer1PhaseBase & heat_transfer =
        getComponentByName<HeatTransfer1PhaseBase>(_heat_transfer_names[i]);

    _Hw_1phase_names.push_back(heat_transfer.getWallHeatTransferCoefficient1PhaseName());
  }
}