HeatStructureInterface.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 "HeatStructureInterface.h"
#include "HeatConductionModel.h"
#include "GeometricalComponent.h"
 
InputParameters
HeatStructureInterface::validParams()
{
  InputParameters params = emptyInputParameters();
 
  params.addParam<FunctionName>("initial_T", "Initial temperature [K]");
  params.addParam<Real>(
      "scaling_factor_temperature", 1.0, "Scaling factor for solid temperature variable.");
 
  params.addPrivateParam<std::string>("component_type", "heat_struct");
 
  return params;
}
 
HeatStructureInterface::HeatStructureInterface(GeometricalComponent * geometrical_component)
  : _geometrical_component_hsi(*geometrical_component)
{
}
 
std::shared_ptr<HeatConductionModel>
HeatStructureInterface::buildModel()
{
  auto & factory = _geometrical_component_hsi.getMooseApp().getFactory();
 
  const std::string class_name = "HeatConductionModel";
  InputParameters params = factory.getValidParams(class_name);
  params.set<THMProblem *>("_thm_problem") = &_geometrical_component_hsi.getTHMProblem();
  params.set<HeatStructureInterface *>("_hs") = this;
  params.applyParameters(_geometrical_component_hsi.parameters());
  return factory.create<HeatConductionModel>(
      class_name, _geometrical_component_hsi.name(), params, 0);
}
 
void
HeatStructureInterface::init()
{
  _hc_model = buildModel();
}
 
void
HeatStructureInterface::check() const
{
  auto & moose_app = _geometrical_component_hsi.getMooseApp();
  bool ics_set = _geometrical_component_hsi.getTHMProblem().hasInitialConditionsFromFile() ||
                 _geometrical_component_hsi.isParamValid("initial_T");
  if (!ics_set && !moose_app.isRestarting())
    _geometrical_component_hsi.logError("Missing initial condition for temperature.");
}
 
void
HeatStructureInterface::addVariables()
{
  _hc_model->addVariables();
  if (_geometrical_component_hsi.isParamValid("initial_T"))
    _hc_model->addInitialConditions();
}
 
void
HeatStructureInterface::addMooseObjects()
{
  if (useCylindricalTransformation())
    _hc_model->addHeatEquationRZ();
  else
    _hc_model->addHeatEquationXYZ();
}
 
FunctionName
HeatStructureInterface::getInitialT() const
{
  if (_geometrical_component_hsi.isParamValid("initial_T"))
    return _geometrical_component_hsi.getParam<FunctionName>("initial_T");
  else
    _geometrical_component_hsi.mooseError(
        _geometrical_component_hsi.name(),
        ": The parameter 'initial_T' was requested but not supplied");
}