Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "TemperatureWall3EqnMaterial.h"
      11             : 
      12             : registerMooseObject("ThermalHydraulicsApp", TemperatureWall3EqnMaterial);
      13             : 
      14             : InputParameters
      15           0 : TemperatureWall3EqnMaterial::validParams()
      16             : {
      17           0 :   InputParameters params = Material::validParams();
      18           0 :   params.addClassDescription("Computes the wall temperature from the fluid temperature, "
      19             :                              "the heat flux and the heat transfer coefficient");
      20           0 :   params.addRequiredParam<MaterialPropertyName>("T", "Fluid temperature");
      21           0 :   params.addRequiredParam<MaterialPropertyName>("Hw", "Heat transfer coefficient");
      22           0 :   params.addRequiredParam<MaterialPropertyName>("q_wall", "Wall heat flux");
      23             : 
      24           0 :   return params;
      25           0 : }
      26             : 
      27           0 : TemperatureWall3EqnMaterial::TemperatureWall3EqnMaterial(const InputParameters & parameters)
      28             :   : Material(parameters),
      29           0 :     _T_wall(declareProperty<Real>("T_wall")),
      30           0 :     _q_wall(getMaterialProperty<Real>("q_wall")),
      31           0 :     _Hw(getMaterialProperty<Real>("Hw")),
      32           0 :     _T(getMaterialProperty<Real>("T"))
      33             : {
      34           0 : }
      35             : 
      36             : void
      37           0 : TemperatureWall3EqnMaterial::computeQpProperties()
      38             : {
      39           0 :   if (_q_wall[_qp] == 0)
      40           0 :     _T_wall[_qp] = _T[_qp];
      41             :   else
      42             :   {
      43             :     mooseAssert(_Hw[_qp] != 0,
      44             :                 "The wall heat transfer coefficient is zero, yet the wall heat flux is nonzero.");
      45           0 :     _T_wall[_qp] = _q_wall[_qp] / _Hw[_qp] + _T[_qp];
      46             :   }
      47           0 : }